Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	450,100	T→C	Y181C (TAC→TGC)	cyoB ←	cytochrome bo3 ubiquinol oxidase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	450,100	0	T	C	100.0%	151.6 / NA	44	Y181C (TAC→TGC)	cyoB	cytochrome bo3 ubiquinol oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/14);  total (30/14)

TGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAA  >  NC_000913/450023‑450175                                                                              |                                                                            tGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACTGACTCCCGGACTGCAct                                                                          >  1:91394/1‑81 (MQ=255) tGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCAct                                                                          >  1:91399/1‑81 (MQ=255) tGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCAct                                                                          >  1:91395/1‑81 (MQ=255)      ccGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTAt                                                                      <  1:91363/80‑1 (MQ=255)        ggTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTcc                                                                   >  2:91401/1‑81 (MQ=255)         gTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTccc                                                                  >  2:91396/1‑81 (MQ=255)         gTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTccc                                                                  >  1:91402/1‑81 (MQ=255)          tAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTccc                                                                  >  1:91393/1‑80 (MQ=255)                  gtACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGt                                                         >  1:91409/1‑81 (MQ=255)                  gtACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGt                                                         >  1:91412/1‑81 (MQ=255)                    aCCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTg                                                        >  2:91400/1‑80 (MQ=255)                     ccTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGa                                                      >  1:91416/1‑81 (MQ=255)                          aCCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGgc                                                 >  2:91407/1‑81 (MQ=255)                            cGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGgcca                                               <  1:91364/81‑1 (MQ=255)                            cGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGgc                                                 >  2:91419/1‑79 (MQ=255)                            cGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGgc                                                 >  1:91406/1‑79 (MQ=255)                              gATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGgccagc                                             >  2:91423/1‑81 (MQ=255)                               aTAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGgccagcc                                            <  2:91365/81‑1 (MQ=255)                                   cTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCagccagcc                                        <  2:91366/81‑1 (MQ=255)                                     ggAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCgg                                      >  2:91404/1‑81 (MQ=255)                                         aCTCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTg                                  >  1:91403/1‑81 (MQ=255)                                           tCCATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGACCCGCCAGCCGGTCTgcg                                >  1:91413/1‑81 (MQ=255)                                            ccATATCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTgcg                                >  2:91411/1‑80 (MQ=255)                                                atCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAAt                           >  2:91426/1‑81 (MQ=255)                                                 tCCAGTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAAt                           <  1:91367/80‑1 (MQ=255)                                                     gTAATCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGcc                      <  2:91368/81‑1 (MQ=255)                                                       aaTCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCa                    >  1:91431/1‑81 (MQ=255)                                                         tCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACg                  >  1:91414/1‑81 (MQ=255)                                                         tCGACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACg                  <  1:91369/81‑1 (MQ=255)                                                           gACACCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGcc                <  2:91370/81‑1 (MQ=255)                                                            acacCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCg               <  1:91371/81‑1 (MQ=255)                                                            acacCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCg               <  1:91372/81‑1 (MQ=255)                                                              acCGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCgag             <  1:91373/81‑1 (MQ=255)                                                                cGACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCga              >  2:91437/1‑78 (MQ=255)                                                                 gACTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCgagag           >  1:91398/1‑80 (MQ=255)                                                                   cTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGaaa        <  2:91374/81‑1 (MQ=255)                                                                   cTCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGaaa        <  1:91375/81‑1 (MQ=255)                                                                    tCCCGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAAc       >  1:91435/1‑81 (MQ=255)                                                                     cccGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACg      >  1:91422/1‑81 (MQ=255)                                                                      ccGGACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGt     >  2:91434/1‑81 (MQ=255)                                                                        ggACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGtt    >  2:91410/1‑80 (MQ=255)                                                                        ggACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGtt    >  2:91432/1‑80 (MQ=255)                                                                        ggACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTa   <  1:91376/81‑1 (MQ=255)                                                                         gACTGCACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTaa  >  2:91397/1‑81 (MQ=255)                                                                              |                                                                            TGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAA  >  NC_000913/450023‑450175

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