Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,202 G→A 26.9% I259I (ATC→ATT 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,2020GA26.9% 40.2 / 13.7 26I259I (ATC→ATTetfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base A (5/2);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

GGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAG  >  NC_002947/4748071‑4748333
                                                                                                                                   |                                                                                                                                   
ggAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc                                                                                                                                 >  2:297688/1‑136 (MQ=255)
          tgAAGAGATCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGAt                                                                                                                       >  2:237815/2‑136 (MQ=255)
          cGAAGAGATCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGAt                                                                                                                       <  1:230781/136‑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)
                                                                                                 cTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGc                                <  1:339183/136‑1 (MQ=255)
                                                                                                 cTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGc                                <  2:329573/136‑1 (MQ=255)
                                                                                                     cATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCa                            >  1:301020/1‑136 (MQ=255)
                                                                                                                        ttCGCGCCGTAGATACCGACGGCGATGTACAGCTGTGGCGCTACGTTCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGt         <  2:301020/135‑1 (MQ=255)
                                                                                                                        aTCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGt         <  1:154426/136‑1 (MQ=255)
                                                                                                                        aTCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGt         <  1:180586/136‑1 (MQ=255)
                                                                                                                           gcgcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCgg      <  1:297688/136‑1 (MQ=255)
                                                                                                                           gcgcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCgg      <  1:201151/136‑1 (MQ=255)
                                                                                                                           gcgcCGGAGATACCGACGGCGATGTACAGCTGGGGCGCAACGAACTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCgg      <  2:360814/136‑1 (MQ=255)
                                                                                                                             gcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGcc    <  1:412680/136‑1 (MQ=255)
                                                                                                                             gcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGcc    <  1:407941/136‑1 (MQ=255)
                                                                                                                               cGGAGATACCGACGGCGAAGTACAGCTCTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCaa  >  1:461135/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                   
GGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAG  >  NC_002947/4748071‑4748333

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: