Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,873,639 T→C 8.0% L140P (CTG→CCG)  cigR → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,873,6390TC 8.0% 140.6 / 2.3 50L140P (CTG→CCG) cigRmembrane protein
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (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 = 9.97e-01

CATTGCCTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAA  >  NC_002947/2873495‑2873779
                                                                                                                                                |                                                                                                                                             
cATTGCCTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTc                                                                                                                                           <  4:269699/149‑1 (MQ=255)
       tCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTgg                                                                                                                                    >  3:404932/1‑149 (MQ=255)
       tCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTgg                                                                                                                                    >  7:162539/1‑149 (MQ=255)
       tCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTgg                                                                                                                                    >  1:4173/1‑149 (MQ=255)
          aGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGa                                                                                                                                 <  3:169524/149‑1 (MQ=255)
          aGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGa                                                                                                                                 >  8:307978/1‑149 (MQ=255)
                 cAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGAccagaat                                                                                                                          >  1:75020/1‑144 (MQ=255)
                  agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTgg                                                                                                                                             <  3:285229/129‑1 (MQ=255)
                  agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTgg                                                                                                                                             >  4:285229/1‑129 (MQ=255)
                  agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGat                                                                                                                          >  5:319085/1‑148 (MQ=255)
                  agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGAccagaat                                                                                                                          <  7:307978/148‑6 (MQ=255)
                       aaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGACCAGAatcaggtc                                                                                                                    <  2:75020/149‑5 (MQ=255)
                        aCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTAt                                                                                                                   <  2:4173/149‑1 (MQ=255)
                        aCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTAt                                                                                                                   <  3:311818/149‑1 (MQ=255)
                          cTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTAt                                                                                                                   >  8:388332/1‑147 (MQ=255)
                          cTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGa                                                                                                                 <  7:153623/149‑1 (MQ=255)
                               ccGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGGGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgtgt                                                                                                            >  1:335079/1‑149 (MQ=255)
                               ccGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTAt                                                                                                                   >  3:127803/1‑142 (MQ=255)
                               ccGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTAt                                                                                                                   <  4:127803/142‑1 (MQ=255)
                                  cgGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGa                                                                                                         <  3:386977/149‑1 (MQ=255)
                                     cAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACg                                                                                                      <  6:319085/149‑1 (MQ=255)
                                          cccTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑Ca                                                                                                >  2:203140/1‑149 (MQ=255)
                                           ccTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CAt                                                                                               >  4:29986/1‑149 (MQ=255)
                                               ccACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGa                                                                                           <  1:203140/149‑1 (MQ=255)
                                                  ccTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACt                                                                                         >  3:216669/1‑148 (MQ=255)
                                                           gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAg                                                                               <  4:404932/149‑1 (MQ=255)
                                                                  agCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAg                                                                        <  2:335079/149‑1 (MQ=255)
                                                                           gtcgtcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCt                                                               <  8:162539/149‑1 (MQ=255)
                                                                            tcgtcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTc                                                              >  3:46929/1‑149 (MQ=255)
                                                                              gtcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAg                                                                               >  6:54168/1‑130 (MQ=255)
                                                                              gtcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAg                                                                               <  5:54168/130‑1 (MQ=255)
                                                                               tcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGt                                                           >  4:212583/1‑149 (MQ=255)
                                                                                      gggCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGt                                                           <  6:309321/142‑1 (MQ=255)
                                                                                      gggCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGt                                                           >  5:309321/1‑142 (MQ=255)
                                                                                             tGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGa                                                                                                                 <  1:49914/82‑1 (MQ=255)
                                                                                             tGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGa                                                                                                                 >  2:49914/1‑82 (MQ=255)
                                                                                               ccTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTc                                                          >  4:247885/1‑134 (MQ=255)
                                                                                               ccTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTc                                                          <  3:247885/134‑1 (MQ=255)
                                                                                               ccTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGc                                           >  6:88965/1‑149 (MQ=255)
                                                                                                cTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGatcat                                                                                                                       >  7:377011/1‑73 (MQ=255)
                                                                                                cTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGatcat                                                                                                                       <  8:377011/73‑1 (MQ=255)
                                                                                                 tCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGt                                         >  5:180466/1‑149 (MQ=255)
                                                                                                        gacgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGc                                   >  4:137323/1‑148 (MQ=255)
                                                                                                             gTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGTTCCCTGCAGGCGCTGGCGTGCCTGCGATggg                             <  7:388332/149‑1 (MQ=255)
                                                                                                                 gAGTGGATGCAATCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGATGTGTTGAACGGGGCCATTCGACTAAGGCCGCAGGTGCGAGCAg                                                                     <  1:232690/106‑1 (MQ=255)
                                                                                                                 gAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAg                                                                     >  2:232690/1‑105 (MQ=255)
                                                                                                                       aTGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                   <  5:88965/149‑1 (MQ=255)
                                                                                                                                  cgcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCACTGGCGTGCCTGCGTTGGGCTGCTGCGCAGCCCATCgcag        <  4:216669/149‑4 (MQ=255)
                                                                                                                                   gcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCATCgcagg       <  4:46929/149‑5 (MQ=255)
                                                                                                                                        ccTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCaa  >  7:61049/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CATTGCCTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGG‑CATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAA  >  NC_002947/2873495‑2873779

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