Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,138,686 A→T 14.5% intergenic (‑34/+284) USA300HOU_2020 ← / ← xtmA hypothetical membrane protein/terminase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,138,6860AT14.5% 73.8 / 4.7 26intergenic (‑34/+284)USA300HOU_2020/xtmAhypothetical membrane protein/terminase small subunit
Reads supporting (aligned to +/- strand):  ref base A (15/7);  new base T (2/2);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATCCGAAAGCAATAAATTCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTTGAAAAAAATATCATAC  >  CP000730/2138569‑2138818
                                                                                                                     |                                                                                                                                    
atcatcCGAAAGCAATAAATTCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCTATTGTTTACAAAACACGAAAttt                                                                                                               >  2:230990/1‑141 (MQ=255)
atcatcCGAAAGCAATAAATTCAAACTAATGCATTGTTTATGACCCATCGAATTTATTGACCAAATAGCTAGAGAAACAAACAACCCAAAATTGAAGATAAATGATATAGTAATAGCTATGGTTTACATAACACGAAAttt                                                                                                               >  2:230991/1‑141 (MQ=255)
  catcCGAAAGCAATAAATTCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAAttttt                                                                                                             >  1:78486/1‑141 (MQ=255)
  catcCGAAAGCAATAAATTCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAAttttt                                                                                                             >  1:144628/1‑141 (MQ=255)
              aTAAATTCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTtttattta                                                                                                 >  2:102954/1‑141 (MQ=255)
              aTAAATTCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTtttattta                                                                                                 >  1:69985/1‑141 (MQ=255)
                    tCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTAtattat                                                                                           >  1:208998/1‑141 (MQ=255)
                          tAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCAttt                                                                                     >  2:85547/1‑141 (MQ=255)
                                         taCCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCtt                                                                      <  2:227235/141‑1 (MQ=255)
                                                   aTTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCAtttt                                                            >  2:166195/1‑141 (MQ=255)
                                                   aTTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCAtttt                                                            >  2:132978/1‑141 (MQ=255)
                                                    tttATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTtta                                                           >  1:23864/1‑141 (MQ=255)
                                                               aaTAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTAc                                                >  2:230151/1‑141 (MQ=255)
                                                                aTAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACa                                               >  2:21333/1‑141 (MQ=255)
                                                                  aGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACaat                                             <  1:195686/141‑1 (MQ=255)
                                                                        agaAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTTTTTTTTTTTTTTATTATATATTACAATaatcaa                                       >  2:146682/1‑141 (MQ=255)
                                                                          aaaTAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATaatcaatc                                     <  1:110059/141‑1 (MQ=255)
                                                                            aTAAACACCCCAACATTTAAGATAAGTGATATAGTAATAGCTATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTg                                   <  1:230991/141‑1 (MQ=255)
                                                                            aTAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCTATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTg                                   <  1:230990/141‑1 (MQ=255)
                                                                                   acCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGtt                            <  1:241715/141‑1 (MQ=255)
                                                                                    cccAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTg                           <  1:85547/141‑1 (MQ=255)
                                                                                                   aaatGATATAGTAATAGCAATTGTTTACAAAACACGAAATGTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTTGaaaaa            >  2:227859/1‑141 (MQ=255)
                                                                                                       gATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATTATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTTGAAAAAAata        >  2:37849/1‑141 (MQ=255)
                                                                                                       gATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTTGAAAAAAata        >  1:57494/1‑141 (MQ=255)
                                                                                                         tataGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTTGAAAAAAATATc      <  1:59050/141‑1 (MQ=255)
                                                                                                             gTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGAAAAAAAAAAAAAAAAATATCATAc  <  2:95027/141‑1 (MQ=255)
                                                                                                                     |                                                                                                                                    
ATCATCCGAAAGCAATAAATTCAAACTAATGCATTGTTTATTACCCATCGAATTTATTGACCAAATAGCTAGAGAAATAAACAACCCAAAATTTAAAATAAATGATATAGTAATAGCAATTGTTTACAAAACACGAAATTTTTCATTTTTATTTATATTATCCATTTTTCTCCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTTGAAAAAAATATCATAC  >  CP000730/2138569‑2138818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑