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,7790CT91.7% 46.3 / ‑3.7 18W440* (TGG→TAG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (7/9);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Not supported by required number of reads on both strands.

GCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCA  >  CP009273/444646‑444899
                                                                                                                                     |                                                                                                                        
gcgcAAACGTCGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATcc                                                                                                           <  1:7002/146‑1 (MQ=255)
gCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATcc                                                                                                           >  1:349704/1‑149 (MQ=255)
            ggTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGCAGAAGCCGATGATCCAGCACCAGAACGCGCGTTTACCCCCGGTTTCGTTCAGTTTGAAACCGAACGCGTTAGGCCACCAGTAGGGCATCCCTGCgg                                                                                                     >  2:176195/1‑142 (MQ=255)
                       tGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGGCCAcg                                                                                    >  1:243958/1‑149 (MQ=255)
                            cccAGCGCATACCGTGGCATAAAGCCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGGTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACgccgcc                                                                               <  1:134054/149‑1 (MQ=255)
                                  gcATCCAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatgat                                                                         <  1:190520/149‑1 (MQ=255)
                                  gcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatgat                                                                         <  1:229168/149‑1 (MQ=255)
                                  gcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatgat                                                                         <  1:355784/149‑1 (MQ=255)
                                       cAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatga                                                                          >  2:64721/1‑143 (MQ=255)
                                               tAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTa                                                                  <  2:162687/143‑1 (MQ=255)
                                                      aaCGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGt                                                           <  2:305340/143‑1 (MQ=255)
                                                               gCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCAGATGATCACGTTATGGAAGTGCGCAATCAGGAACa                                            >  1:205734/1‑149 (MQ=255)
                                                                ccGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCCCGTGATGGAAGCGCGCaagccg                                                 >  2:269778/1‑139 (MQ=255)
                                                                         cagaaccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCt                                        <  2:134125/143‑1 (MQ=255)
                                                                              ccagaacGCGCGTCTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGaa                             <  1:71985/149‑1 (MQ=255)
                                                                                                       tCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcag    >  1:260542/1‑149 (MQ=255)
                                                                                                               tttGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTATTATGCAGAACGAAGTCCGCGCCCGGCACGGCcagcc  >  2:33361/1‑142 (MQ=255)
                                                                                                               cttGAACCCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagca  <  2:315707/142‑1 (MQ=255)
                                                                                                                                     |                                                                                                                        
GCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCA  >  CP009273/444646‑444899

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: