Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,169 T→C 37.5% K270K (AAA→AAG 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,1690TC37.5% 11.7 / 22.8 24K270K (AAA→AAGetfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (7/8);  new base C (6/3);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

CAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGC  >  NC_002947/4748041‑4748303
                                                                                                                                |                                                                                                                                      
cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCCCGATCACTTTGGAGTCTTTCATGc                                                                                                                                 >  2:407941/1‑136 (MQ=255)
cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGc                                                                                                                                 >  2:412680/1‑136 (MQ=255)
  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  <  2:329573/136‑1 (MQ=255)
                                                                                                                               cTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGc  <  1:339183/136‑1 (MQ=255)
                                                                                                                                |                                                                                                                                      
CAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGC  >  NC_002947/4748041‑4748303

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: