Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,990,934 G→T 21.7% D66Y (GAC→TAC)  PP_5245 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,990,9340GT21.7% 56.2 / 4.8 23D66Y (GAC→TAC) PP_5245AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (9/9);  new base T (2/3);  total (11/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 = 5.47e-01

GATGTCACCGGTCCCATGGATGCATTCGCCATCGCCAACCGGTTCTTGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAGGCCGTGCGCCAG  >  NC_002947/5990802‑5991053
                                                                                                                                    |                                                                                                                       
gATGTCACCGGTCCCATGGATGCATTCGCCATCGCCAACCGGTTCTTGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGAcc                                                                                                                      <  2:370638/136‑1 (MQ=255)
  tGTCACCGGTCCCATGGATGCATTCGCCATCGCCAACCGGTTCTTGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCgg                                                                                                                    >  1:42826/1‑136 (MQ=255)
        cGGTCCCATGGATGCATTCGCCATCGCCAACCGGTTCTTGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTg                                                                                                              >  1:436865/1‑136 (MQ=255)
                                   caACCGGTTCTTGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTg                                                                                   >  1:275068/1‑136 (MQ=255)
                                         gTTCTTGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTg                                                                             <  1:76945/136‑1 (MQ=255)
                                              tGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGcc                                                                        <  2:436865/136‑1 (MQ=255)
                                              tGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGcc                                                                        <  1:527658/136‑1 (MQ=255)
                                                      gcaCGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGTCCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGcc                                                                <  1:471716/136‑1 (MQ=255)
                                                      ccaCGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGcc                                                                <  1:517309/135‑1 (MQ=255)
                                                       caCGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACTGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGccc                                                               >  1:85879/1‑136 (MQ=255)
                                                       caCGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGccc                                                               >  1:295414/1‑136 (MQ=255)
                                                                    gctgctGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCTACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGAc                                                  >  2:476157/1‑136 (MQ=255)
                                                                    gctgctGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCTACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGAc                                                  >  2:470431/1‑136 (MQ=255)
                                                                       gctgACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATc                                               <  1:115484/136‑1 (MQ=255)
                                                                            ccTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCTACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCt                                          <  1:151117/136‑1 (MQ=255)
                                                                              ttGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCTGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGc                                        <  2:37431/136‑1 (MQ=255)
                                                                                   cGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAggg                                   >  2:497888/1‑136 (MQ=255)
                                                                                                        cgGCTCATGCGGCCTGAAGGTAGTGGCCTACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAg              <  1:476157/136‑1 (MQ=255)
                                                                                                        cgGCTCATGCGGCCTGAAGGTAGTGGCCTACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAg              <  1:470431/136‑1 (MQ=255)
                                                                                                         ggCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAgg             >  1:465857/1‑136 (MQ=255)
                                                                                                              aTGCGGCCTGAAGGTAGTGGCCGACATTAGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAGGCCGTg        >  2:278771/1‑136 (MQ=255)
                                                                                                              aTGCGGCCTGAAGGTAGTGGCCGACAAGAGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAGGCCGTg        >  2:314006/1‑136 (MQ=255)
                                                                                                                    ccTGAAGGTAGTGGCCGACTGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAGGCCGTGCGCCAg  <  2:85879/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                       
GATGTCACCGGTCCCATGGATGCATTCGCCATCGCCAACCGGTTCTTGCCTGCAGCACGGCAATATCGGCTGCTGACCTTGGCCGAGGGGCAAGCAGTGGTACGCGGCTCATGCGGCCTGAAGGTAGTGGCCGACCGGCGCCTGGAAGACTTGCCGGAGGCCGTCGACCTGTTGCTGGTGCCGGGCGGCCCTGGGGCCTATGACATCGCCCTGCCAGGGATCGAGCGCTGGCTGCCCAAGGCCGTGCGCCAG  >  NC_002947/5990802‑5991053

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: