Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,621,033 A→G 8.3% K200R (AAG→AGG)  tctD → transcriptional regulator TctD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,621,0330AG 8.3% 135.3 / 2.8 48K200R (AAG→AGG) tctDtranscriptional regulator TctD
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base G (2/2);  total (24/24)
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.85e-01

CTCCCGCGAGCAGAGTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCTGGTGGTGCTGATGCTG  >  NC_002947/1620887‑1621170
                                                                                                                                                  |                                                                                                                                         
cTCCCGCGAGCAGAGTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGc                                                                                                                                         >  6:303717/1‑149 (MQ=255)
       gAGCAGAGTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAGGCTCGACgg                                                                                                                                  >  4:22776/1‑149 (MQ=255)
       gAGCAGAGTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAGGCTCGACgg                                                                                                                                  >  4:72299/1‑149 (MQ=255)
          cAGAGTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAGGCTCGACGGGCa                                                                                                                               <  3:72299/149‑1 (MQ=255)
          cAGAGTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAGGCTCGACGGGCa                                                                                                                               <  3:22776/149‑1 (MQ=255)
           agagTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCAc                                                                                                                              <  1:365788/149‑1 (MQ=255)
             agTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGAc                                                                                                                                    >  8:147019/1‑141 (MQ=255)
             agTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGAc                                                                                                                                    <  7:147019/141‑1 (MQ=255)
                    tGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGc                                                                                                                      >  2:316959/1‑148 (MQ=255)
                    tGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGc                                                                                                                      >  4:241274/1‑148 (MQ=255)
                     ggAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCtgtg                                                                                                                        >  6:61559/1‑145 (MQ=255)
                     ggAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCtgtg                                                                                                                        <  5:61559/145‑1 (MQ=255)
                      gAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCAt                                                                                                                   >  5:62165/1‑149 (MQ=255)
                      gAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCAt                                                                                                                   >  6:68250/1‑149 (MQ=255)
                             tGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACt                                                                                                            >  7:277938/1‑149 (MQ=255)
                                    tcgcCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCcgcg                                                                                                     >  5:46890/1‑149 (MQ=255)
                                          gggccgggTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTgg                                                                                               <  3:240307/149‑1 (MQ=255)
                                           ggccgggTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTggg                                                                                              <  5:68250/149‑1 (MQ=255)
                                                    atgaGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGc                                                                                      <  5:303717/148‑1 (MQ=255)
                                                     tgaGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCtt                                                                                    <  8:277938/149‑1 (MQ=255)
                                                     tgaGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCt                                                                                     <  6:62165/148‑1 (MQ=255)
                                                        gCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCtt                                                                                    >  6:262959/1‑146 (MQ=255)
                                                        gCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGa                                                                                  <  5:262959/148‑1 (MQ=255)
                                                                cagcTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGa                                                                         >  3:364688/1‑149 (MQ=255)
                                                                      gctgCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGt                                                                    <  1:316959/148‑1 (MQ=255)
                                                                      gctgCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGt                                                                    >  1:62852/1‑148 (MQ=255)
                                                                       ctgCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGa                                                                  <  6:96896/149‑1 (MQ=255)
                                                                        tgCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGa                                                                  >  2:306971/1‑148 (MQ=255)
                                                                           ccAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCAcc                                                                             <  4:283871/134‑1 (MQ=255)
                                                                           ccAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCAcc                                                                             >  3:283871/1‑134 (MQ=255)
                                                                            cAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACg                                                             <  3:241274/149‑1 (MQ=255)
                                                                               gtgtTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCa                                                          <  1:306971/149‑1 (MQ=255)
                                                                                      gCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTgcg                                                   >  4:296124/1‑149 (MQ=255)
                                                                                                    gccagccCTGACGCCATCGAAATCTACATTCACCGGCTGCGCaaaaa                                                                                                                                           <  6:19854/47‑1 (MQ=255)
                                                                                                    gccagccCTGACGCCATCGAAATCTACATTCACCGGCTGCGCaaaaa                                                                                                                                           >  5:19854/1‑47 (MQ=255)
                                                                                                     ccagccCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTggg                                    <  2:62852/149‑1 (MQ=255)
                                                                                                     ccagccCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCACCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTggg                                    <  7:146502/149‑1 (MQ=255)
                                                                                                      cagccCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTggg                                    >  6:164217/1‑148 (MQ=255)
                                                                                                          ccTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGATATTTGCTTGAGCACCGCGATGCGTGACAACGGAAGCCTGCGCGGGCGCTTGCTGGGTAAtt                               >  3:320/1‑149 (MQ=255)
                                                                                                              aCGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTg                             <  3:296124/147‑1 (MQ=255)
                                                                                                                gCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGCCAGCCTGCGCCGGCGCTTGCTGGGTAATTTGgcgc                         >  2:178916/1‑149 (MQ=255)
                                                                                                                   aTCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGt                      >  8:178901/1‑149 (MQ=255)
                                                                                                                         aTCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCtgg                <  1:326747/149‑1 (MQ=255)
                                                                                                                         aTCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCtgg                <  4:364688/149‑1 (MQ=255)
                                                                                                                           cTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCtggtg              >  6:106988/1‑149 (MQ=255)
                                                                                                                            tACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCtggtgg             >  8:172983/1‑149 (MQ=255)
                                                                                                                             aCATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGACGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCTggtggt            >  2:63255/1‑149 (MQ=255)
                                                                                                                               aTTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCTGGTGGTGc          <  2:178999/149‑1 (MQ=255)
                                                                                                                                    ccGGCTGCGCAAATAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCTGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCTGGCGCTTGCTGTGTAATTTGGCGCTGTTGCTGGTGGTGCTGATg     <  2:282631/149‑1 (MQ=255)
                                                                                                                                       gCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCTGGTGGTGCTGATGCTg  <  7:178901/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
CTCCCGCGAGCAGAGTGTGCTGGAAGCGCTGATCGCTCGCCCGGGCCGGGTGATGAGCAAGGAGCAGCTGGCTGCCCAGGTGTTTGGCCTGGACGAAGAAGCCAGCCCTGACGCCATCGAAATCTACATTCACCGGCTGCGCAAAAAGCTCGACGGGCACCCTGTGGCCATCGTCACTTTCCGCGGTCTGGGCTATTTGCTTGAGCACCGCGATGCGTGACAACGGCAGCCTGCGCGGGCGCTTGCTGGGTAATTTGGCGCTGTTGCTGGTGGTGCTGATGCTG  >  NC_002947/1620887‑1621170

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