Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,732,541 G→C 23.6% A89G (GCC→GGC)  sucB ← 2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,732,5410GC23.6% 24.7 / 6.4 17A89G (GCC→GGC) sucB2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit
Reads supporting (aligned to +/- strand):  ref base G (7/6);  new base C (2/2);  total (9/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 = 7.22e-01

CCTTGGTGACGCGACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAACCTTGTCAGTCTCG  >  NC_002947/4732392‑4732690
                                                                                                                                                     |                                                                                                                                                     
ccTTGGTGACGCGACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCgg                                                                                                                                                       <  2:68634/150‑1 (MQ=255)
        aCGCGACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCcgcg                                                                                                                                               >  1:349595/1‑150 (MQ=255)
              ccGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGgc                                                                                                                                         >  1:269789/1‑150 (MQ=255)
                    ttACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCCGCTGCGCCTGCCGCGGCAGCCGCacct                                                                                                                                   >  2:64221/1‑147 (MQ=255)
                      aCCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCCGCTGCGCCTGCCGCGGCAGCCGCacctgc                                                                                                                                 <  1:64221/150‑6 (MQ=255)
                        cGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGgcag                                                                                                                               >  2:567995/1‑150 (MQ=255)
                                 cAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGcccc                                                                                                                      >  2:166147/1‑149 (MQ=255)
                                         gCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCgg                                                                                                                                                       <  1:451774/109‑1 (MQ=255)
                                         gCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCgg                                                                                                                                                       >  2:451774/1‑109 (MQ=255)
                                         gCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATc                                                                                                              <  2:349595/150‑1 (MQ=255)
                                                                       ttGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGgtg                                                                                <  1:1709/150‑1 (MQ=255)
                                                                                        tgagATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGAc                                                               <  2:228569/147‑1 (MQ=255)
                                                                                            aTCGGGTCATCTTCGCCGGCGTCAGCGGCCGCAGCGGCCGGCGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTGGGACAGGACGGTGTCGCCCTCGCCCTTGACGATg                                                           >  1:343049/1‑150 (MQ=255)
                                                                                               gggTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCg                                                        >  1:317338/1‑150 (MQ=255)
                                                                                                   cATCTTCGCCGGCGTCGGCGGCAGCAGCGGCCGGAGCAGCGCCAGGTTCGGCTGCCGCGGCACGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCa                                                    <  1:166147/150‑1 (MQ=255)
                                                                                                               cGTCAGCGGCAGCAGCGGCCGCAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCgg                                        <  2:317338/150‑1 (MQ=255)
                                                                                                                                  cGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAg                     <  1:567995/150‑1 (MQ=255)
                                                                                                                                         gcTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAAcc              >  2:396637/1‑150 (MQ=255)
                                                                                                                                               gctGCGCCTGCCGCGGCAGCCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCgg                                        <  1:423365/116‑1 (MQ=255)
                                                                                                                                               gctGCGCCTGCCGCGGCAGCCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCgg                                        >  2:423365/3‑118 (MQ=255)
                                                                                                                                                     gCTGCCGCGGCAGGCGCAGCGGCAGCACCACCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAACCTTGTCAGTCTCg  >  1:562680/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTTGGTGACGCGACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAACCTTGTCAGTCTCG  >  NC_002947/4732392‑4732690

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: