Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,873,647 G→A 8.0% V143M (GTG→ATG)  cigR → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,873,6470GA 8.0% 139.2 / 2.0 50V143M (GTG→ATG) cigRmembrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base A (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

CATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACAC  >  NC_002947/2873508‑2873793
                                                                                                                                           |                                                                                                                                                  
cATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGAAGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCg                                                                                                                                           <  2:180185/149‑1 (MQ=255)
 tttCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCgg                                                                                                                                          >  1:91117/2‑149 (MQ=255)
  ttCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCCGGTCGCCACGGCGAccaga                                                                                                                                         <  7:50235/149‑4 (MQ=255)
     agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGat                                                                                                                                       >  3:122838/1‑148 (MQ=255)
      gaagaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCAATACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCggggaac                                                                                                                                     >  8:298562/1‑145 (MQ=255)
       aagaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCTGGatcct                                                                                                                                    >  8:293563/1‑147 (MQ=255)
       aagaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATCCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGatcat                                                                                                                                    <  7:293563/149‑1 (MQ=255)
         gaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCt                                                                                                                                  <  5:269197/149‑1 (MQ=255)
            ccTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATg                                                                                                                               <  6:74672/149‑1 (MQ=255)
             cTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTAt                                                                                                                                >  5:15293/1‑147 (MQ=255)
             cTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGa                                                                                                                              <  5:66208/149‑1 (MQ=255)
                gCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCATGGCGACCGGGATCATCTATGAg                                                                                                                             <  4:91162/147‑1 (MQ=255)
                 cccGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCTCCGTGGCGACCGGGATCATCTATGAGgtg                                                                                                                          >  2:88084/1‑149 (MQ=255)
                  ccGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgtgt                                                                                                                         >  1:139557/1‑149 (MQ=255)
                  ccGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgtgt                                                                                                                         >  8:252421/1‑149 (MQ=255)
                  ccGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTAATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgtgt                                                                                                                         >  5:242785/1‑149 (MQ=255)
                            ccccTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGAGCATCTATGAGGTGGTGACCgggg                                                                                                                >  8:142288/1‑148 (MQ=255)
                                              gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCg                                                                                                                                           >  3:41312/1‑103 (MQ=255)
                                              gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCg                                                                                                                                           <  4:41312/103‑1 (MQ=255)
                                               ccAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGTGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                         <  1:95462/136‑1 (MQ=255)
                                               ccAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                         >  2:95462/1‑136 (MQ=255)
                                                cAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGt                                                                                           >  7:276839/1‑149 (MQ=255)
                                                  agaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTAATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGc                                                                                         <  6:242785/149‑1 (MQ=255)
                                                       cTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCa                                                                                    <  1:88084/149‑1 (MQ=255)
                                                         ggATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACg                                                                                                                   <  8:250232/116‑1 (MQ=255)
                                                                  tcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGt                                                                         <  1:300977/149‑1 (MQ=255)
                                                                      ttGGGGCAGCTGCCTCACTACGACGGGTAAGAGTGGATGCAAGCGGGCGCACACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCt                                                                     >  6:241879/1‑149 (MQ=255)
                                                                        ggggCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGCCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGAAGTCCGCTCGGTCCCTg                                                                    >  7:31026/1‑148 (MQ=255)
                                                                           gCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAgg                                                                >  2:280813/1‑149 (MQ=255)
                                                                           gCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAgg                                                                >  7:269642/1‑149 (MQ=255)
                                                                              gCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGgcgc                                                             <  7:252421/149‑1 (MQ=255)
                                                                                 gCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTgg                                                          >  4:91680/1‑149 (MQ=255)
                                                                                            acgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCt                                                                     >  4:30428/1‑127 (MQ=255)
                                                                                              gggTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATg                                             <  6:77580/149‑1 (MQ=255)
                                                                                                 tATGAGTGGATGCAAGCGGTCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGc                                          <  2:63261/149‑1 (MQ=255)
                                                                                                   tGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCt                                         <  3:87454/148‑1 (MQ=255)
                                                                                                    gAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTgc                                       >  2:102942/1‑149 (MQ=255)
                                                                                                         gATGCAAGCGGGCGCAGACCTGATTCAGGGCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagc                                  <  6:17330/149‑1 (MQ=255)
                                                                                                             cAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCagcagc                               <  1:102942/148‑1 (MQ=255)
                                                                                                             cAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGcc                              <  6:15293/149‑1 (MQ=255)
                                                                                                               cgCGGGCGCAGACCTGATTCTGGGCGCCGTGGCGACCGGGATCATCTATGAGGGGTTGACCGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCa                            <  7:142288/148‑1 (MQ=255)
                                                                                                                  gggCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACg                         <  8:269642/149‑1 (MQ=255)
                                                                                                                          aCCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATc                  >  3:173444/1‑148 (MQ=255)
                                                                                                                           ccTGATGCTGGGCGCCGTGGCGACCGGTATCATCTATGAGGTGTTGAACGGGGCATTCGATAACGGCCGCAGGTGCGAGCAGTCCGCGCGGTCCGTGCAGGCGGGGGCGTGCGGGCTTTGGGCTGCGCAGCAGCCCAACGGCCTATCaa                <  2:91117/149‑1 (MQ=255)
                                                                                                                            cggATTGTGGTGGCCGTGGGCACCGGTATCTTATATTCGGTGTTGAACGGGGCTTTCCACTACTCCGGAAGGTGCTAGCTTTCGGGGCGGTCCCTGCTGGCGGGGGGGTGCTTGGTTTGTGCTGCGAATCAGCCCAACGGCCTATCAAc               <  8:31026/147‑1 (MQ=255)
                                                                                                                              gATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGc             >  4:63704/1‑149 (MQ=255)
                                                                                                                               atcccggTCGCCACGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCt            >  1:19442/6‑149 (MQ=255)
                                                                                                                               aTTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGCGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCt            <  7:298562/149‑1 (MQ=255)
                                                                                                                                 cccggTCGCCACGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGc          >  3:114257/4‑149 (MQ=255)
                                                                                                                                         ccGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAacac  >  3:206219/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACAC  >  NC_002947/2873508‑2873793

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