Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	448,202	A→G	G153G (GGT→GGC)	cyoC ←	cytochrome bo3 ubiquinol oxidase subunit 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	448,202	0	A	G	100.0%	160.4 / NA	47	G153G (GGT→GGC)	cyoC	cytochrome bo3 ubiquinol oxidase subunit 3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/25);  total (22/25)

AACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAA  >  NC_000913/448067‑448338                                                                                                                                        |                                                                                                                                         aacaacAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAg                                                                                                                                      <  2:428447/140‑1 (MQ=255) aacaacAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAg                                                                                                                                      <  1:373547/140‑1 (MQ=255)   caacaGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGaa                                                                                                                                    <  1:264753/140‑1 (MQ=255)     acaGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGt                                                                                                                                  <  2:555216/140‑1 (MQ=255)       aGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGa                                                                                                                                <  1:153350/140‑1 (MQ=255)        gTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGAc                                                                                                                               >  1:420693/1‑140 (MQ=255)           aaCACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTg                                                                                                                            <  1:175188/140‑1 (MQ=255)                    aTCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACcgtg                                                                                                                   >  1:6024/1‑140 (MQ=255)                            cacaTCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGa                                                                                                           >  2:312809/1‑140 (MQ=255)                             acaTCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGAc                                                                                                          <  1:799376/140‑1 (MQ=255)                              caTCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGAcc                                                                                                         >  2:163000/1‑140 (MQ=255)                                         tGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCa                                                                                                        >  2:191308/1‑130 (MQ=255)                                         tGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCa                                                                                                        <  1:191308/130‑1 (MQ=255)                                         tGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGaa                                                                                              >  2:7963/1‑140 (MQ=255)                                         tGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGaa                                                                                              <  1:7963/140‑1 (MQ=255)                                           ccAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACg                                                                                            >  2:311465/1‑140 (MQ=255)                                            cAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGc                                                                                           <  1:776476/140‑1 (MQ=255)                                                    gCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGa                                                                                   <  1:312809/140‑1 (MQ=255)                                                    gCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGa                                                                                   >  1:482391/1‑140 (MQ=255)                                                     cTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGTACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGaa                                                                                  >  1:746610/1‑140 (MQ=255)                                                     cTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGaa                                                                                  >  2:29920/1‑140 (MQ=255)                                                      tCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAg                                                                                 >  2:563664/1‑140 (MQ=255)                                                      tCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAg                                                                                 <  2:517530/140‑1 (MQ=255)                                                            acaTGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTg                                                                           <  1:151795/140‑1 (MQ=255)                                                            acaTGATGCGGGTACGGTTAGTGCTGATCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCTAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTg                                                                           <  2:740934/140‑1 (MQ=255)                                                              atgatgCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCg                                                                         >  1:45446/1‑140 (MQ=255)                                                              atgatgCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCg                                                                         >  1:534471/1‑140 (MQ=255)                                                                gatgCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGAt                                                                       >  1:338559/1‑140 (MQ=255)                                                                        tACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGAccc                                                               <  1:563664/140‑1 (MQ=255)                                                                          cGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATc                                                                      >  2:195535/1‑131 (MQ=255)                                                                          cGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATc                                                                      <  1:195535/131‑1 (MQ=255)                                                                           ggTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATg                                                            <  2:338559/140‑1 (MQ=255)                                                                                gTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCAGCGATCCGGACCCATGCCGtt                                                       >  2:755925/1‑140 (MQ=255)                                                                                  gCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTaa                                                     <  1:417064/140‑1 (MQ=255)                                                                                   cTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAAc                                                    <  2:534471/140‑1 (MQ=255)                                                                                       tCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATc                                                <  1:571339/140‑1 (MQ=255)                                                                                           gCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGt                                            >  2:149705/1‑140 (MQ=255)                                                                                                gACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATgg                                       >  2:777488/1‑140 (MQ=255)                                                                                                gACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATgg                                       >  1:231113/1‑140 (MQ=255)                                                                                                gACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATg                                        >  1:474540/1‑139 (MQ=255)                                                                                                gACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATg                                        <  2:474540/139‑1 (MQ=255)                                                                                                     gCGATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTc                                  >  2:797828/1‑140 (MQ=255)                                                                                                       gATTTGCACCATCAGCACCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCAt                                <  2:420693/140‑1 (MQ=255)                                                                                                                       aCCGCCATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGAt                <  2:308308/140‑1 (MQ=255)                                                                                                                           ccATCCAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAAt            >  2:384947/1‑140 (MQ=255)                                                                                                                                cAGATAAGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGc       <  1:29920/140‑1 (MQ=255)                                                                                                                                     aaGGCCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCaa  <  1:797828/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         AACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAA  >  NC_000913/448067‑448338

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