Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,609 T→A 12.5% L128Q (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,6090TA12.5% 152.0 / 13.6 56L128Q (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 = 6.00e-01

TTGTGCGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGG  >  NC_002947/3528461‑3528749
                                                                                                                                                    |                                                                                                                                            
ttGTGCGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGTAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCctg                                                                                                                                             <  2:1047753/150‑1 (MQ=255)
      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                                                                                                               <  2:896128/134‑1 (MQ=255)
                                              gTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGc                                                                                                               >  1:896128/1‑134 (MQ=255)
                                              gTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGCGGCGCTACTCCTACCGgcc                                                                                               >  2:189207/1‑149 (MQ=255)
                                               ttACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGgcgc                                                                                              >  1:477581/1‑150 (MQ=255)
                                               ttACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg                                                                                                                                 >  1:931481/1‑112 (MQ=255)
                                               ttACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg                                                                                                                                 <  2:931481/115‑4 (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                                                                                                                                 <  2:207335/70‑4 (MQ=38)
                                                                                            agagAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg                                                                                                                                 >  1:207335/1‑67 (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                                        <  1:547513/108‑1 (MQ=255)
                                                                                                                                               aGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCg                                        >  2:547513/1‑108 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TTGTGCGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGG  >  NC_002947/3528461‑3528749

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: ‑