Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,172,551 G→T 12.9% A196E (GCG→GAG)  PP_5410 ← DeoR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,172,5510GT12.9% 77.2 / 3.0 31A196E (GCG→GAG) PP_5410DeoR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (15/12);  new base T (2/2);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCCGCGACCTGGCCCTGGAACGATT  >  NC_002947/6172403‑6172689
                                                                                                                                                    |                                                                                                                                          
gcgGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGc                                                                                                                                           <  2:391121/150‑1 (MQ=255)
      ttGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTATGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGAc                                                                                                                                     >  1:81378/1‑150 (MQ=255)
          aTTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCgg                                                                                                                                 >  1:258734/1‑150 (MQ=255)
                tCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGc                                                                                                                           <  1:144759/150‑1 (MQ=255)
                            gctgGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCg                                                                                                               <  2:258734/150‑1 (MQ=255)
                                 cAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGaa                                                                                                          >  1:447694/1‑150 (MQ=255)
                                    gCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGt                                                                                                       >  2:29381/1‑150 (MQ=255)
                                                      ggTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAg                                                                                     <  2:447694/150‑1 (MQ=255)
                                                       gtgtGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAg                                                                                                     <  1:168971/133‑1 (MQ=255)
                                                       gtgtGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAg                                                                                                     >  2:168971/1‑133 (MQ=255)
                                                               acTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCTCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGcc                                                                            >  1:158429/1‑150 (MQ=255)
                                                               acTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGcc                                                                            >  2:66455/1‑150 (MQ=255)
                                                                    cccAGGGCAGCTCGAGGTTAGTGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCgg                                                                       >  1:423507/1‑150 (MQ=255)
                                                                     ccAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGgc                                                                      >  2:316704/1‑150 (MQ=255)
                                                                                     ttAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTgg                                                      >  1:445757/1‑150 (MQ=255)
                                                                                              ttGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCg                                             >  1:61685/1‑150 (MQ=255)
                                                                                                         tttGTCCGACTCGACCATCGCGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAg                                  <  2:325338/149‑1 (MQ=255)
                                                                                                               cGACTCGACCATCACGATCACTTCACGGGCAACTTCCTCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGt                                                                                   >  1:328350/1‑95 (MQ=255)
                                                                                                               cGACTCGACCATCACGATCACTTCACGGGCAACTTCCTCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGt                                                                                   <  2:328350/95‑1 (MQ=255)
                                                                                                               cGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTg                            >  1:159047/1‑150 (MQ=255)
                                                                                                                    cGACCATCACGATCACTTCACGGGCAACTTCCTCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCcg                       <  2:158429/150‑1 (MQ=255)
                                                                                                                     gACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCcgc                      <  2:423507/150‑1 (MQ=255)
                                                                                                                        cATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAgcgc                                    <  1:367918/133‑1 (MQ=255)
                                                                                                                        cATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAgcgc                                    >  2:367918/1‑133 (MQ=255)
                                                                                                                            aCGATCACTTCACGGGCAACTTCCGCCATGACACGGCTCAGTCTAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCCGCGACCTgg               <  2:81378/150‑1 (MQ=255)
                                                                                                                                   cTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCCGCGACCTGGCCCTGGa        <  1:205465/150‑1 (MQ=255)
                                                                                                                                    ttCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCCGCGACCTGGCCCTGGaa       <  1:29381/150‑1 (MQ=255)
                                                                                                                                     tCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCa                                                                                                                        <  2:33321/36‑1 (MQ=255)
                                                                                                                                     tCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCa                                                                                                                        >  1:33321/1‑36 (MQ=255)
                                                                                                                                     tCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCCGCGACCTGGCCCTGGAAc      >  1:308673/1‑150 (MQ=255)
                                                                                                                                         gtgCAACTTCAGCCATGACCCGGCTCAGCACAAGCAACTCGTTGAAGGGAGTGGTGCCCCGGCGCAGGACGATGCCGTCGGCGCCGATGAATAATTGGTCGGAATCGTAGGAGCGCAGTACCTGCTCCGCGACCTGGCCCTGGAACGAtt  >  2:319019/3‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCGCCGATGAATAATTGGTCGAAATCGTAGGAGCGCAGTACCTGCTCCGCGACCTGGCCCTGGAACGATT  >  NC_002947/6172403‑6172689

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