Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,419,305 C→T 25.0% L78L (CTG→TTG)  PP_3905 → LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,419,3050CT25.0% 31.2 / 4.9 16L78L (CTG→TTG) PP_3905LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base T (2/2);  total (9/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

GGTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATGCTGAAGCTCAAGGTTGCCCGGCGCGCCG  >  NC_002947/4419189‑4419437
                                                                                                                    |                                                                                                                                    
ggTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGt                                                                                                                   >  1:109329/1‑136 (MQ=255)
 gTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTa                                                                                                                  <  2:490973/136‑1 (MQ=255)
     aTGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGccc                                                                                                              >  2:326840/1‑136 (MQ=255)
                             gCAGCCCACACCGAGATCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGcc                                                                                    >  1:277724/1‑136 (MQ=14)
                             gCAGCCCACACCGAGATCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGcc                                                                                    >  1:512614/1‑136 (MQ=14)
                                         cGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCg                                                                          <  1:411078/136‑1 (MQ=255)
                                          gagaTCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGc                                                                        <  2:277724/136‑1 (MQ=18)
                                          gagaTCAAAAATGGACACGATCACTATCGGCGCATGGATAGGCCACCTCGGCCGAGGCCTGGCGCCTCGCGAATTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGc                                                                        <  2:512614/136‑1 (MQ=18)
                                                                   aaTCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGaaa                                                <  1:192484/135‑1 (MQ=255)
                                                                            aTGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGt                                       >  2:27662/1‑136 (MQ=255)
                                                                                aTAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACg                                   >  1:565355/1‑136 (MQ=255)
                                                                                aTAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACg                                   >  2:253519/1‑136 (MQ=255)
                                                                                    gCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAAt                               <  2:347262/136‑1 (MQ=255)
                                                                                     ccATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATg                              <  1:383639/136‑1 (MQ=255)
                                                                                           cGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATGCTGAAg                        >  1:27418/1‑136 (MQ=255)
                                                                                                                 gAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATGCTGAAGCTCAAGGTTGCCCGGCGcgccg  >  1:533607/1‑136 (MQ=255)
                                                                                                                    |                                                                                                                                    
GGTTCATGGCAGCGAGCGTCACGACCAATGCAGCCCACCCCCGAGATCAAATCATGGACACGATCACTATCGGCGCATGGATAGGCCATCTCGGCCGAGGCCTGGCGCCTCGCGAACTGCAATGCCTGTTGGAAGTAGCCCAGGGCTTCACCACCAAGGAAATCGCCAAGCACTTCGGCATCAGCGAGAGCGGCGTAGAGAAACGCATCGGTGACGCAATGCTGAAGCTCAAGGTTGCCCGGCGCGCCG  >  NC_002947/4419189‑4419437

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: