Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 444,779 C→T 77.2% 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,7790CT77.2% 41.5 / 10.1 28W440* (TGG→TAG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (11/10);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCC  >  CP009273/444635‑444914
                                                                                                                                                |                                                                                                                                       
tCAATCTGCTGCCGCAAACGACGGGTCATGCCCATGAAGCACCGCGCATACCGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGcc                                                                                                                                         <  1:378150/145‑1 (MQ=255)
 cAATGTGCTGGCTCAAACGCCGGGTCATGCCCATGGAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAATCAGAACGCGCGTTGACCCCTGGTTTCGTTCAGTTTGATACCGAACGCTTTAGGcca                                                                                                                                        <  1:272704/145‑1 (MQ=255)
      tgctgcCTCAAACGACGGGTCATGCCCATGAAGTCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGt                                                                                                                                   <  1:182864/145‑1 (MQ=255)
      tgctgGCTCAAACGCGGGGTCATCCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGt                                                                                                                                   <  2:106807/145‑1 (MQ=255)
                       ggtcttcgccACGAAGCCCAGGGCATGCAGTGGCAGAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAag                                                                                                                  <  1:335880/137‑1 (MQ=255)
                          cATGCCCATGAAGCCCAGCGCATACAGTTGCATAAAGGCAACAAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAagcag                                                                                                               <  2:32911/145‑1 (MQ=255)
                                 aTGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGa                                                                                                        <  2:296959/145‑1 (MQ=255)
                                    aaGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCa                                                                                                     >  2:304713/1‑145 (MQ=255)
                                    aaGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGAACa                                                                                                     <  2:190817/145‑1 (MQ=255)
                                       cccAGCGCATACAGTGGCATAAAGCCGACGAAGAAGCCGATGATCCAGAACCAAAACGCGCGTGTACCCCAGGGATCGTTCAGTCTGAAAGCGAGCGCGTTAGGCTGCTAGTAGGTCATCCCAGCGAAGCAGCCGAGAACAAcgc                                                                                                  >  1:325635/1‑145 (MQ=255)
                                               aTACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatg                                                                                          >  1:156835/1‑145 (MQ=255)
                                                tACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGGTCAGTTTGAAACCGCACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatga                                                                                         >  2:49944/1‑145 (MQ=255)
                                                  cAGTGGCAGAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGACCGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATc                                                                                       <  2:218495/145‑1 (MQ=255)
                                                        cATAAGGGCGACGAAGAAGCCGATGATCCAGGACCAGGACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTa                                                                                 <  1:182475/145‑1 (MQ=255)
                                                           aaaGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATgg                                                                              >  1:102905/1‑145 (MQ=255)
                                                                cAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCAGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATggccggg                                                                         >  1:146479/1‑140 (MQ=255)
                                                                  aCGAAGTAGCCTCTGTCCCCGAACCATACCGCNCGTTCACCCGAGGTTTCTTTCAGTTTAAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCNGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTgcg                                                                       <  1:271982/145‑1 (MQ=255)
                                                                   cGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTgcgc                                                                      >  2:391820/1‑145 (MQ=255)
                                                                               tgatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAg                                                          <  2:344259/145‑1 (MQ=255)
                                                                                  tCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCt                                                       <  2:102905/145‑1 (MQ=255)
                                                                                   ccagaaccagaaCGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTg                                                      >  2:421454/1‑145 (MQ=255)
                                                                                    cagaaccagaacGCGCGTTCACCCCAGGTTGCGACCAGTTTGAAACCGAACGCTTTAGGCAACCAGTAGGACATCCCTGCGCAGCGGCCGGAGACCGGGCCGCCAATGTGCAGGTTATGAAAGAGCGCATCCAGGAGCAGACTGt                                                     >  1:80059/1‑145 (MQ=255)
                                                                                      gaaccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTa                                                   >  2:155823/1‑145 (MQ=255)
                                                                                          cagcacGCGGGTTTACCCCAGGTTTCGGTCAGTTTGAAACCGAACGCTGTAGGCCACCAGTAGGGCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGAGGAGCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCa                                               <  1:407402/145‑1 (MQ=255)
                                                                                                gcgcgTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCTATCAGGAACAGGCTGTTATGCAGAACGa                                         <  2:254745/145‑1 (MQ=255)
                                                                                                               gTTGCGGTCAGTTTGAAACCGAACGCTTTAGGCTACCAGCAGGTCATCCCTGCGCGGCGGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTa                          <  2:203428/145‑1 (MQ=255)
                                                                                                                    gTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAACAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGGGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGACCGGTAAGGca                     >  1:82426/1‑144 (MQ=255)
                                                                                                                        aGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTCCGCAATCAGGAACACACTGTTATGCAGAACGAAGTCCGCCCCCGGTACGGCcagca                 >  1:232987/1‑145 (MQ=255)
                                                                                                                                     aCGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTCTGGAAGCGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCAt    >  2:236661/1‑145 (MQ=255)
                                                                                                                                       gCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATcc  >  2:64338/1‑145 (MQ=255)
                                                                                                                                                |                                                                                                                                       
TCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCC  >  CP009273/444635‑444914

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