Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,178 G→A 47.6% A267A (GCC→GCT 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,1780GA47.6% ‑4.4 / 22.0 20A267A (GCC→GCTetfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (6/3);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCA  >  NC_002947/4748043‑4748307
                                                                                                                                       |                                                                                                                                 
ggCGGATCAGACCAGGTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCTGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTAGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                                                                                                   <  2:390423/136‑1 (MQ=255)
ggCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                                                                                                   <  2:227756/136‑1 (MQ=255)
ggCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                                                                                                   <  2:186486/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                                                                  >  1:403618/1‑136 (MQ=255)
                                                                 cGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGcc                                                                  >  2:318974/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)
                                                                                                                             cTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGc      <  1:339183/136‑1 (MQ=255)
                                                                                                                             cTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGc      <  2:329573/136‑1 (MQ=255)
                                                                                                                                 cATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCa  >  1:301020/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                 
GGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCA  >  NC_002947/4748043‑4748307

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: