Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,162 T→G 50.0% K273Q (AAA→CAA) ‡ etfA ← electron transfer flavoprotein subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,748,1620TG50.0% 0.4 / 31.1 22K273Q (AAA→CAA) ‡etfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base G (8/3);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCC  >  NC_002947/4748041‑4748279
                                                                                                                         |                                                                                                                     
cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCCCGATCACTTTGGAGTCTTTCATGc                                                                                                         >  2:407941/1‑136 (MQ=255)
cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGc                                                                                                         >  2:412680/1‑136 (MQ=255)
  ggCGGATCAGACCAGGTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCTGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTAGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                                                                       <  2:390423/136‑1 (MQ=255)
  ggCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                                                                       <  2:186486/136‑1 (MQ=255)
  ggCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                                                                       <  2:227756/136‑1 (MQ=255)
                    ttCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATcgcg                                                                                     <  2:291647/136‑1 (MQ=255)
                              ggAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc                                                                           >  2:297688/1‑136 (MQ=255)
                                        tgAAGAGATCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGAt                                                                 >  2:237815/2‑136 (MQ=255)
                                        cGAAGAGATCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGAt                                                                 <  1:230781/136‑1 (MQ=255)
                                                  accgacgaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATcg                                                                                       >  1:294148/3‑104 (MQ=255)
                                                  accgacgaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATcg                                                                                       <  2:294148/102‑1 (MQ=255)
                                                     gACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCtccggt                                                    >  1:417825/1‑131 (MQ=255)
                                                     gACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCtccggt                                                    >  1:414346/1‑131 (MQ=255)
                                                            ccGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGa                                             >  2:180586/1‑136 (MQ=255)
                                                                   cGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGcc                                      >  2:318974/1‑136 (MQ=255)
                                                                   cGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGcc                                      >  1:403618/1‑136 (MQ=255)
                                                                                     gtgcTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGcca                    >  2:493401/3‑134 (MQ=255)
                                                                                      tgcTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGccag                   >  2:369310/2‑133 (MQ=255)
                                                                                        cttcttCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGCCGACCTGccagtc                 <  2:417825/136‑6 (MQ=255)
                                                                                         ttcttcGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCg                <  1:402828/136‑1 (MQ=255)
                                                                                                  cTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACg       <  1:441569/136‑1 (MQ=255)
                                                                                                       tgATCGCCACGATCACTATGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGcc  >  1:403409/1‑136 (MQ=255)
                                                                                                                         |                                                                                                                     
CAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCC  >  NC_002947/4748041‑4748279

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: