Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	915,224	C→T	intergenic (‑367/+128)	lysO ← / ← aqpZ	L‑lysine/thialysine exporter/water channel AqpZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	915,224	0	C	T	100.0%	55.9 / NA	19	intergenic (‑367/+128)	lysO/aqpZ	L‑lysine/thialysine exporter/water channel AqpZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/12);  total (7/12)

ATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGC  >  NC_000913/914994‑915468                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     atTCTTGCGCTTTAACTACAAACTGCATCGATGTTGGCGGCGTCAAAATCGCCGAGGCGTTACCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGAGTGGCCGCCCTCTGCGCGCTCCTTGCGcc                                                                                                                                                                                                                                    >  2:242995/1‑249 (MQ=14)            ttAATTACAAACTGCCCCGATGTGGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCGGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATg                                                                                                                                                                                                                        <  1:55174/250‑1 (MQ=255)                              gATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                                                                                                                                     >  1:147829/1‑251 (MQ=255)                                                aTCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGatgatg                                                                                                                                                                                    >  2:125904/1‑250 (MQ=255)                                                                  gAAGGCCGGGGCAGCCCCCATGGATGGGGGCTGAGGGCGCGTTTTACAGGAATGTTCCCCCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                                                                                                                                 <  1:147146/251‑1 (MQ=255)                                                                                                                aaGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGGGGGCCCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa                                                                                                                   <  1:228382/250‑1 (MQ=255)                                                                                                                       tttcccccccgcccaccccGGTAGCCGTAAGGGATAAGGCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc                                                                                                             <  2:41649/234‑1 (MQ=255)                                                                                                                              cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc                                                                                                       >  2:46550/1‑249 (MQ=255)                                                                                                                                  cccGCCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAgg                                                                                                 <  1:222857/251‑1 (MQ=255)                                                                                                                                                  tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGCCGGGGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTACCTCATGAAAGGAGCGCAACAAGAGGCTGAATCAACATATAACAACATCTTAAAAAAAGGACAGACATGAAGCCACGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGac                                                                                 >  2:206396/1‑251 (MQ=25)                                                                                                                                                     gggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGaccac                                                                              <  2:147829/251‑1 (MQ=255)                                                                                                                                                          aaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATa                                                                         <  2:2710/251‑1 (MQ=255)                                                                                                                                                             gcGAGGGCACCGCGAAGCGGCGAGTTTGTGCGCAAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAA‑GGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgc                                                                     <  1:130041/250‑1 (MQ=38)                                                                                                                                                                    cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGa                                                               <  2:44834/251‑1 (MQ=255)                                                                                                                                                                           cAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTgg                                                        >  2:223074/1‑251 (MQ=255)                                                                                                                                                                                  gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCaccacc                                                 <  2:97326/251‑1 (MQ=255)                                                                                                                                                                                                   cccGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                                                                                                                              >  2:79050/1‑109 (MQ=255)                                                                                                                                                                                                                  ggcggcggcGATTGGCCGCCCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgcc                 <  2:4602/251‑1 (MQ=255)                                                                                                                                                                                                                                 ccgccCTCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGc  <  1:125904/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     ATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGC  >  NC_000913/914994‑915468

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