Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,240,670 Δ1 bp coding (55/441 nt) BW25113_RS06225 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,240,6670G.78.6% 28.5 / 10.3 14coding (52/441 nt)BW25113_RS06225hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base . (6/5);  total (8/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

TCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTT  >  NZ_CP009273/1240583‑1240754
                                                                                    |                                                                                       
tCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCg                                                                                    <  2:581634/90‑1 (MQ=255)
                                tgtgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                                   <  2:474282/90‑1 (MQ=255)
                                  tgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAAc                                                  >  2:226858/1‑90 (MQ=255)
                                         gAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAAAATTTAGaa                                          >  2:454150/1‑90 (MQ=255)
                                                 atatTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCaa                                  <  2:554196/90‑1 (MQ=255)
                                                           cAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGc                        <  1:284294/90‑1 (MQ=255)
                                                             gTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTTTTTAAAATTTAGAATTGATCAATCCTTTTGGCAg                      >  2:570934/1‑90 (MQ=255)
                                                                  cAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTaa                 <  2:405269/90‑1 (MQ=255)
                                                                    gTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGc               >  2:239202/1‑90 (MQ=255)
                                                                    gTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGc               <  2:557332/90‑1 (MQ=255)
                                                                     tCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAACATTTAGAATTGATCAAACCTTCTGGCAGATTAAAcc              >  2:377910/1‑90 (MQ=255)
                                                                       aGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAAAATTTTGAATTGATCAATCCTTCTTGCAGATTAAACCaa            >  1:221623/1‑90 (MQ=255)
                                                                           gTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGGGCTGCGTTTTGAAAATGTAAAATAGGTAAAACCTTTTGGCAGATTGAGGCaaaaaa        >  1:357267/1‑88 (MQ=255)
                                                                                  aTGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGtt  >  2:204547/1‑90 (MQ=255)
                                                                                    |                                                                                       
TCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTT  >  NZ_CP009273/1240583‑1240754

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