Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 444,779 C→T 100% W440* (TGG→TAG)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273444,7790CT86.8% 28.9 / ‑0.6 17W440* (TGG→TAG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (5/9);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Not supported by required number of reads on both strands.

CATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCC  >  CP009273/444667‑444915
                                                                                                                |                                                                                                                                         
cATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGa                                                                                                             <  1:297867/143‑1 (MQ=255)
      gCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCAcgccgc                                                                                                 <  2:419625/149‑1 (MQ=255)
                             acgccaCGAAGAAGCCGATGACCCAGAACCAGCGCGCGCGTTTACCCCAGGTTTCGTTCAGTGTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTg‑c                                                                         <  2:77579/144‑1 (MQ=255)
                             aGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGACCCCGTCATgg                                                                                >  1:391312/1‑143 (MQ=255)
                               gCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTAACAGCTTCTCAGCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTg‑cgc                                                                       >  2:316838/1‑149 (MQ=255)
                                      agaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCGCGCCGCCGATGAGCACGTTATGGAAGTG‑CGCAATCAgg                                                                >  2:304099/1‑149 (MQ=255)
                                           ccGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAg                                                                 <  1:186498/143‑1 (MQ=255)
                                             gatgatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTACAGTTTGAAACCGAACGCTTTAGGCAACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGc                                                         >  2:78043/1‑149 (MQ=255)
                                                 atCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGCTGtt                                                     >  2:73797/1‑149 (MQ=255)
                                                          cagaacGCGCGTTTACCCCAGGTTTTGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGCTGTTATGCAGAAc                                            <  2:267435/149‑1 (MQ=255)
                                                                       tACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTcc                                     <  1:170139/143‑1 (MQ=255)
                                                                         ccccAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGCACAGGCTGTTATGCAGAACGACGTCcgc                                   >  1:156034/1‑143 (MQ=255)
                                                                                tttCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGt                            <  1:261212/143‑1 (MQ=255)
                                                                                   cGTTCAGTTTGACACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGCG‑CGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACg                         <  1:291352/143‑1 (MQ=255)
                                                                                   cGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑GGCAATCAGGAACAGGCTGTTATGCAGAACGCAGTCCGCGCCCGGTACGGCcagc                   >  2:333845/1‑149 (MQ=255)
                                                                                     ttCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTCT‑GCAGTGACTTAGTCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag                 >  2:426869/1‑149 (MQ=255)
                                                                                                    aaCGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATccc  <  2:163291/149‑1 (MQ=255)
                                                                                                                |                                                                                                                                         
CATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTG‑CGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCC  >  CP009273/444667‑444915

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑