Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	451,411	C→T	W67* (TGG→TAG)	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,411	0	C	T	100.0%	101.1 / NA	26	W67* (TGG→TAG)	cyoA	cytochrome bo3 ubiquinol oxidase subunit 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/13);  total (13/13)

AGACCACAGCTTCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTA  >  NC_000913/451334‑451480                                                                              |                                                                      aGACCACAGCTTCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAgg                                                                    <  2:87085/81‑1 (MQ=255)    ccACAGCTTCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGa                                                                 <  1:87087/81‑1 (MQ=255)    ccACAGCTTCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGa                                                                 <  1:87086/81‑1 (MQ=255)     cacaGCTTCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGa                                                                 >  1:87122/1‑80 (MQ=255)       caGCTTCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAAc                                                              >  1:87139/1‑81 (MQ=255)            tCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCaa                                                           >  2:87123/1‑79 (MQ=255)                     ttGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAg                                                <  2:87088/81‑1 (MQ=255)                      tGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATc                                                   >  1:87131/1‑77 (MQ=255)                        gAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGAtt                                             <  2:87090/81‑1 (MQ=255)                        gAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGAtt                                             <  2:87089/81‑1 (MQ=255)                        gAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGAtt                                             >  1:87144/1‑81 (MQ=255)                           tgtgACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCg                                          <  1:87091/81‑1 (MQ=255)                            gtgACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCgg                                         >  1:87140/1‑81 (MQ=255)                               aCCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGaa                                      >  1:87133/1‑81 (MQ=255)                                ccAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAAt                                     <  1:87092/81‑1 (MQ=255)                                  aGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATaa                                   >  1:87126/1‑81 (MQ=255)                                   gTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAAc                                  >  1:87138/1‑81 (MQ=255)                                   gTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAAc                                  <  1:87093/81‑1 (MQ=255)                                      cGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGa                                >  1:87137/1‑80 (MQ=255)                                            gTACTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCaa                         >  1:87147/1‑81 (MQ=255)                                              aCTTAGCATCTTTATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACa                       >  2:87142/1‑81 (MQ=255)                                                          tATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCa               <  1:87095/77‑1 (MQ=255)                                                          tATTGCTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAAt            <  2:87094/80‑1 (MQ=255)                                                               cTCGCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGt       <  1:87096/80‑1 (MQ=255)                                                                  gCACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGt   <  1:87097/81‑1 (MQ=255)                                                                   cACGGTACTTCTAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTa  >  2:87132/1‑81 (MQ=255)                                                                              |                                                                      AGACCACAGCTTCCACTTTATTGGAGTGTGACCAGTTCGGGCTGTACTTAGCATCTTTATTGCTCGCACGGTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTA  >  NC_000913/451334‑451480

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