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.

CATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCC  >  CP009273/444667‑444915
                                                                                                                |                                                                                                                                         
cATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGa                                                                                                             <  1:297867/143‑1 (MQ=255)
      gCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCAcgccgc                                                                                                 <  2:419625/149‑1 (MQ=255)
                             acgccaCGAAGAAGCCGATGACCCAGAACCAGCGCGCGCGTTTACCCCAGGTTTCGTTCAGTGTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTgc                                                                         <  2:77579/144‑1 (MQ=255)
                             aGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGACCCCGTCATgg                                                                                >  1:391312/1‑143 (MQ=255)
                               gCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTAACAGCTTCTCAGCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTgcgc                                                                       >  2:316838/1‑149 (MQ=255)
                                      agaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCGCGCCGCCGATGAGCACGTTATGGAAGTGCGCAATCAgg                                                                >  2:304099/1‑149 (MQ=255)
                                           ccGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAg                                                                 <  1:186498/143‑1 (MQ=255)
                                             gatgatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTACAGTTTGAAACCGAACGCTTTAGGCAACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGc                                                         >  2:78043/1‑149 (MQ=255)
                                                 atCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGtt                                                     >  2:73797/1‑149 (MQ=255)
                                                          cagaacGCGCGTTTACCCCAGGTTTTGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAAc                                            <  2:267435/149‑1 (MQ=255)
                                                                       tACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTcc                                     <  1:170139/143‑1 (MQ=255)
                                                                         ccccAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGCACAGGCTGTTATGCAGAACGACGTCcgc                                   >  1:156034/1‑143 (MQ=255)
                                                                                tttCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGt                            <  1:261212/143‑1 (MQ=255)
                                                                                   cGTTCAGTTTGACACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGCGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACg                         <  1:291352/143‑1 (MQ=255)
                                                                                   cGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGGGCAATCAGGAACAGGCTGTTATGCAGAACGCAGTCCGCGCCCGGTACGGCcagc                   >  2:333845/1‑149 (MQ=255)
                                                                                     ttCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTCTGCAGTGACTTAGTCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag                 >  2:426869/1‑149 (MQ=255)
                                                                                                    aaCGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATccc  <  2:163291/149‑1 (MQ=255)
                                                                                                                |                                                                                                                                         
CATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCC  >  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: