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.94.1% 60.3 / ‑1.0 17coding (52/441 nt)BW25113_RS06225hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (7/9);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

ATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAA  >  NZ_CP009273/1240582‑1240744
                                                                                     |                                                                             
aTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATgggcg                                                                           >  1:714216/1‑88 (MQ=255)
   cTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGAt                                                                       <  2:988441/90‑1 (MQ=255)
   cTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGAt                                                                       <  2:1609/90‑1 (MQ=255)
                      gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTg                                                    <  1:1182857/90‑1 (MQ=255)
                        tgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTgcg                                                  <  1:548185/90‑1 (MQ=255)
                           ttttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTTGCt                                                     >  2:1459171/1‑84 (MQ=255)
                           ttttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCt                                                     <  1:1459171/84‑1 (MQ=255)
                              tGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTAtt                                            <  1:451992/90‑1 (MQ=255)
                                 tgtgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                         <  2:554648/90‑1 (MQ=255)
                                   tgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGGTTTTGCTGCGCTTATTTaaa                                        >  1:1434637/1‑89 (MQ=255)
                                           aaaaTAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAAt                               <  1:1451251/90‑1 (MQ=255)
                                                    atTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGGTTTTGCTGCGCTTTTTTAAAATTTAGAATTGATTAATc                      >  1:1177363/1‑90 (MQ=255)
                                                              gTAGCCAGTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAAAATTTAGAATTGATTAATCCTTCTGGCAg            >  1:366036/1‑90 (MQ=255)
                                                                     gTCAGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTTGCAGATTAAGc     >  2:1457795/1‑90 (MQ=255)
                                                                        aGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGGTGCGCCTTTTTAAAAATTAGAATTGATTAAACCTTTTGGCAGATTTAGCCaa  >  1:1250716/1‑90 (MQ=255)
                                                                        aGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCaa  >  1:1379425/1‑90 (MQ=255)
                                                                        aGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                         <  1:97416/51‑1 (MQ=255)
                                                                        aGGCGTTAAAAATGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                         >  2:97416/1‑51 (MQ=255)
                                                                                     |                                                                             
ATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAA  >  NZ_CP009273/1240582‑1240744

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: