Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,041	C→T	A93A (GCC→GCT)	yneK →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,041	0	C	T	100.0%	91.6 / NA	29	A93A (GCC→GCT)	yneK	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/15);  total (14/15)

CAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCC  >  NC_000913/1615918‑1616163                                                                                                                            |                                                                                                                           cAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGacac                                                                                                             >  2:554264/1‑139 (MQ=255)             aaaaaaGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATAtt                                                                                                 <  1:554264/139‑1 (MQ=255)                       tttAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAgg                                                                                       >  2:1217390/1‑139 (MQ=255)                        ttAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgt                                                                                      <  2:442637/139‑1 (MQ=255)                         taataaCCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgta                                                                                     <  1:1217390/139‑1 (MQ=255)                               ccTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgt                                                                                      <  2:134113/132‑1 (MQ=255)                               ccTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgt                                                                                      >  1:134113/1‑132 (MQ=255)                                 ttGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATAtt                                                                                                 <  2:1414814/119‑1 (MQ=255)                                 ttGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATAtt                                                                                                 >  1:1414814/1‑119 (MQ=255)                                                      agagGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgt                                                                                      <  2:736221/109‑1 (MQ=255)                                                      agagGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgt                                                                                      >  1:736221/1‑109 (MQ=255)                                                      agagGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtatttat                                                             <  1:1361479/134‑1 (MQ=255)                                                      agagGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtatttat                                                             >  2:1361479/1‑134 (MQ=255)                                                        agGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATa                                                      >  2:773416/1‑139 (MQ=255)                                                             ttcttTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGacac                                                                                                             <  1:1330908/79‑1 (MQ=255)                                                             ttcttTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGacac                                                                                                             >  2:1330908/1‑79 (MQ=255)                                                                       ggTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTa                                                                         >  1:1312084/1‑105 (MQ=255)                                                                       ggTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTa                                                                         <  2:1312084/105‑1 (MQ=255)                                                                                                aGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGc                                              >  2:126956/1‑107 (MQ=255)                                                                                                aGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGc                                              <  1:126956/107‑1 (MQ=255)                                                                                                aGGTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTcc              >  1:812801/1‑139 (MQ=255)                                                                                                 ggTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTccc             <  1:773416/139‑1 (MQ=255)                                                                                                  gTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGa                                                                                                        <  1:1234912/47‑1 (MQ=255)                                                                                                  gTTCAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGa                                                                                                        >  2:1234912/1‑47 (MQ=255)                                                                                                     cAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTTTTTATTTAGATAATTTTCGCCatgatg                                       <  2:1436219/109‑1 (MQ=255)                                                                                                     cAGTTTAAACTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTTTTTATTTAGATAATTTTCGCCatgatg                                       >  1:1436219/1‑109 (MQ=255)                                                                                                            aaCTCACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGcc  <  2:812801/139‑1 (MQ=255)                                                                                                                cACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCta                                                                  <  2:130402/71‑1 (MQ=255)                                                                                                                cACTTATCCCGCTAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCta                                                                  >  1:130402/1‑71 (MQ=255)                                                                                                                            |                                                                                                                           CAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCC  >  NC_000913/1615918‑1616163

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