Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 670,635 A→T L324Q (CTG→CAG)  holA ← DNA polymerase III subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913670,6350AT95.7% 62.0 / ‑6.0 23L324Q (CTG→CAG) holADNA polymerase III subunit delta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (11/11);  minor base G (0/1);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACC  >  NC_000913/670396‑670877
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
cATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGTTGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAACTCCGTCCGCCCGGGGTTTATGGCACAAAAGAAGAGATAACCCTTCCTGCTCCGCCcac                                                                                                                                                                                                                                         >  2:190143/1‑251 (MQ=255)
                                gCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTg                                                                                                                                                                                                           >  1:29941/1‑249 (MQ=255)
                                      gACGATGCGGAGTAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGGTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGTTTTATGGCACAACAGAAGAGATAACCCTTCCTGTTCTGCTCACACGGACTGACCGTAATCTTGTTTTTGGGTGAGTTCCGt                                                                                                                                                                                                   <  2:121106/251‑1 (MQ=255)
                                             cGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAAGTCGCCAGCGTGTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTGCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGATAGAACCCTTGCTGCTCTGCTCACACTTACGGACCGTAAGCTGGGTTGAGGGTGACTTCCGTTCGTGTc                                                                                                                                                                                            <  2:282603/251‑1 (MQ=255)
                                              ggaggcaacatattaGAGATGATTTTGAACCCCCTCAGACACATCAAATTCGACAGCCTTTCCACGGGCTTTCGACGCCCATAGTGCACCGGATCAACGCTGCCGCCCAACAGAGCCTGTACAGATTTCATCTCAACAGTCGATAAATACGTCCGCCAGGGGTTCATGGCACAACAGAAGAGATCACCCTTCCTGCTCTGCCCTAACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGGGTcc                                                                                                                                                                                           >  2:181153/11‑250 (MQ=255)
                                                                                                cGCCAGCGCTTCCACGGGATTCAGATGACCATAGTGAACCGGATCAAAGGTGCCGCAAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATACATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTACTCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGACTGACGTAACTGCGTCTGACTTAAGCGATTTAACGcc                                                                                                                                         >  2:303697/1‑251 (MQ=255)
                                                                                                  ccAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCTGAGCAAAGGTGCCTCCAAACAGAGCCTGTAAAGTTTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCGGCTCGGGCCACTCTGAGTTACCGTTATCTTGTGTGAGGTTGAGTGCGGTTCGTGTCTGGATTTGTACTGCCGGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTc                                                                                                                                       <  1:285772/251‑1 (MQ=255)
                                                                                                              cGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCGCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGcc                                                                                                                           <  2:212245/251‑1 (MQ=255)
                                                                                                               gggTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGccc                                                                                                                          <  1:238839/251‑1 (MQ=255)
                                                                                                                   tttAGATGGCCAGAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAATATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCGCTGCCCACGCTGACTGACCGTAATCTGGTTTGAGGGTGAGTTCCGTTCGTGGCAGGAGTTGCACGGCCTGACGTAACTGCGTCGGACGTAAGCGATTTAACGCCTCGCCCATCATGCCCcgg                                                                                                                       <  2:269485/250‑1 (MQ=255)
                                                                                                                             cATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCa                                                                                                            <  1:106471/251‑1 (MQ=255)
                                                                                                                                       cgtaTCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACGTACTGACCGTAATCTTGTTTGAGGGTGATTTCCTTTCGTGTCAGGTGTTGCACGTCGGGACGTAACTGTTTCTGACTTAAGCGTTTTTACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGc                                                                                                  <  1:196259/248‑1 (MQ=255)
                                                                                                                                               aGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCaa                                                                                           >  2:143298/1‑250 (MQ=255)
                                                                                                                                                          gaCAGAGCCTGTAAAGATTTCATATCAATAGTCGATAAATACGTCCGCAGGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATATTGTTTGAGGGAGAGTTCCGTTCGTGTCAGGAGCTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATTCCCCGGCGCTTCTGTCATCCCCGATACTTAtacacacacgaccgcc                                                                               >  2:178341/2‑236 (MQ=255)
                                                                                                                                                            ccgagCCTGTAAAGAGTTCATAGCAACCGTCGATAAATACGTCCGCCAGGGGCTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGGCTGACCGTAATCTTGTTTGAGGGTGAGTGCCGTTCGTGGCAGGAGTTGCACGGCCTGACGTAAGCGCGTCTGACTTAAGCGACTTAACGCGTCGCCCAACATGCCCCGGCGGTTCTGCCGTACCCGATGCTTATCAAACAGCGCACGCAGt                                                                             <  2:106624/249‑1 (MQ=255)
                                                                                                                                                                      aaaGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCCTACCCGATGCTTATCAAACAACCCACGCAGTGGCGTATggg                                                                   >  2:168667/1‑251 (MQ=255)
                                                                                                                                                                                   cAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGACGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTTTCAAACAACGCACGCACTGGCGTATGGGCAGACTGGCGttt                                                      >  2:148221/1‑251 (MQ=255)
                                                                                                                                                                                          cGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTa                                               <  1:143298/251‑1 (MQ=255)
                                                                                                                                                                                           gATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTaa                                              <  1:281783/251‑1 (MQ=255)
                                                                                                                                                                                                              gggttATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCACACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCTGTTCGTGTCCGGCGTTGCACGGCCTGACGTAACTGAGTCTGACTTAAGCGAATTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGGCATACCCGATGCTTATAAAACAACGAACGAAGGGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGAAAAAACAggtaacg                           >  1:144736/4‑245 (MQ=255)
                                                                                                                                                                                                                        acaacaGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAAtgtg                 <  2:62075/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaCCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGACCCGGCGGGTCTGCCATACCCGCTGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAAAAACAGTTCCCGCTGTAATGTGCGTAAACAAATAAc   >  2:66467/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aCCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAAcc  >  1:20939/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
CATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACC  >  NC_000913/670396‑670877

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