Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,093,493	A→G	H186R (CAC→CGC)	yqgE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,093,493	0	A	G	100.0%	119.5 / NA	36	H186R (CAC→CGC)	yqgE	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/20);  total (16/20)

TGGAACAAGAAATTCTCGATAACGCGTGGCTAACGGCCCCGGCAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACACGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGT  >  NC_000913/3093349‑3093640                                                                                                                                                 |                                                                                                                                                    tGGAACAAGAAATTCTCGATAACGCGTGGCTAACGGCCCCGGCAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCtg                                                                                                                                               >  2:422650/1‑151 (MQ=255)       aaGAAATTCTCGATAACGCGTGGCTAACGGCCCCGGCAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGt                                                                                                                                         <  1:428228/151‑1 (MQ=255)         gAAATTCTCGATAACGCGTGGCTAACGGCCCCGGCAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTgg                                                                                                                                       >  2:336745/1‑151 (MQ=255)                        gcgTGGCTAACGGCCCCGGCAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCtt                                                                                                                        >  2:434611/1‑151 (MQ=255)                                aaCGGCCCCGGCAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTc                                                                                                                 <  1:660877/150‑1 (MQ=255)                                           cAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGc                                                                                                     <  1:264632/151‑1 (MQ=255)                                                 tGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGc                                                                                               >  2:523768/1‑151 (MQ=255)                                                  gAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGc                                                                                               <  1:562158/150‑1 (MQ=255)                                                             tCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGc                                                                                   >  1:287144/1‑151 (MQ=255)                                                             tCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGc                                                                                   >  1:399427/1‑151 (MQ=255)                                                                       gATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCAtt                                                                          >  1:5156/1‑150 (MQ=255)                                                                           gCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCgg                                                                     <  1:452931/151‑1 (MQ=255)                                                                               aCCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCAcc                                                                 >  1:394152/1‑151 (MQ=255)                                                                                  gCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCAcc                                                                 >  2:661144/1‑148 (MQ=255)                                                                                   gtGGCGCGAGGCGGCAGGACTGGGTGGAGTGGATATTCTCACCATGCGTGGTGTGGCAGGCCGCGCCTGATGAAAAAAAACTTACTCGACTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTc                                                             >  1:591479/2‑151 (MQ=255)                                                                                     ggCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTc                                                             <  1:504123/149‑1 (MQ=255)                                                                                     ggCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGc                                                           <  1:169870/151‑1 (MQ=255)                                                                                       cgcgAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGccc                                                         <  1:392701/151‑1 (MQ=255)                                                                                           aGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTg                                                     >  2:107350/1‑151 (MQ=255)                                                                                                      cTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTa                                           >  2:334245/1‑150 (MQ=255)                                                                                                              tgtgGATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAgg                                  >  2:435086/1‑151 (MQ=255)                                                                                                                 ggATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGAc                                <  1:152646/150‑1 (MQ=255)                                                                                                                 ggATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGAc                                <  1:322622/150‑1 (MQ=255)                                                                                                                  gATATTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACgg                              <  2:394284/151‑1 (MQ=255)                                                                                                                   atatTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGt                             <  2:376932/151‑1 (MQ=255)                                                                                                                     atTCTCACCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTAc                           <  2:653558/151‑1 (MQ=255)                                                                                                                           aCCATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGa                     <  2:287143/151‑1 (MQ=255)                                                                                                                            ccATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGAc                    <  1:334242/151‑1 (MQ=255)                                                                                                                             cATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACt                   <  1:19427/151‑1 (MQ=255)                                                                                                                             cATGCCTGGTGTGGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACt                   >  1:224131/1‑151 (MQ=255)                                                                                                                                       gtgGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTa          <  2:114128/150‑1 (MQ=255)                                                                                                                                       gtgGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTAt         <  1:633260/151‑1 (MQ=255)                                                                                                                                       gtgGCAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCA‑CAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATc        <  1:359515/151‑1 (MQ=255)                                                                                                                                           cAGGACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCg       >  1:265203/1‑149 (MQ=255)                                                                                                                                             ggACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCg   >  1:73365/1‑151 (MQ=255)                                                                                                                                              gACGCGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGt  <  1:669713/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TGGAACAAGAAATTCTCGATAACGCGTGGCTAACGGCCCCGGCAGATCTGAATATTCTGTTCAAAACGCCGATTGCCGACCGCTGGCGCGAGGCGGCAAAACTGATTGGTGTGGATATTCTCACCATGCCTGGTGTGGCAGGACACGCCTGATGAGTGGAACCTTACTCGCCTTCGACTTCGGCACCAAAAGCATTGGCGTAGCGGTCGGCCAACGCATTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGT  >  NC_000913/3093349‑3093640

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