Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,998,666 C→T 9.5% E127K (GAA→AAA)  PP_3527 ← TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,998,6660CT 9.5% 112.9 / 3.1 42E127K (GAA→AAA) PP_3527TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (13/25);  new base T (2/2);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TGTGCAGGTCGATGACCTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCGGATATCCCCGTACAAC  >  NC_002947/3998519‑3998805
                                                                                                                                                   |                                                                                                                                           
tgtgCAGGTCGATGACCTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCg                                                                                                                                            <  7:58257/149‑1 (MQ=255)
tgtgCAGGTCGATGACCTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCg                                                                                                                                            <  3:64756/149‑1 (MQ=255)
    cAGGTCGATGACCTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTg                                                                                                                                         >  4:259826/1‑148 (MQ=255)
     aGGTCGATGACCTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTTGGTg                                                                                                                                         >  6:11023/1‑147 (MQ=255)
                cTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCa                                                                                                                            <  1:106850/149‑1 (MQ=255)
                cTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCa                                                                                                                            >  8:162039/1‑149 (MQ=255)
                  cGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCAc                                                                                                                           >  8:109138/1‑148 (MQ=255)
                   gATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCAc                                                                                                                           <  6:64123/147‑1 (MQ=255)
                              ggAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTg                                                                                                              >  3:217011/1‑149 (MQ=255)
                                    cGCCTTCCTGAACGCCACGGCGCAGGGTGGTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCa                                                                                                        <  1:21239/149‑1 (MQ=255)
                                      ccTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATg                                                                                                      <  7:109138/149‑1 (MQ=255)
                                           cTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGa                                                                                                  <  3:259826/148‑1 (MQ=255)
                                              aaCGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTTGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCg                                                                                               >  3:157674/1‑148 (MQ=255)
                                                cGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTTGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAg                                                                                            <  5:11023/149‑1 (MQ=255)
                                                cGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAg                                                                                            <  2:145212/149‑1 (MQ=255)
                                                   cACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGAc                                                                                                                  <  8:129228/124‑1 (MQ=255)
                                                   cACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGAc                                                                                                                  >  7:129228/1‑124 (MQ=255)
                                                           aGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTccac                                                                                 <  5:270939/149‑1 (MQ=255)
                                                            gggTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATgg                                                                                                     >  6:132805/1‑128 (MQ=255)
                                                            gggTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATgg                                                                                                     <  5:132805/128‑1 (MQ=255)
                                                            gggTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTccacc                                                                                <  5:38352/149‑1 (MQ=255)
                                                              gtgtTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAg                                                                              <  2:152708/149‑1 (MQ=255)
                                                              gtgtTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAg                                                                              >  3:151515/1‑149 (MQ=255)
                                                                    cccAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGAACGAAGGTGAACTCCACCAGgcggcg                                                                        >  7:63669/1‑149 (MQ=255)
                                                                              ccAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGc                                                              >  1:267579/1‑149 (MQ=255)
                                                                               cAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTg                                                               <  2:152244/147‑1 (MQ=255)
                                                                                                   cgcgAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACtccag                                          <  4:151515/148‑1 (MQ=255)
                                                                                                       aaTGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACtccaggtcca                                     <  4:217011/149‑1 (MQ=255)
                                                                                                        aTGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACtccaggtcca                                     <  5:195167/148‑1 (MQ=255)
                                                                                                                      cTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGc                                                                  >  5:44479/1‑105 (MQ=255)
                                                                                                                      cTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGc                                                                  <  6:44479/105‑1 (MQ=255)
                                                                                                                      cTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGAACGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCt                      <  8:63669/149‑1 (MQ=255)
                                                                                                                           cGTACTCGCCATAGTGGATGTTTTTGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCg                  <  4:157674/148‑1 (MQ=255)
                                                                                                                             tACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCGGat               <  7:162039/149‑1 (MQ=255)
                                                                                                                                 cGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTc                                                            <  6:249462/100‑1 (MQ=255)
                                                                                                                                 cGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTc                                                            >  5:249462/1‑100 (MQ=255)
                                                                                                                                     aTAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCGGATATCCCCg        >  2:200847/1‑148 (MQ=255)
                                                                                                                                        gTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCgg                 <  2:207995/136‑1 (MQ=255)
                                                                                                                                        gTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCgg                 >  1:207995/1‑136 (MQ=255)
                                                                                                                                          ggATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCGGATATCCCCGTACaa   <  3:33471/148‑1 (MQ=255)
                                                                                                                                           gATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCGGATATCCCCGTACAAc  <  2:267579/148‑1 (MQ=255)
                                                                                                                                                  tCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGt                                                     <  8:262467/90‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TGTGCAGGTCGATGACCTCGATGCCGGCGCGGAACACGCCTTCCTGAACGCCACGGCGCAGGGTGTTGCCCAAGGCCTCCAGCACCAGGCTGCTCATTTCGCGAATGGCCGGCGACTGCTTGACGTACTCGCCATAGTGGATGTTTTCGGTGCAGATGATGCGCACGAAATCGACGTTGCGGTCATGGTGATCGAAGGTGAACTCCACCAGGCGGCGGATCGCCTGCTCAGCCGGTAACGACTCCAGGTCCAGGCTCTGTTCGGTCTTGCGGATATCCCCGTACAAC  >  NC_002947/3998519‑3998805

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