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%	161.3 / NA	43	W190* (TGG→TAG)	cyoB	cytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/22);  total (21/22)

ATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAGT  >  NC_000913/449926‑450214                                                                                                                                                    |                                                                                                                                              atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGt                                                                                                                                            <  1:158198/151‑1 (MQ=255)   cAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTaa                                                                                                                                          >  1:9139/1‑151 (MQ=255)          ttCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGacac                                                                                                                                   <  1:76110/151‑1 (MQ=255)          ttCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGaca                                                                                                                                    >  2:231787/1‑150 (MQ=255)            cgcgCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCg                                                                                                                                 >  2:143095/1‑151 (MQ=255)                           ccAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAc                                                                                                                  <  2:191322/151‑1 (MQ=255)                              ggTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTAct                                                                                                                 >  2:305670/1‑149 (MQ=255)                               gTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTActc                                                                                                                <  2:194064/149‑1 (MQ=255)                                tAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTAt                                                                                                             <  2:136992/151‑1 (MQ=255)                                tAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTAt                                                                                                             <  1:167418/151‑1 (MQ=255)                                     aCTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCg                                                                                                        <  2:205551/151‑1 (MQ=255)                                     aCTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCg                                                                                                        <  2:168407/151‑1 (MQ=255)                                       tGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGAt                                                                                                      <  2:311201/151‑1 (MQ=255)                                              ttGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTg                                                                                               <  2:94766/151‑1 (MQ=255)                                                 aaCATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAGTCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGAt                                                                                            <  2:296648/151‑1 (MQ=255)                                                 aaCATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGAt                                                                                            <  2:46820/151‑1 (MQ=255)                                                 aaCATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGAt                                                                                            <  1:277389/151‑1 (MQ=255)                                                       gTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcca                                                                                      >  1:128178/1‑151 (MQ=255)                                                       gTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcc                                                                                       <  2:186813/150‑1 (MQ=255)                                                       gTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcc                                                                                       <  2:72263/150‑1 (MQ=255)                                                       gTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcc                                                                                       <  2:98893/150‑1 (MQ=255)                                                                 gTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGt                                                                            >  2:62655/1‑151 (MQ=255)                                                                          tCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCaaa                                                                   >  1:102964/1‑151 (MQ=255)                                                                          tCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCaaa                                                                   >  2:76261/1‑151 (MQ=255)                                                                           cTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAAt                                                                  >  2:255624/1‑151 (MQ=255)                                                                           cTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAAt                                                                  >  1:38411/1‑151 (MQ=255)                                                                            ttCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAAtt                                                                 <  2:44358/151‑1 (MQ=255)                                                                                                gTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACg                                             >  2:215190/1‑151 (MQ=255)                                                                                                 ttGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAAc                                              <  1:218459/149‑1 (MQ=255)                                                                                                     tACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTaa                                         >  2:159657/1‑150 (MQ=255)                                                                                                       ccGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCa                                      >  1:71208/1‑151 (MQ=255)                                                                                                        cGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAg                                     >  1:267863/1‑151 (MQ=255)                                                                                                              gCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATca                                >  2:161763/1‑150 (MQ=255)                                                                                                                 tcgtACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCacac                             <  1:271581/150‑1 (MQ=255)                                                                                                                 tcgtACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCacac                             <  1:157755/150‑1 (MQ=255)                                                                                                                    tACCTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaa                          >  1:196176/1‑150 (MQ=255)                                                                                                                      ccTATACCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACcaacaa                       >  1:141706/1‑151 (MQ=255)                                                                                                                         ataCCGGATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACgg                    >  1:167399/1‑151 (MQ=255)                                                                                                                               gATAGCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAAc               <  2:280958/150‑1 (MQ=255)                                                                                                                                   gCTGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGaa          >  1:133686/1‑151 (MQ=255)                                                                                                                                     tGGAGACTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGc        <  1:131589/151‑1 (MQ=255)                                                                                                                                           cTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAgt  >  1:226103/1‑151 (MQ=255)                                                                                                                                           cTCCATATCTAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAgt  >  2:133673/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                              ATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGCAAATTCGCCCACGCCGAGAGAAACGTTAACCAGAATCACACCAACAACGGTAAACCAGAAGCTTAAGT  >  NC_000913/449926‑450214

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