Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,955,621 C→T 8.1% G171G (GGC→GGT elbB → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,955,6210CT 8.1% 136.3 / 2.1 49G171G (GGC→GGTelbBputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base T (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

ATTGGCCGAAGCCTTTGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACG  >  NC_002947/5955474‑5955765
                                                                                                                                                   |                                                                                                                                                
aTTGGCCGAAGCCTTTGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCa                                                                                                                                                 <  7:66230/149‑1 (MQ=255)
      cGAAGCCTTTGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCAt                                                                                                                                            >  5:70840/1‑148 (MQ=255)
      cGAAGCCTTTGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCAt                                                                                                                                            >  4:121653/1‑148 (MQ=255)
           ccTTTGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGTGCCGCCCAtcttttc                                                                                                                                      <  4:71880/149‑7 (MQ=255)
            cTTTGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAAt                                                                                                                                     >  4:253320/1‑149 (MQ=255)
             tttGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATg                                                                                                                                    >  5:312625/1‑149 (MQ=255)
             tttGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATg                                                                                                                                    >  3:152625/1‑149 (MQ=255)
                gCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCg                                                                                                                                  <  2:312139/148‑1 (MQ=255)
                 ccGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGa                                                                                                                                 <  1:342607/148‑1 (MQ=255)
                 ccGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGa                                                                                                                                 <  3:226577/148‑1 (MQ=255)
                 ccGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGa                                                                                                                                 <  2:20566/148‑1 (MQ=255)
                 ccGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGAt                                                                                                                                >  6:68754/1‑149 (MQ=255)
                 ccGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGAt                                                                                                                                <  6:312625/149‑1 (MQ=255)
                 ccGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGAt                                                                                                                                <  4:152625/149‑1 (MQ=255)
                  cgccGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCCCCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGCTg                                                                                                                               >  7:302744/4‑149 (MQ=255)
                        cTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTacac                                                                                                                          >  1:222388/1‑148 (MQ=255)
                               cccGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGt                                                                                                                  <  3:253320/149‑1 (MQ=255)
                                ccGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTc                                                                                                                 >  4:131544/1‑149 (MQ=255)
                                 ggTAGGCCTGATGTGCTTCTCGCCGGCGCTGGCGGCGAAGATTTTCGGGCCGGGTGTGGTCTGCAGCATCGGGCGGGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATTAAGAATGCGATGTACACGACATTGTCg                                                                                                                <  4:249647/148‑1 (MQ=255)
                                   tacgCCTTAGCTGCTTCTCGCCGGCGCTGTCGGCGAAGATTTACGGGCCGGGGGTTGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGaa                                                                                                              <  8:203808/146‑1 (MQ=255)
                                   tAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGaa                                                                                                              <  2:222388/149‑1 (MQ=255)
                                   tAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGaa                                                                                                              <  8:197828/149‑1 (MQ=255)
                                           aTCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTCCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACAccc                                                                                                       <  5:368385/148‑1 (MQ=255)
                                               gCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGac                                                                                                           >  7:1172/1‑140 (MQ=255)
                                               gCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGac                                                                                                           <  8:1172/140‑1 (MQ=255)
                                                  tctcGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAg                                                                                               >  5:101749/1‑149 (MQ=255)
                                                    tcGCCGGCGCTAGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCt                                                                                             >  4:232482/1‑149 (MQ=255)
                                                        cGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGTACCCATGAAGAATGCGa                                                                                                                                 <  7:208837/109‑1 (MQ=255)
                                                        cGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGTACCCATGAAGAATGCGa                                                                                                                                 <  7:203362/109‑1 (MQ=255)
                                                        cGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGTACCCATGAAGAATGCGa                                                                                                                                 >  8:203362/1‑109 (MQ=255)
                                                        cGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGTACCCATGAAGAATGCGa                                                                                                                                 >  8:208837/1‑109 (MQ=255)
                                                         ggCGCTAGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGa                                                                                        <  3:232482/149‑1 (MQ=255)
                                                                     gAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTTTCGAAAAGATGGGCTCCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCg                                                                             >  5:299019/1‑148 (MQ=255)
                                                                        gATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACa                                                                         >  7:46072/1‑149 (MQ=255)
                                                                          tttACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAATATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATg                                                                       <  6:98306/149‑1 (MQ=255)
                                                                               gggCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGc                                                                  >  8:28281/1‑149 (MQ=255)
                                                                                ggCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGTCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCTATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGcc                                                                 <  6:10747/149‑1 (MQ=255)
                                                                                ggCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGcc                                                                 <  5:68754/149‑1 (MQ=255)
                                                                                              tGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGc                                                     >  6:364977/1‑147 (MQ=255)
                                                                                                cACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACATGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAgatc                                                   >  4:17486/1‑144 (MQ=255)
                                                                                                     tCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGgcggcg                                            <  8:46072/149‑1 (MQ=255)
                                                                                                                gCCGGCACTCCGTCGGCTGGCGAAAAGCTGGGCGGCACCCTTGAAGAATGCGATGTACACGACATTGTCGAAGACACCCATCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGAt                                    <  6:299019/146‑1 (MQ=255)
                                                                                                                       cTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCtg                          <  6:70840/149‑1 (MQ=255)
                                                                                                                          cggcggCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCtgg                         <  3:121653/147‑1 (MQ=255)
                                                                                                                           ggcggcTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCtggtgg                      <  7:28281/149‑1 (MQ=255)
                                                                                                                                   tCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTg              >  3:340451/1‑149 (MQ=255)
                                                                                                                                     gAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCt            >  2:12565/1‑149 (MQ=255)
                                                                                                                                       aaaGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTg           >  3:75627/1‑148 (MQ=255)
                                                                                                                                           aTGGGCGGCACCCATGAAAAATGCGATGTACACGACATTCTCGAAGCCACCCAGCGCAAGCTCGTGAGTACGCCGGCGTAGATGGAGGCGCAGTCGATCAGCTCGGCGGCGGGCGGGCTGTACGAGCTGGGGGCCCGAGTGGTGGAGCt      >  6:278974/1‑149 (MQ=255)
                                                                                                                                               gcggcACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACg  >  7:55956/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATTGGCCGAAGCCTTTGCCGACGCCTGCAAACCCGTAGGCCTGATCTGCATCTCGCCGGCGCTGGCGGCGAAGATTTACGGGCCGGGGGTGGTCTGCACCATCGGGCGCGATGCCGGCACTGCGGCGGCTGTCGAAAAGATGGGCGGCACCCATGAAGAATGCGATGTACACGACATTGTCGAAGACACCCAGCGCAAGCTCGTGACTACCCCGGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACG  >  NC_002947/5955474‑5955765

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