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% 59.3 / ‑6.7 23V503L (GTA→TTA) cyoBcytochrome bo3 ubiquinol oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (17/5);  minor base G (0/1);  total (17/6)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAA  >  NC_000913/448996‑449273
                                                                                                                                           |                                                                                                                                          
aaCGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGCCGCCAGTCAGGTCACGATTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGc                                                                                                                                         >  1:629628/1‑143 (MQ=255)
        cGGGGAAGAGGTTGCACACTCCAGCGTACGGCCACCCCACGGGTCGCCAGGCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCgctgg                                                                                                                                 <  2:793250/143‑1 (MQ=255)
                 ggTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAatca                                                                                                                        >  2:60043/1‑143 (MQ=255)
                   ttGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCCGCAatcatc                                                                                                                      >  1:806382/1‑143 (MQ=255)
                                     ggCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGc                                                                                                    >  1:758003/1‑143 (MQ=255)
                                       ccACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCCTGGTGTGGAACTGCgg                                                                                                  >  2:531575/1‑143 (MQ=255)
                                              aCGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACTCAGCGAAATCAGTAATGGACAACTGGCTGAAAACATCAGAATGGTGTGGAACTGCGAGTCAATc                                                                                           >  1:24192/1‑143 (MQ=255)
                                                     gCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATATACATAACGACCCGTAGAATCCTGCGCCCAACCAGTAATGCACCGCAGGCTGCAATCATCAGCATGGTGTGGACCTGCGGCTCAATCTGCTGGc                                                                                    >  1:373744/1‑143 (MQ=255)
                                                                cACGGTGCTCGTCGCCGTCGCGAATAGAAGCGTACATCTGAATAACGAGGCAGAGACTACCCAGCGCAGTCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGGGGAACTGCGGGTCAATCTGCTGGCTCAAACGACgg                                                                         <  2:793650/143‑1 (MQ=255)
                                                                                 tcgcgAATAGAAGCGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGcc                                                        >  1:161259/1‑143 (MQ=255)
                                                                                   gcgAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCa                                                      <  1:691578/143‑1 (MQ=255)
                                                                                   gcgAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCa                                                      >  1:705285/1‑143 (MQ=255)
                                                                                                tACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGAGAAAAATTATGAATTTGTCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTgg                                         >  1:184519/1‑143 (MQ=255)
                                                                                                   aTCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCATTCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCAt                                      >  1:754284/1‑143 (MQ=255)
                                                                                                     cTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATaa                                    >  1:530182/1‑143 (MQ=255)
                                                                                                              gAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGCGGAACTACGGGTCAATCTGCTGGCTCAAACGACGCGTCACGCCCATGACACCCAGCACACACAATCACAAAAAAGCAAccc                           >  2:671703/1‑141 (MQ=255)
                                                                                                                  cAGAGAACACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACgaagaa                       <  1:426370/143‑1 (MQ=255)
                                                                                                                      gaATACCCAGCGCAATCCGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCg                   <  1:60043/143‑1 (MQ=255)
                                                                                                                           cccAGCGCAATCAGTAATGCACCGCTGGCTGCAATAATCAGCATGGTGGGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATGAAGGCAACGAAGAAGGCgatgat              >  1:33975/1‑143 (MQ=255)
                                                                                                                              aGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCAAGCGGAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCAGATGATgca           >  1:113692/1‑143 (MQ=255)
                                                                                                                               gcgcAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCCACGAAGCAGCAGATGATccag          >  2:448258/1‑143 (MQ=255)
                                                                                                                                   agtCAGTAGGGCCCCGCTGGCTGCGATCATCAGCATGGTGTGGGACTGCGGGTCAATCGGCGGGCTCAAACGGCGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATccagaacc      <  2:654718/141‑1 (MQ=255)
                                                                                                                                       aGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACCGCGGGTCAATCTGCTGGCTCAAACGACGGGGCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAGGCCGCTGATccggaacccgcc  >  1:250744/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                          
AACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAA  >  NC_000913/448996‑449273

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: