Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,592,118 G→A 30.8% A440A (GCG→GCA benK → benzoate MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,592,1180GA30.8% 13.5 / 5.8 13A440A (GCG→GCAbenKbenzoate MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (3/6);  new base A (2/2);  total (5/8)
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.40e-01

GCTGATGGGCATCAACCTGCCACTGCAACTCAATTTCATCGCCTTTGCCATCCCTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGAGGTACAGCAGTTCTTCGAGCAGGCCGCAGGCTTTTGTAATGCGGCCGGCAACCCACGCCTCAAAC  >  NC_002947/3591992‑3592250
                                                                                                                              |                                                                                                                                    
gcTGATGGGCATCAACCTGCCACTGCAACTCAATTTCATCGCCTTTGCCATCCCTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCACAGGCCTGa                                                                                                                             >  2:368472/1‑136 (MQ=255)
gcTGATGGGCATCAACCTGCCACTGCAACTCAATTTCATCGCCTTTGCCATCCCTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCACAGGCCTGa                                                                                                                             >  2:475665/1‑136 (MQ=255)
                 tGCCACTGCAACTCAATTTCATCGCCTTTGCCATCCCTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGCCCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGAc                                                                                                            >  1:455223/1‑136 (MQ=255)
                                    cATCGCCTTTGCCATCCCTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCg                                                                                         >  2:91993/1‑136 (MQ=255)
                                      tCGCCTTTGCCATCCCTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTg                                                                                       <  1:476081/136‑1 (MQ=255)
                                         cttttGCCATCCGTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAAc                                                                                    <  1:66494/134‑1 (MQ=255)
                                                         cgcCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCACAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCg                                                                    <  1:368472/136‑1 (MQ=255)
                                                         cgcCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCACAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCg                                                                    <  1:475665/136‑1 (MQ=255)
                                                                          ccATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGAGGTACAGCAGttcttc                                                   <  2:15228/136‑1 (MQ=255)
                                                                          ccATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGAGGTACAGCAGttcttc                                                   <  2:19960/136‑1 (MQ=255)
                                                                                 cGTGCACCTGGCTCTTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGAGGTACAGCAGTTCTTCGAGCAgg                                            <  1:405117/136‑1 (MQ=255)
                                                                                                        cccACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGAGGTACAGCAGTTCTTCGAGCAGGCCGCAGGCTTTTGTAATGcggcc                     <  1:91993/135‑1 (MQ=255)
                                                                                                                    cccTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGAGGTACAGCAGTTCTTCGAGCAGGCCGCAGGCTTTTGTAATGCGGCCGGCAACCCACGc         >  2:371027/1‑136 (MQ=255)
                                                                                                                           ggcgCTGGCCTGAACCGAAGAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGTGGTACAGCAGTTCTTCGAGCAGGCCGCAGGCTTTTGTAATGCGGCCGGCAACCCACGCCTCAAAc  <  1:182262/136‑1 (MQ=255)
                                                                                                                              |                                                                                                                                    
GCTGATGGGCATCAACCTGCCACTGCAACTCAATTTCATCGCCTTTGCCATCCCTGGCGCCATCGCTGCGCTGGCCATGGCCGTGCACCTGGCCAGTGGTCGGCGCCACGCCCAGACCCTGGCGGCGCAGGCCTGAACCGAACAGAACAAGACGAGGTAACCCCCATGACCGTGAACATTTCCCATACTGCCGAGGTACAGCAGTTCTTCGAGCAGGCCGCAGGCTTTTGTAATGCGGCCGGCAACCCACGCCTCAAAC  >  NC_002947/3591992‑3592250

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