Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,247,856 C→T 19.4% A146V (GCC→GTC)  PP_1978 → TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,247,8560CT19.4% 44.4 / 3.9 22A146V (GCC→GTC) PP_1978TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (7/10);  new base T (2/2);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTGCCCGCTGCCTAACGCGTTGGGCGAAATGCAGCGCCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCATGGCGATGACTACCCTG  >  NC_002947/2247723‑2247986
                                                                                                                                     |                                                                                                                                  
gCTTGCCCGCTGCCTAACGCGTTGGGCGAAATGCAGCGCCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCg                                                                                                                                  <  2:491234/136‑1 (MQ=255)
   tGCCCGCTGCCTAACGCGTTGGGCGAAATGCAGCGCCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGAcc                                                                                                                               >  2:407526/1‑136 (MQ=255)
        gCTGCCTAACGCGTTGGGCGAAATGCAGCGCCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGcc                                                                                                                          >  1:510179/1‑136 (MQ=255)
                         gCGAAATGCAGCGCCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGgct                                                                                                         >  2:441211/1‑136 (MQ=255)
                                gCAGCGCCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGgctggcgctg                                                                                                  >  1:115086/1‑136 (MQ=255)
                                    cgcCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGccc                                                                                              <  1:304031/136‑1 (MQ=255)
                                                   ttCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCa                                                                               <  1:407526/136‑1 (MQ=255)
                                                          tggTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCt                                                                        <  1:348140/135‑1 (MQ=255)
                                                          cggTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCt                                                                        <  1:275146/135‑1 (MQ=255)
                                                             tGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTc                                                                     <  1:441211/136‑1 (MQ=255)
                                                                          caTCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGCACGACCTGGCCCTGATGGTTGTGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTg                                                        >  2:164599/1‑136 (MQ=255)
                                                                                        cTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCt                                                         <  1:147625/121‑1 (MQ=255)
                                                                                        cTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCt                                                         >  2:147625/1‑121 (MQ=255)
                                                                                          gCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATcgg                                        <  2:510179/136‑1 (MQ=255)
                                                                                                      gACCGCCCTGAGGATGCCGACAAAGCGCTGGTCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAAc                            >  2:238481/1‑136 (MQ=255)
                                                                                                      gACCGCCCTGAGGATGCCGACAAAGCGCTGGTCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAAc                            >  2:241570/1‑136 (MQ=255)
                                                                                                                aGGATGCCGACAAAGCGCTGGTCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCAt                  <  1:241570/136‑1 (MQ=255)
                                                                                                                aGGATGCCGACAAAGCGCTGGTCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCAt                  <  1:238481/136‑1 (MQ=255)
                                                                                                                        gACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCAt                  <  1:55409/128‑1 (MQ=255)
                                                                                                                        gACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCAt                  >  2:55409/1‑128 (MQ=255)
                                                                                                                        gACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCATGGCGATGa          >  1:386103/1‑136 (MQ=255)
                                                                                                                                gcTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCATGGCGATGACTACCCTg  <  1:463679/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GCTTGCCCGCTGCCTAACGCGTTGGGCGAAATGCAGCGCCTGCCGGAGGTCTTCCGCGAGGTGCTGGCCGAGCACATCGAGTTGATGACTGCCGCGATGGTCGACCGCCCTGAGGATGCCGACAAAGCGCTGGCCGACCTGGCCCTGATGGTTGGGGGGCTGGCGCTGGCCCGGGCCTTGGGGACCAGTGACCTGTCTGACCGTATTCTGCGCGCCGCCAAATCGGCGGTGCTTTAACCCTGGAGCATGGCGATGACTACCCTG  >  NC_002947/2247723‑2247986

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