Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	914,159	T→C	E233E (GAA→GAG)	lysO ←	L‑lysine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	914,159	0	T	C	100.0%	94.2 / NA	26	E233E (GAA→GAG)	lysO	L‑lysine exporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/13);  total (13/13)

GATGTGGCACCGCATAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGTTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCG  >  NC_000913/914082‑914231                                                                              |                                                                         gATGTGGCACCGCATAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCa                                                                       >  1:174036/1‑81 (MQ=255)  aTGTGGCACCGCATAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCAc                                                                      <  1:174009/81‑1 (MQ=255)              aTAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCAt                                                          >  1:174042/1‑81 (MQ=255)              aTAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCAt                                                          >  1:174045/1‑81 (MQ=255)               tAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCAt                                                          >  1:174041/1‑80 (MQ=255)               tAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCAt                                                          >  1:174049/1‑80 (MQ=255)                   cccAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTaaaa                                                     <  2:174010/81‑1 (MQ=255)                     cAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTaaaaaa                                                   >  1:174039/1‑81 (MQ=255)                        tGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAAcgc                                                <  1:174011/81‑1 (MQ=255)                          cAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAAcgccg                                              <  2:174012/81‑1 (MQ=255)                             tAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCgcgc                                           <  2:174013/81‑1 (MQ=255)                              agagCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCt                                          >  1:174051/1‑81 (MQ=255)                               gagCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTc                                         <  2:174014/81‑1 (MQ=255)                               gagCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTc                                         >  2:174025/1‑81 (MQ=255)                                     cTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGAtt                                   <  1:174015/81‑1 (MQ=255)                                       gcggcgAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTAc                                 >  2:174048/1‑81 (MQ=255)                                         ggcgAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCg                               >  2:174040/1‑81 (MQ=255)                                          gcgAATCAGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCg                               >  2:174050/1‑80 (MQ=255)                                                 aGCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCaaaa                       <  2:174016/81‑1 (MQ=255)                                                  gCCCAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAg                      <  2:174017/81‑1 (MQ=255)                                                   cccAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCaa                         >  1:174046/1‑77 (MQ=255)                                                    ccAGGGATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGAt                    <  2:174018/81‑1 (MQ=255)                                                       gggATCAACATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCg                 <  2:174019/81‑1 (MQ=255)                                                               cATAATAGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATaa         <  2:174020/81‑1 (MQ=255)                                                                ataataGCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATaa         >  1:174043/1‑80 (MQ=255)                                                                      gCAATCAGCTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCg  <  2:174021/81‑1 (MQ=255)                                                                              |                                                                         GATGTGGCACCGCATAAGCCCAGTGCAGTAGAGCGGCTGCGGCGAATCAGCCCAGGGATCAACATAATAGCAATCAGTTCACGGGCCAGATCATTAAAAAACGCCGCGCTCCCGATTACCGGACCAAAAGATTCGGTCAATAAAATACCG  >  NC_000913/914082‑914231

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