Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	772,049	T→C	F198L (TTC→CTC)	cydA →	cytochrome d terminal oxidase, subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	772,049	0	T	C	100.0%	145.4 / NA	39	F198L (TTC→CTC)	cydA	cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/18);  total (21/18)

GGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGA  >  NC_000913/771902‑772188                                                                                                                                                    |                                                                                                                                            ggCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTc                                                                                                                                           >  2:347863/1‑150 (MQ=255)    tGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATcc                                                                                                                                       >  2:308874/1‑151 (MQ=255)       aTGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCg                                                                                                                                    >  1:134557/1‑151 (MQ=255)          cAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTa                                                                                                                                 <  1:511314/151‑1 (MQ=255)            aaaCCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATc                                                                                                                               <  1:136569/151‑1 (MQ=255)             aaCCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCa                                                                                                                              <  2:531684/151‑1 (MQ=255)                  aaTCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAgcgc                                                                                                                          >  1:44704/1‑150 (MQ=255)                          ccGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTatat                                                                                                                 >  2:299033/1‑151 (MQ=255)                           cGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATg                                                                                                                <  1:345496/151‑1 (MQ=255)                                  aaCTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGaa                                                                                                           >  1:97252/1‑149 (MQ=255)                                        gAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGt                                                                                                    <  1:158583/150‑1 (MQ=255)                                        gAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTg                                                                                                   <  2:69001/151‑1 (MQ=255)                                          aaCTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGAc                                                                                                 <  1:28892/151‑1 (MQ=255)                                              aTGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGActtcg                                                                                             <  2:176257/151‑1 (MQ=255)                                                  gTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGActtcgcctt                                                                                         >  1:213543/1‑151 (MQ=255)                                                         gaTGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTaa                                                                                   >  2:95239/1‑150 (MQ=255)                                                                 gCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCtt                                                                          <  2:431355/151‑1 (MQ=255)                                                                   ttCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTg                                                                        >  1:492166/1‑151 (MQ=255)                                                                     ctcCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGCCGTGACTTCGCCTTCGCTAAACGCTCCTTTGc                                                                       >  2:538170/1‑150 (MQ=255)                                                                       ccGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCt                                                                      <  2:356021/149‑1 (MQ=255)                                                                             tGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGc                                                              >  2:382115/1‑151 (MQ=255)                                                                                       cccGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCgg                                                     >  2:306322/1‑150 (MQ=255)                                                                                       cccGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGt                                                    <  2:573506/151‑1 (MQ=255)                                                                                         cGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTa                                                   >  1:108100/1‑150 (MQ=255)                                                                                          ggTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATg                                                 >  2:364439/1‑151 (MQ=255)                                                                                                 cAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGgctgct                                           <  2:492156/150‑1 (MQ=255)                                                                                                          ttCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTtctgtctg                                 >  1:109439/1‑151 (MQ=255)                                                                                                           gcGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTctgtctgt                                >  2:455395/2‑151 (MQ=255)                                                                                                                      gTAGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCt                       <  2:409397/149‑1 (MQ=255)                                                                                                                        aGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATg                                                 <  1:594846/121‑1 (MQ=255)                                                                                                                        aGCGTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATg                                                 >  2:594831/1‑121 (MQ=255)                                                                                                                           gTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGtt                         <  2:256356/142‑1 (MQ=255)                                                                                                                           gTCTGGTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGtt                         >  1:256363/1‑142 (MQ=255)                                                                                                                               ggTTATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGtgatga              <  2:566292/149‑1 (MQ=255)                                                                                                                                  tATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATcc          <  1:336149/150‑1 (MQ=255)                                                                                                                                  tATGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCg         >  1:70299/1‑151 (MQ=255)                                                                                                                                   aTGTGACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCgg        <  1:45647/151‑1 (MQ=255)                                                                                                                                       gACTGGCGCGATGCTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGGATCCGGCTAc    >  2:307667/1‑151 (MQ=255)                                                                                                                                            gcgcGATGCTCATCCTCGGAATCAGCGCATGGTATATGCTGCAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTGTCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGa  >  2:10081/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                            GGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGTCGTGACTTCGCCTTCGCTAAACGCTCCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGA  >  NC_000913/771902‑772188

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