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.C88.2% 42.8 / 1.6 17coding (75/120 nt)BW25113_RS26130hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (6/9);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAA  >  NZ_CP009273/1417693‑1417862
                                                                                            |                                                                                
ttGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCGCCGCCTCCGGCGTGGATAATGCAG‑‑CCt                                                                                   <  2:640992/90‑1 (MQ=18)
      cTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTcccccc                                                                            >  2:903442/1‑90 (MQ=33)
                 cgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCGACTGGCAGac                                                                 <  1:210061/90‑1 (MQ=33)
                  ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGac                                                                 <  2:123074/90‑1 (MQ=255)
                            cAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                       <  1:371900/90‑1 (MQ=255)
                                        tcagtcagTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAAACACCGCTGAACGGGATTATTTCa                                           >  2:737395/1‑90 (MQ=11)
                                             cagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                      >  2:427390/1‑90 (MQ=255)
                                                 gCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCggg                                  >  1:403749/1‑90 (MQ=255)
                                                   gTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGaa                                <  1:352265/90‑1 (MQ=14)
                                                   gTGAAGCCACCGCCGCCTCCGGCGTTGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGaa                                <  1:831554/90‑1 (MQ=14)
                                                    tGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAAc                               <  2:65650/90‑1 (MQ=255)
                                                     gAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACg                              >  1:529143/1‑90 (MQ=255)
                                                         ccaccaCCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                          <  2:529143/90‑1 (MQ=255)
                                                              accGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  1:427390/90‑1 (MQ=255)
                                                                      cGGCGTGGATAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAgg             >  2:675947/1‑90 (MQ=255)
                                                                          gTGGATAATGCAGCCCCCCCCCCCGTCTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  2:693476/90‑1 (MQ=21)
                                                                               tAATGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTg    >  2:691861/1‑90 (MQ=255)
                                                                                 aTGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGaa  <  1:499605/90‑1 (MQ=255)
                                                                                 aTGCAG‑‑CCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGaa  <  2:972251/90‑1 (MQ=255)
                                                                                            |                                                                                
TTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAG‑‑CCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAA  >  NZ_CP009273/1417693‑1417862

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