Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,612 T→A 12.5% L129Q (CTG→CAG)  dnaEB → error‑prone DNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,528,6120TA12.5% 155.7 / 13.1 56L129Q (CTG→CAG) dnaEBerror‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base T (30/19);  new base A (3/4);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCT  >  NC_002947/3528467‑3528760
                                                                                                                                                 |                                                                                                                                                    
gACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGCACCTACGAGCGCGGATCACCCGTGCCCGGCGGCGGGCAGCGAAGGGGGCGTACCAGCTGTTTCGTGCTGACTTGCTGCTGCACCACCAAGGCCTGCTGGcc                                                                                                                                                  >  2:152016/1‑150 (MQ=255)
    ggCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTtgt                                                                                                                                              <  1:778477/150‑1 (MQ=255)
     gCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTtgtg                                                                                                                                             <  2:812997/150‑1 (MQ=255)
      cccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCA‑CAAGGCCAGCAGGCCTtggtgg                                                                                                                                           <  1:142092/150‑5 (MQ=255)
       ccAATCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtggg                                                                                                                                           <  1:1078434/150‑1 (MQ=255)
        cAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgggt                                                                                                                                          >  2:314182/1‑150 (MQ=255)
        cAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgggt                                                                                                                                          <  1:426146/150‑1 (MQ=255)
                       ggTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGGTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTgg                                                                                                                           <  2:536466/150‑1 (MQ=255)
                       ggTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTgg                                                                                                                           <  1:663519/150‑1 (MQ=255)
                        gTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGc                                                                                                                          <  2:873052/150‑1 (MQ=255)
                        gTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtgg                                                                                                                                           <  1:571943/133‑5 (MQ=255)
                        gTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGAGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtgg                                                                                                                                           >  2:571943/1‑129 (MQ=255)
                           gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGAt                                                                                                                       >  2:1057075/1‑150 (MQ=255)
                                 cTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCt                                                                                                                 >  2:1129563/1‑150 (MQ=255)
                                  tGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGcc                                                                                                                                                  <  2:127183/116‑1 (MQ=255)
                                  tGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGcc                                                                                                                                                  >  1:127183/1‑116 (MQ=255)
                                  tGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTa                                                                                                                >  1:302208/1‑150 (MQ=255)
                                      gggTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCgg                                                                                                            >  2:548987/1‑150 (MQ=255)
                                       ggTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGgc                                                                                                           >  1:703080/1‑150 (MQ=255)
                                        gTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGc                                                                                                                          >  1:896128/1‑134 (MQ=255)
                                        gTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGc                                                                                                                          <  2:896128/134‑1 (MQ=255)
                                        gTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGCGGCGCTACTCCTACCGgcc                                                                                                          >  2:189207/1‑149 (MQ=255)
                                         ttACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGgcgc                                                                                                         >  1:477581/1‑150 (MQ=255)
                                         ttACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg                                                                                                                                            <  2:931481/115‑4 (MQ=255)
                                         ttACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg                                                                                                                                            >  1:931481/1‑112 (MQ=255)
                                           aCCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGt                                                                                                       <  1:1053328/150‑1 (MQ=255)
                                             cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg                                                                                                     >  2:887001/1‑150 (MQ=255)
                                               gAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCt                                                                                                   >  1:655174/1‑150 (MQ=255)
                                                       gcgcgcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGCGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAgtg                                                                                           <  2:655174/150‑1 (MQ=255)
                                                         gcgcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAgtgtg                                                                                         >  1:594133/1‑150 (MQ=255)
                                                               aTCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCt                                                                                   >  2:867197/1‑150 (MQ=255)
                                                                 cACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCt                                                                                   >  2:955396/1‑148 (MQ=255)
                                                                 cACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCt                                                                                   <  1:955396/148‑1 (MQ=255)
                                                                    ccGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAgcg                                                                              >  2:714161/1‑150 (MQ=255)
                                                                           cggcggcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCCCCACCAAGGCCTGCTGGGCTTGTGGGTGGCCGCTGACAGTTGCGATACTGCTACCGGCCCGTGGGCGCGCCGTGGGTTCGCTGAGCCCCTGtgg                                                                       >  2:66660/1‑150 (MQ=255)
                                                                                gcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGtggctggc                                                                  >  2:803539/1‑150 (MQ=255)
                                                                                      agagAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg                                                                                                                                            >  1:207335/1‑67 (MQ=38)
                                                                                      agagAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg                                                                                                                                            <  2:207335/70‑4 (MQ=38)
                                                                                          aaGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGc                                                                                                                                                   >  1:619599/1‑59 (MQ=255)
                                                                                          aaGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGc                                                                                                                                                   <  2:619599/59‑1 (MQ=255)
                                                                                             ggggCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCAt                                                     >  2:763664/1‑150 (MQ=255)
                                                                                              gggCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATc                                                    <  1:1057075/150‑1 (MQ=255)
                                                                                                  gTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCgggg                                                >  2:297323/1‑150 (MQ=255)
                                                                                                      cAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGt                                            >  1:683288/1‑150 (MQ=255)
                                                                                                       aGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAgc                                                                               <  2:632707/114‑1 (MQ=255)
                                                                                                       aGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAgc                                                                               >  1:632707/1‑114 (MQ=255)
                                                                                                       aGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTg                                           >  1:755493/1‑150 (MQ=255)
                                                                                                               cGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGcc                                   <  2:683288/150‑1 (MQ=255)
                                                                                                                    tgaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTgcg                              >  2:463614/1‑150 (MQ=255)
                                                                                                                       cTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCt                           <  1:459015/150‑1 (MQ=255)
                                                                                                                        ttGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTg                          >  1:91625/1‑150 (MQ=255)
                                                                                                                        ttGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTg                          >  1:264266/1‑150 (MQ=255)
                                                                                                                                  accaccAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGCGCAGTGACGATGCCGTGCGCATGCAACtgctgc                >  2:474807/1‑150 (MQ=255)
                                                                                                                                     accaAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCggg             <  2:477581/150‑1 (MQ=255)
                                                                                                                                         aGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCg                                                   >  2:547513/1‑108 (MQ=255)
                                                                                                                                         aGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCg                                                   <  1:547513/108‑1 (MQ=255)
                                                                                                                                               gctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggnggc   <  2:703080/150‑1 (MQ=255)
                                                                                                                                                ctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCt  >  2:284693/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCT  >  NC_002947/3528467‑3528760

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