Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,813,033 G→A 26.7% A23T (GCC→ACC)  PP_3369 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,813,0330GA26.7% 18.4 / 5.0 15A23T (GCC→ACC) PP_3369LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base A (2/2);  total (8/7)
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

AGGCGGGCATTGAAGCCGGCATGGTCATAAAGGGTTACCCTACCCCTTTATGCACCCGCCACTGCCGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTCGACCGACGTCCCCGCGCTTATGCCCTCACCCACGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGAT  >  NC_002947/3812892‑3813181
                                                                                                                                             |                                                                                                                                                    
aGGCGGGCATTGAAGCCGGCATGGTCATAAAGGGTTACCCTACCCCTTTATGCACCCGCCACTGCCGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGAcc                                                                                                                                              >  2:529916/1‑150 (MQ=255)
    gggCATTGAATCCGGCATGGTCATAAAGGGTTACCCTACCCCTTTATGCACCCGCCACTGCCGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTgg                                                                                                                                          <  2:569604/150‑1 (MQ=255)
               ccGGCATGGTCATAAAGGGTTACCCTACCCCTTTATGCACCCGCCACTGCCGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCAc                                                                                                                               >  1:225736/1‑150 (MQ=255)
                           tAAAGGGTTACCCTACCCCTTTATGCACCCGCCACTGCCGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGcc                                                                                                                   >  2:178736/1‑150 (MQ=255)
                                           cccTTTATGCACCCGCCACTGCCGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCAtcgcttcgc                                                                                                   >  2:128207/1‑150 (MQ=255)
                                                                 cGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTACCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCttga                                                                             >  2:492905/1‑147 (MQ=255)
                                                                 cGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTACCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCttga                                                                             >  2:492923/1‑147 (MQ=255)
                                                                 cGAGTCATTCATGGTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGCCCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCaa                                                                             <  1:178736/150‑1 (MQ=255)
                                                                     tcattcatGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTACCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCttgagctg                                                                         <  1:492923/150‑8 (MQ=255)
                                                                     tcattcatGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTACCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTTGAGCTg                                                                         <  1:492905/150‑1 (MQ=255)
                                                                        ttcatGTTCGGTTGGGAAGACCTGCGCCATTTCTCCGCCGTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTc                                                                      <  2:420541/150‑1 (MQ=255)
                                                                                  aTTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTCGACCGACGTc                                                            <  1:529916/150‑1 (MQ=255)
                                                                                                cgcCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTCGACCGACGTCCCCGCGCTTATGcc                                              <  2:601536/150‑1 (MQ=255)
                                                                                                    aTTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTCGCACGACGTCCCCGCGCTTATGCCCTcc                                          >  2:192125/1‑149 (MQ=255)
                                                                                                              cTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTCGACCGACGTCCCCGCGCTTATGCCCTCACCCACGAAgg                                >  2:203306/1‑150 (MQ=255)
                                                                                                                                            tgcCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTCGACCGACGTCCCCGCGCTTATGCCCTCACCCACGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGAt  <  1:128207/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AGGCGGGCATTGAAGCCGGCATGGTCATAAAGGGTTACCCTACCCCTTTATGCACCCGCCACTGCCGAGTCATTCATGTTCGATTGGGAAGACCTGCGCCATTTCTCCGCCTTCAGCTCGGCCGGCTCGCTGTCGGCTGCTGCCAAGACCCTGGGTGTCGACCACGCCACCGTGGCCCGGCGCATCGCTTCGCTCGAGGCTTCGCTCCAGCTCAAGCTGGTCGACCGACGTCCCCGCGCTTATGCCCTCACCCACGAAGGCCAACGTGTTGCGGAACTCGCCGAGCAGAT  >  NC_002947/3812892‑3813181

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