Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,328,639	T→C	G158G (GGA→GGG)	yciK ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,328,639	0	T	C	100.0%	123.8 / NA	35	G158G (GGA→GGG)	yciK	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/18);  total (17/18)

CACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTG  >  NC_000913/1328506‑1328770                                                                                                                                      |                                                                                                                                    cACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGa                                                                                                                               <  2:46055/140‑1 (MQ=255)     cATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTc                                                                                                                           >  2:635929/1‑140 (MQ=255)        tGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCaa                                                                                                                        >  1:575466/1‑140 (MQ=255)                gcgTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTtgat                                                                                                                <  2:854160/140‑1 (MQ=255)                 cgTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCGCCTGACGTCCAACGCTtgatg                                                                                                               <  2:775020/140‑1 (MQ=255)                      cGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTa                                                                                                          <  2:575466/140‑1 (MQ=255)                           gggTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGAc                                                                                                     >  1:129249/1‑140 (MQ=255)                           gggTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGAc                                                                                                     <  2:115781/140‑1 (MQ=255)                            ggTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTtgatg                                                                                                               >  2:727926/1‑129 (MQ=255)                            ggTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTtgatg                                                                                                               <  1:727926/129‑1 (MQ=255)                            ggTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGAcc                                                                                                    >  1:12980/1‑140 (MQ=255)                                 aTGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGa                                                                                               <  1:550728/140‑1 (MQ=255)                                      gTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCgg                                                                                          >  2:745232/1‑140 (MQ=255)                                      gTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCgg                                                                                          >  2:160593/1‑140 (MQ=255)                                       ttGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGc                                                                                         >  2:363560/1‑140 (MQ=255)                                         gACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGt                                                                                       <  2:393040/140‑1 (MQ=255)                                                aGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGAttt                                                                                <  1:363560/140‑1 (MQ=255)                                                 ggCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTg                                                                               >  2:297800/1‑140 (MQ=255)                                                    gctgctGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGt                                                                            >  1:646567/1‑140 (MQ=255)                                                    gctgctGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGt                                                                            >  1:567878/1‑140 (MQ=255)                                                            tatTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGtaat                                                                         >  1:644308/1‑135 (MQ=255)                                                            tatTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGtaat                                                                         <  2:644308/135‑1 (MQ=255)                                                               tCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAg                                                                 <  1:50489/140‑1 (MQ=255)                                                                           aCCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGt                                                     <  2:567878/140‑1 (MQ=255)                                                                                     tcCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAg                                           >  2:461637/1‑140 (MQ=255)                                                                                            ggTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCAt                                    <  1:635929/140‑1 (MQ=255)                                                                                                 cAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACg                               >  2:506118/1‑140 (MQ=255)                                                                                                   aaTTTCGACGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGtt                             <  1:506118/140‑1 (MQ=255)                                                                                                          aCGCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGc                      >  2:691939/1‑140 (MQ=255)                                                                                                            gCTGCATATGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCAt                    <  2:815044/140‑1 (MQ=255)                                                                                                                 atatGCACCCCAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGt               <  2:713812/140‑1 (MQ=255)                                                                                                                        ccccAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCa        <  2:129249/140‑1 (MQ=255)                                                                                                                          ccAGTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGa      >  2:471140/1‑140 (MQ=255)                                                                                                                             gTTGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCt   >  1:680182/1‑140 (MQ=255)                                                                                                                              ttGGCTCGCCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTg  <  1:160593/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    CACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTG  >  NC_000913/1328506‑1328770

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