Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 449,135 C→A 100% V503L (GTA→TTA)  cyoB ← cytochrome bo3 ubiquinol oxidase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913449,1350CA100.0% 82.6 / NA 30V503L (GTA→TTA) cyoBcytochrome bo3 ubiquinol oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/18);  total (12/18)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAAC  >  NC_000913/448998‑449268
                                                                                                                                         |                                                                                                                                     
cGGAGGCGGGGAAGAGGTGGCCCACTCCAGCGGACGGCCACCCCCCGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCAc                                                                                                                                  <  1:293069/143‑1 (MQ=255)
          gAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGc                                                                                                                        >  1:212628/1‑143 (MQ=255)
                      ccctcaggCGTGCGGCCCCCCCACGGCTCGCCAGTCGGGTCAGGGTTTTGTGCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCAt                                                                                                            <  2:377741/136‑1 (MQ=255)
                         tCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAGACGTACATCTGAATAACGAGGCCGAGAATACCCAGCGCCATCCGTAATGCACCGCTGCCCGCAAACCACAGCacgct                                                                                                         >  2:71057/1‑139 (MQ=255)
                           cAGCGTACGGCCACCCCCCTGGGCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGgtgt                                                                                                       <  2:241921/143‑1 (MQ=255)
                            aGCGTGCCGCCCCCCCACGGGGCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGtgtg                                                                                                      <  2:564395/143‑1 (MQ=255)
                                   ggCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCCTGGTGTGGAACTGc                                                                                               >  1:645795/1‑143 (MQ=255)
                                   ggCCACCCCACGGGTCGCCAGTCAGATCACGGTTATGGTCGCGGTCGCGAATAGAAACGTACATCGGACTAACGAGGCAGAGGATAACCAGCGCAATAAGTAATGCACCGCTGGCTGCAATCAACAGCATGGGATGGAACAGc                                                                                               >  2:28877/1‑143 (MQ=255)
                                     cccccccTCGGGTCGCCAGGCAGGTCAGGGTTCTGGGCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAAGCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCgg                                                                                             <  2:180423/140‑1 (MQ=255)
                                                gTCGCCAGGCCGGCCACGGGTCTGGTCGCGGTCGCGAATAGAACCGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATctgct                                                                                  <  2:481908/143‑1 (MQ=255)
                                                 gcGCCAGTCAGGTCCCGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATctgctg                                                                                 <  2:81465/142‑1 (MQ=255)
                                                   gCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCaa                                                                                        <  2:83978/134‑1 (MQ=255)
                                                   gCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCaa                                                                                        >  1:83978/1‑134 (MQ=255)
                                                     cAGTCAGGTCACGGTTCCGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTc                                                                             <  1:252169/143‑1 (MQ=255)
                                                              cACGGTTCTGGTCGCGGTGCCCACTAGGAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACgg                                                                    <  2:703955/143‑1 (MQ=255)
                                                                             gggcgcgAAGAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCTCCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAg                                                     <  2:633154/140‑1 (MQ=255)
                                                                                 gcgAAGAGAGACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCa                                                 <  2:435737/143‑1 (MQ=255)
                                                                                         aaaCGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATAc                                         >  2:448650/1‑143 (MQ=255)
                                                                                           gcGTACACCTGAGCACCGAGGCAGGGAACACCCGGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCTTACAg                                       <  2:393029/142‑1 (MQ=255)
                                                                                           aCGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGACGCCCAGCGCATACAg                                       >  1:275765/1‑143 (MQ=255)
                                                                                             ccaCATCTCAATAACGGGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTg                                     <  1:41217/141‑1 (MQ=255)
                                                                                                     gAATAACGAGGGAGAGAATACCCCGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAg                             <  1:286062/143‑1 (MQ=255)
                                                                                                      aaTAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAgg                            >  1:283405/1‑143 (MQ=255)
                                                                                                         aaCGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGACGCCCAGCGCATACAGAGGCAGCAAGGGaa                         >  2:483447/1‑143 (MQ=255)
                                                                                                         aaCGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCaa                         <  2:652202/143‑1 (MQ=255)
                                                                                                         aaCGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCaa                         <  2:212628/143‑1 (MQ=255)
                                                                                                               gCAGAGAATACCCAGCGCACTCAGTAATGCACCGCTGGCTGCAACCATCAGCATCGTGTAGAACTGCGGCTCAATCTGCTGGCTCAGACGGCGGGCCATGCCCATGAAGCCCAGCGCACACGGTGGCATAAAGGcgacgacga                   >  2:402868/1‑143 (MQ=255)
                                                                                                               gCAGAGAATACCCAGCGCAATCAGCAATGCACCGCAGGCTGCAAGCCTCAGCATGGGGTGGAACTGCGGGGCAGACCGTCGGCTCGGCCGACGGGGCGAGCCCATGAGGCCCGGCCGACGCAGTGGCAGAGAGGCagagacgg                   >  2:634248/1‑136 (MQ=255)
                                                                                                                            aGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATcca      >  1:199666/1‑143 (MQ=255)
                                                                                                                                cacccAGTAATGCACCGCTGTCTGCAATCATCACCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGCCGGGTCACGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATccagaac  <  1:200380/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
CGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAAC  >  NC_000913/448998‑449268

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: