Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	450,073	C→T	W190* (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
*	NC_000913	450,073	0	C	T	100.0%	132.5 / NA	35	W190* (TGG→TAG)	cyoB	cytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/17);  total (18/17)

AATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAGTTG  >  NC_000913/449925‑450216                                                                                                                                                     |                                                                                                                                                aatCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAg                                                                                                                                               <  2:198298/151‑1 (MQ=255)  atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGt                                                                                                                                              >  1:90671/1‑151 (MQ=255)  atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGt                                                                                                                                              <  1:404434/151‑1 (MQ=255)  atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGt                                                                                                                                              <  1:259107/151‑1 (MQ=255)       tACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATc                                                                                                                                          <  2:378643/150‑1 (MQ=255)                   caGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTcc                                                                                                                              >  1:104036/1‑150 (MQ=255)                     gTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCgg                                                                                                                           >  1:471329/1‑151 (MQ=255)                        aTGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGAc                                                                                                                         <  2:436219/150‑1 (MQ=255)                               ggTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAct                                                                                                                   >  1:228540/1‑149 (MQ=255)                               ggTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAct                                                                                                                   >  1:318832/1‑149 (MQ=255)                                 tAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTAt                                                                                                               <  1:16330/151‑1 (MQ=255)                                 tAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTAt                                                                                                               <  1:421139/151‑1 (MQ=255)                                        tGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGAt                                                                                                        <  1:119732/151‑1 (MQ=255)                                        tGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGAt                                                                                                        <  1:188583/151‑1 (MQ=255)                                        tGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGAt                                                                                                        <  2:404830/151‑1 (MQ=255)                                                  aaCATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGAt                                                                                              <  2:386667/151‑1 (MQ=255)                                                        gTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcca                                                                                        >  2:448952/1‑151 (MQ=255)                                                          cATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgccag                                                                                       >  1:356925/1‑150 (MQ=255)                                                                   tGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTc                                                                             <  2:4903/151‑1 (MQ=255)                                                                   tGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTc                                                                             <  2:246002/151‑1 (MQ=255)                                                                           tCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCaaa                                                                     >  1:382399/1‑151 (MQ=255)                                                                                         aCGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCga                                                       >  1:131492/1‑151 (MQ=255)                                                                                                ggTTCATACAGGTAAGCGACGTACCTATAGCGGATAGCTGGAGATTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGAGGATAGGCCAGCCACCCGGTCTGCGCAAATTCGCCCACGCCGAGAGaaa                                                 >  1:179883/2‑150 (MQ=255)                                                                                                agTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAc                                                                                                                    >  2:198178/1‑83 (MQ=255)                                                                                                agTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAc                                                                                                                    <  1:198175/83‑1 (MQ=255)                                                                                                               gCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATca                                  >  2:200230/1‑150 (MQ=255)                                                                                                                  tcgtACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCacac                               <  1:128581/150‑1 (MQ=255)                                                                                                                     tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaa                            >  2:45202/1‑150 (MQ=255)                                                                                                                     tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGAGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGATAACCAGAATCACACcaa                            >  1:250017/1‑150 (MQ=255)                                                                                                                       ccTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaacaa                         >  2:176314/1‑151 (MQ=255)                                                                                                                          ataCCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACgg                      >  2:467827/1‑151 (MQ=255)                                                                                                                                gATAGCTGGAGACTCCATATCTAGTAATCGACACCGACCCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAAc                 <  2:440120/150‑1 (MQ=255)                                                                                                                                  tAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAg              >  2:62462/1‑151 (MQ=255)                                                                                                                                           aCTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTaa      <  2:19819/150‑1 (MQ=255)                                                                                                                                              ccATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAgttg  >  2:134968/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                AATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAGTTG  >  NC_000913/449925‑450216

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