Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,621,542	G→A	D71N (GAC→AAC)	ydeE →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,542	0	G	A	100.0%	109.7 / NA	35	D71N (GAC→AAC)	ydeE	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

CCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATTCTGTTTTTG  >  NC_000913/1621423‑1621671                                                                                                                        |                                                                                                                                  ccATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTcaacaataa                                                                                                                            <  2:946754/127‑4 (MQ=255) ccATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAATaa                                                                                                                            >  1:946754/1‑127 (MQ=255)   aTTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                         >  2:299040/1‑128 (MQ=255)   aTTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                         <  1:299040/128‑1 (MQ=255)   aTTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTAc                                                                                                              >  1:866641/1‑139 (MQ=255)            aGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTa                                                                                                     >  2:540064/1‑139 (MQ=255)              tCGCCAGTACAGCCTGATTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTAcc                                                                                                   >  1:672749/1‑139 (MQ=255)              tCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGG‑TTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                           >  1:74744/1‑114 (MQ=255)              tCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGG‑TTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                           <  2:74744/114‑1 (MQ=255)                               gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAg                                                                                                                              <  2:961125/95‑1 (MQ=255)                               gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAg                                                                                                                              >  1:961125/1‑95 (MQ=255)                                 gtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                                <  2:497535/139‑1 (MQ=255)                                   cGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTAtt                                                                              <  1:1043729/139‑1 (MQ=255)                                     aTCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                            <  1:540064/139‑1 (MQ=255)                                          aTCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTa                                                                       <  2:672749/139‑1 (MQ=255)                                          aTCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTa                                                                       <  2:866641/139‑1 (MQ=255)                                                tATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                           >  1:860368/1‑81 (MQ=255)                                                tATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                           <  2:860368/81‑1 (MQ=255)                                                tATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                         <  2:806414/83‑1 (MQ=255)                                                tATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                         >  1:806414/1‑83 (MQ=255)                                                 aTGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGt                                                                <  2:151425/139‑1 (MQ=255)                                                                               tcgtTTTTAGCCTCGGTTTTGGTATCCTGACGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTATTGAATAACGTGACGCTGGTGGTCCTCttttt                                  >  2:989208/1‑139 (MQ=255)                                                                                 gtTTTTTGCCTCGGTTTTGGTATCCTTACGGATAAATTTAACAAGAAACCCTATATGTTACTGGCCAATTCCGCCTTCGCCAGCGGTTTTATTGCCATTTCTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                                <  1:989208/139‑1 (MQ=255)                                                                                 gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                                <  1:855245/91‑1 (MQ=255)                                                                                 gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                                >  2:855245/1‑91 (MQ=255)                                                                                 gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTa                                                                       <  2:730925/100‑1 (MQ=255)                                                                                 gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTa                                                                       >  1:730925/1‑100 (MQ=255)                                                                                 gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGaa                                                             <  2:267615/110‑1 (MQ=255)                                                                                 gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGaa                                                             >  1:267615/1‑110 (MQ=255)                                                                                             ggTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATg                                                                                                                  >  2:258208/1‑45 (MQ=255)                                                                                             ggTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATg                                                                                                                  <  1:258208/45‑1 (MQ=255)                                                                                             ggTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACt                    >  1:79956/1‑139 (MQ=255)                                                                                                                 tAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGc                                                                                       >  2:350281/1‑52 (MQ=255)                                                                                                                 tAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGc                                                                                       <  1:350281/52‑1 (MQ=255)                                                                                                                 tAAGTTCAACAAGAAACGCTATATGTTACTGGAAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATTCTGTTTTTg  <  1:808433/137‑1 (MQ=255)                                                                                                                 tAAGTTCAACAAGAAACGCTATATGTTACTGGAAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATTCTGTTTTTg  >  2:808433/1‑137 (MQ=255)                                                                                                                        |                                                                                                                                  CCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATTCTGTTTTTG  >  NC_000913/1621423‑1621671

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