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.93.3% 53.6 / ‑1.1 15coding (52/441 nt)BW25113_RS06225hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (9/5);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAA  >  NZ_CP009273/1240587‑1240745
                                                                                |                                                                              
tttttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAg                                                                      >  2:507505/1‑90 (MQ=255)
     aaTTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAgg                                                                 <  1:657349/90‑1 (MQ=255)
          ttGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGtt                                                            >  1:105371/1‑90 (MQ=255)
                cGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCt                                                      <  1:512070/90‑1 (MQ=255)
                 gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTg                                                     >  2:960079/1‑90 (MQ=255)
                      ttttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGGTGCGCtt                                                >  2:537091/1‑90 (MQ=255)
                                 ttCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACAttt                                     <  1:507505/90‑1 (MQ=255)
                                       aaataaatATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAAtt                               <  1:493778/90‑1 (MQ=255)
                                                        aGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGc               <  1:473894/90‑1 (MQ=255)
                                                              cAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGGTTTTGGTGCGCTTTTTTAAAATTTTGAATTTATTAAACCTTTTTGCAGATTaa        >  1:789975/1‑90 (MQ=255)
                                                               aGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAAAATTTAGAAATGATCAATCCTTCTGGCAGATTAAg       >  1:326148/1‑90 (MQ=255)
                                                                 tCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTgcgc                                                  >  1:336819/1‑45 (MQ=255)
                                                                 tCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTgcgc                                                  <  2:336819/45‑1 (MQ=255)
                                                                   aGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTTTTTAAAATTTAAAAATGATCAAACCTTTTGGCAGATTAAACCaa   >  2:697034/1‑90 (MQ=255)
                                                                    ggCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAAAATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCaaa  >  1:110661/1‑90 (MQ=255)
                                                                                |                                                                              
TTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAA  >  NZ_CP009273/1240587‑1240745

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