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

GGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCC  >  CP009273/444645‑444914
                                                                                                                                      |                                                                                                                                       
ggCTCAAACCATGGGTCATGCCCATGAAGCCCGGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCCGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAg                                                                                                                                 <  1:49137/143‑1 (MQ=255)
 gCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATcc                                                                                                                          >  2:596253/1‑149 (MQ=255)
     aaaCGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGc                                                                                                                      >  2:37364/1‑149 (MQ=255)
         tcctggTCAAGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCt                                                                                                                        <  1:206237/139‑1 (MQ=255)
             tgTCATGCCCATGAAGCCCAGCGCATACAGTGGCATACAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCATGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagcagc                                                                                                              <  2:197766/148‑1 (MQ=255)
              gTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGGCATCCCTGCGAAGCAGcc                                                                                                             >  2:349148/1‑149 (MQ=255)
                   gCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGa                                                                                                        <  2:164036/149‑1 (MQ=255)
                        tcaaGCCCAGCCCATACAGTGGCATAACGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAg                                                                                                         <  1:165537/141‑1 (MQ=255)
                                      tACAGTGGCATAAAGTCAACGAAAAAGCCGAGGATACAGAACCAGAACGCGCGTTTACCCCAGGTTTGGTTCAGTTTGAAACCGAACGCTTTAGGCTCCCAGTAGGTCATCCCTGCGAAACAACCGACGACCACGccggcgct                                                                                           >  1:52626/1‑141 (MQ=255)
                                      tACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCGCGCCGACGATGAGCAc                                                                                     >  2:81503/1‑149 (MQ=255)
                                           tGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTAt                                                                                >  2:551828/1‑149 (MQ=255)
                                           tGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCAAAGACCACGCCGCCGATGATCACGTGAt                                                                                >  2:93662/1‑149 (MQ=255)
                                            gcaagaaaGGCAACGAGGAAGCCTATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCAc                                                                                     <  1:361088/138‑1 (MQ=255)
                                             gCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTACGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATgg                                                                              >  2:627433/1‑149 (MQ=255)
                                                tAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCCCGTTATGGAAg                                                                           >  2:53230/1‑149 (MQ=255)
                                                          gaagaaGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGGTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCCGCCGACGGCCGCGCCGCCGATGATAAAGTTATGGAAGTGCGCAAGCa                                                                 >  2:216792/1‑149 (MQ=255)
                                                             gaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGCAAAGGAAGAGCGCaa                                                                    >  1:320466/1‑143 (MQ=255)
                                                                gCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCAAGCCGCAGATGATCACGTTATGGAAGTGCGCAAGAAGGaagc                                                           >  2:586062/1‑147 (MQ=255)
                                                                   gacgatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGGGTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGc                                                        <  2:276893/146‑1 (MQ=255)
                                                                              accagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTa                                                   <  1:32814/143‑1 (MQ=255)
                                                                               ccagtgcgcgcgTTTTCCCCGGTTTTCGTTGAGTTTGACAGCGAGCGGTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGaa                                            <  2:66631/143‑1 (MQ=255)
                                                                               ccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTCTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTGAGGGAACAGCGCAATCAGGAACCGGCTggact                                                  >  1:276808/1‑139 (MQ=255)
                                                                                         cgTTTACCCCAGGTTTCGTTCAGGTTGAAACCCAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGaa                                        <  1:242724/143‑1 (MQ=255)
                                                                                                    ggTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACgg                       <  2:250436/149‑1 (MQ=255)
                                                                                                          gTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagca                 <  2:42005/149‑1 (MQ=255)
                                                                                                             cAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAgcagca              <  2:440595/149‑1 (MQ=255)
                                                                                                                   gAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGCTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCgccagcc              >  1:452234/1‑142 (MQ=255)
                                                                                                                               tttAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATcc  <  1:538814/143‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCC  >  CP009273/444645‑444914

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