Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	451,526	G→A	L29L (CTG→TTG)	cyoA ←	cytochrome bo3 ubiquinol oxidase subunit 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	451,526	0	G	A	100.0%	183.2 / NA	52	L29L (CTG→TTG)	cyoA	cytochrome bo3 ubiquinol oxidase subunit 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/24);  total (28/24)

TTATTGCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGT  >  NC_000913/451390‑451655                                                                                                                                         |                                                                                                                                  ttATTGCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAacgc                                                                                                                                >  2:10059/1‑140 (MQ=255)  tATTGCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCa                                                                                                                               >  2:130750/1‑140 (MQ=255)  tATTGCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCa                                                                                                                               >  1:540545/1‑140 (MQ=255)      gCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAAt                                                                                                                           >  2:487040/1‑140 (MQ=255)         cGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTAc                                                                                                                        <  2:540545/140‑1 (MQ=255)         cGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTAc                                                                                                                        >  2:28384/1‑140 (MQ=255)         cGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTAc                                                                                                                        >  1:777749/1‑140 (MQ=255)         cGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTAc                                                                                                                        >  1:404837/1‑140 (MQ=255)         cGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTAc                                                                                                                        >  1:358255/1‑140 (MQ=255)           cACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAg                                                                                                                      <  2:455919/140‑1 (MQ=255)           cACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAg                                                                                                                      >  1:318820/1‑140 (MQ=255)            aCGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGc                                                                                                                     >  2:6845/1‑140 (MQ=255)              ggTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   <  1:4395/140‑1 (MQ=255)              ggTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   <  1:487040/140‑1 (MQ=255)              ggTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   <  2:358255/140‑1 (MQ=255)              ggTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   <  2:88135/140‑1 (MQ=255)              ggTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   <  2:17854/140‑1 (MQ=255)              ggTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   <  2:776826/140‑1 (MQ=255)                tACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACt                                                                                                                 <  1:504329/140‑1 (MQ=255)                      cAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAg                                                                                                                      >  1:316595/1‑129 (MQ=255)                      cAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAg                                                                                                                      <  2:316595/129‑1 (MQ=255)                          cGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   <  1:80090/128‑1 (MQ=255)                          cGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCa                                                                                                                   >  2:80090/1‑128 (MQ=255)                          cGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACt                                                                                                                 >  2:331810/1‑130 (MQ=255)                          cGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACt                                                                                                                 <  1:331810/130‑1 (MQ=255)                          cGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACa                                                                                                       >  1:149124/1‑140 (MQ=255)                                  aCAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGc                                                                                               <  1:130750/140‑1 (MQ=255)                                        aTCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAaataa                                                                                         >  1:144270/1‑140 (MQ=255)                                                          aCGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACCGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACt                                                                       >  1:303908/1‑140 (MQ=255)                                                             aCAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACCGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACtttt                                                                    <  2:303908/140‑1 (MQ=255)                                                              cAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACAGAGCAATACAGTGCCTGCAAATAATGACAACCATCTCAAACTTTTa                                                                   >  2:280479/1‑140 (MQ=255)                                                               aaTCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACGGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTAt                                                                  <  2:50186/140‑1 (MQ=255)                                                               aaTCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACAGAGCAATACAGTGCCTGCAAATAATGACAACCATCTCAAACTTTTAt                                                                  <  1:280479/140‑1 (MQ=255)                                                                aTCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTAtt                                                                 >  1:674598/1‑140 (MQ=255)                                                                 tCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTg                                                                <  1:143657/140‑1 (MQ=255)                                                                         aGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACTAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACtttt                                                                    >  2:533511/1‑128 (MQ=255)                                                                         aGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACTAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACtttt                                                                    <  1:533511/128‑1 (MQ=255)                                                                           gCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGa                                                      >  2:464612/1‑140 (MQ=255)                                                                           gCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGa                                                      <  2:144270/140‑1 (MQ=255)                                                                             cAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGt                                                    <  2:317926/140‑1 (MQ=255)                                                                                aTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGtctc                                                 >  2:803473/1‑140 (MQ=255)                                                                                            tCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCt                                     <  1:10059/140‑1 (MQ=255)                                                                                                tGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAAt                                 >  2:600709/1‑140 (MQ=255)                                                                                                    cGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCa                             >  1:783787/1‑140 (MQ=255)                                                                                                    cGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCa                             >  1:698938/1‑140 (MQ=255)                                                                                                        gCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACggg                         <  1:803473/140‑1 (MQ=255)                                                                                                        gCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACggg                         >  1:645545/1‑140 (MQ=255)                                                                                                        gCTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACggg                         >  1:105580/1‑140 (MQ=255)                                                                                                         cTCCAGACCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGa                        >  2:767938/1‑140 (MQ=255)                                                                                                               aCCAATCTGTCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCtggt                  >  1:666713/1‑140 (MQ=255)                                                                                                                        tCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTaa         >  2:63240/1‑140 (MQ=255)                                                                                                                        tCCTTTGGGATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTaa         <  1:767938/140‑1 (MQ=255)                                                                                                                               ggATCTAACAACGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAgtgtgt  <  2:683482/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  TTATTGCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGT  >  NC_000913/451390‑451655

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑