Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,357,678	Δ1 bp	intergenic (‑10/+124)	ymjA ← / ← puuP	DUF2543 domain‑containing protein YmjA/putrescine:H(+) symporter PuuP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,357,675	0	C	.	100.0%	115.5 / NA	28	intergenic (‑7/+127)	ymjA/puuP	DUF2543 domain‑containing protein YmjA/putrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (15/13);  total (15/13)

ACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAG  >  NC_000913/1357598‑1357747                                                                              |                                                                         aCAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTGccct                                                                       >  2:261664/1‑80 (MQ=255)         tGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGtt                                                                <  1:261637/79‑1 (MQ=255)         tGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGtt                                                                <  2:261636/79‑1 (MQ=255)         tGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTaa                                                              >  2:261672/1‑81 (MQ=255)            gcATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATgc                                                           >  1:261665/1‑81 (MQ=255)            gcATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATgc                                                           <  2:261638/81‑1 (MQ=255)                  ggTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAg                                                     >  1:261690/1‑81 (MQ=255)                   gTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAg                                                     >  2:261680/1‑80 (MQ=255)                      gcgTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATa                                                 <  2:261639/81‑1 (MQ=255)                                  gCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGgatga                                      >  2:261668/1‑80 (MQ=255)                                  gCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGgatga                                      >  2:261667/1‑80 (MQ=255)                                       atCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCAc                                >  1:261669/1‑81 (MQ=255)                                            aaTATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGt                            >  2:261671/1‑80 (MQ=255)                                            aaTATTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTg                           <  2:261640/81‑1 (MQ=255)                                              tatTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGAt                         <  1:261642/81‑1 (MQ=255)                                              tatTTTAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGAt                         <  2:261641/81‑1 (MQ=255)                                                  ttAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCg                      <  2:261643/80‑1 (MQ=255)                                                  ttAGCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCg                      <  1:261644/80‑1 (MQ=255)                                                     gCGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGt                   <  1:261645/80‑1 (MQ=255)                                                      cGGGATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCa                 >  1:261678/1‑81 (MQ=255)                                                       gggATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCa                 >  1:261685/1‑80 (MQ=255)                                                       gggATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCa                 >  2:261676/1‑80 (MQ=255)                                                        ggATATCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATa               <  2:261646/81‑1 (MQ=255)                                                           tatCGTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACAt            >  2:261681/1‑81 (MQ=255)                                                               gTGGTTCATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCttt        <  2:261647/81‑1 (MQ=255)                                                                     cATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAg  >  1:261683/1‑81 (MQ=255)                                                                     cATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAg  >  1:261682/1‑81 (MQ=255)                                                                     cATTATTTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAg  >  2:261677/1‑81 (MQ=255)                                                                      attattTG‑CCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAg  <  1:261648/80‑1 (MQ=255)                                                                              |                                                                         ACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAG  >  NC_000913/1357598‑1357747

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