Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,417,788 (C)6→7 coding (69/120 nt) BW25113_RS26130 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,417,7821.C92.3% 29.7 / ‑2.3 13coding (75/120 nt)BW25113_RS26130hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (5/7);  total (5/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAA  >  NZ_CP009273/1417693‑1417862
                                                                                           |                                                                                
ttGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCt                                                                                   <  1:4680/90‑1 (MQ=33)
                  ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGac                                                                 <  2:226357/90‑1 (MQ=255)
                   gTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGaca                                                                <  1:186956/90‑1 (MQ=255)
                     tGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACAcc                                                              <  2:500103/90‑1 (MQ=255)
                       cacaTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                            <  1:305052/90‑1 (MQ=11)
                                     cTGTCAGTGAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattatt                                              >  2:530794/1‑90 (MQ=11)
                                             cagtGCGTGAAGCCACCACCGCCCCCGGCGTGGATAATGCAGCCCCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                      <  1:71237/90‑1 (MQ=11)
                                               gtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCg                                    >  2:490539/1‑90 (MQ=255)
                                                 gCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCggg                                  >  2:151441/1‑90 (MQ=255)
                                                         ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCa                                                                    <  1:497779/48‑1 (MQ=11)
                                                         ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCa                                                                    >  2:497779/1‑48 (MQ=14)
                                                                         cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  1:478835/90‑1 (MQ=255)
                                                                                 aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGaa  >  1:433563/1‑90 (MQ=255)
                                                                                 aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGaa  <  1:530794/90‑1 (MQ=255)
                                                                                           |                                                                                
TTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAA  >  NZ_CP009273/1417693‑1417862

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