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%	215.1 / NA	66	D71N (GAC→AAC)	ydeE	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/36);  total (30/36)

CTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATTCTGTTTTTGCT  >  NC_000913/1621410‑1621673                                                                                                                                     |                                                                                                                                    cTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                               <  2:579649/139‑1 (MQ=255)          aTGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCtata                                                                                                                      <  2:25999/139‑1 (MQ=255)               cATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTa                                                                                                                 <  2:244135/139‑1 (MQ=255)                 tttACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACt                                                                                                               >  2:1231686/1‑139 (MQ=255)                                   aCAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGc                                                                                             <  1:1228930/139‑1 (MQ=255)                                       ccTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGc                                                                                         >  2:1406874/1‑139 (MQ=255)                                         tGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAAc                                                                                                                                   <  2:1208203/95‑1 (MQ=255)                                         tGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCgg                                                                                       >  2:678456/1‑139 (MQ=255)                                            gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                             <  1:728661/98‑1 (MQ=255)                                            gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                             >  2:728661/1‑98 (MQ=255)                                            gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                           >  1:1000232/1‑100 (MQ=255)                                            gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                           <  2:1000232/100‑1 (MQ=255)                                               tCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTAt                                                                                 <  1:1406874/139‑1 (MQ=255)                                                 gATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                               <  2:1215807/91‑1 (MQ=255)                                                 gATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                               >  1:1215807/1‑91 (MQ=255)                                                    cTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCTTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                            <  2:448275/139‑1 (MQ=255)                                                    cTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                            <  1:678456/139‑1 (MQ=255)                                                       aTCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                               >  2:935484/1‑85 (MQ=255)                                                       aTCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                               <  1:935484/85‑1 (MQ=255)                                                         cGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                                       <  1:1302655/139‑1 (MQ=255)                                                              aTGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGt                                                                  <  1:1231686/139‑1 (MQ=255)                                                                      aCAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCAt                                                                           >  2:497789/1‑122 (MQ=255)                                                                      aCAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCAt                                                                           <  1:497789/122‑1 (MQ=255)                                                                          ttGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                             <  2:1477253/68‑1 (MQ=255)                                                                          ttGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                             >  1:1477253/1‑68 (MQ=255)                                                                             cgcTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCg                                                                                        <  1:834954/102‑1 (MQ=255)                                                                             cgcTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCg                                                                                        >  2:834954/1‑102 (MQ=255)                                                                                tCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGtttt                                                                                   <  1:186080/104‑1 (MQ=255)                                                                                tCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGtttt                                                                                   >  2:186080/1‑104 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                                       <  1:907922/104‑1 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                                       >  2:907922/1‑104 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGt                                                                  >  2:739249/1‑109 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGt                                                                  <  1:739249/109‑1 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGGATAACGTGACg                                                     <  1:901564/122‑1 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGGATAACGTGACg                                                     >  2:901564/1‑122 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCttttt                                    >  2:427085/1‑139 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCttttt                                    >  1:845043/1‑139 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACa                                                                                                                                  <  2:275031/43‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACa                                                                                                                                  >  1:275031/1‑43 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                             <  1:110938/48‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                             >  2:110938/1‑48 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                              <  1:256228/95‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                              >  2:256228/1‑95 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                                  <  2:845043/139‑1 (MQ=255)                                                                                                ttttAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGcc                                <  1:849394/139‑1 (MQ=255)                                                                                                          ggTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTAc                                                                                                                >  1:653612/1‑49 (MQ=255)                                                                                                          ggTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTAc                                                                                                                <  2:653612/49‑1 (MQ=255)                                                                                                                     ccTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCAt                                                                           >  1:653187/1‑75 (MQ=255)                                                                                                                     ccTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCAt                                                                           <  2:653187/75‑1 (MQ=255)                                                                                                                     ccTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                                       >  1:331186/1‑79 (MQ=255)                                                                                                                     ccTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                                       <  2:331186/79‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                            <  2:897309/73‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                            >  1:897309/1‑73 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACttt                                                                     >  2:94122/1‑80 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACttt                                                                     <  1:94122/80‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                                  <  1:239985/115‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                                  >  2:239985/1‑115 (MQ=255)                                                                                                                              tAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                            >  2:1343424/1‑65 (MQ=255)                                                                                                                              tAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCa                                                                            <  1:1343424/65‑1 (MQ=255)                                                                                                                              tAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCa                            >  1:708714/1‑113 (MQ=255)                                                                                                                              tAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCa                            <  2:708714/113‑1 (MQ=255)                                                                                                                              tAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATTCTGTTTTTGCt  >  2:725984/1‑139 (MQ=255)                                                                                                                               aaGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTa                                                                    <  1:156687/72‑1 (MQ=255)                                                                                                                               aaGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTa                                                                    >  2:156687/1‑72 (MQ=255)                                                                                                                               aaGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGc                                           <  2:1369511/97‑1 (MQ=255)                                                                                                                               aaGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGc                                           >  1:1369511/1‑97 (MQ=255)                                                                                                                                     |                                                                                                                                    CTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATTCTGTTTTTGCT  >  NC_000913/1621410‑1621673

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