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

CTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATT  >  NC_000913/1621410‑1621662                                                                                                                                     |                                                                                                                         cTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    <  1:492293/139‑1 (MQ=255) cTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    <  1:1479069/139‑1 (MQ=255)  tGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaa                                                                                                                   <  2:1224561/139‑1 (MQ=255)   gCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                  <  1:1366671/139‑1 (MQ=255)    ccATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAAc                                                                                                                 >  1:294109/1‑139 (MQ=255)     cATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACg                                                                                                                <  1:1650329/139‑1 (MQ=255)            gACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    <  2:367778/128‑1 (MQ=255)            gACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    >  1:367778/1‑128 (MQ=255)              ccATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGtt                                                                                                       >  1:1202277/1‑139 (MQ=255)                aTTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTAc                                                                                                     <  2:1439814/139‑1 (MQ=255)                aTTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTAc                                                                                                     >  1:338352/1‑139 (MQ=255)                 tttACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACt                                                                                                    >  2:1669850/1‑139 (MQ=255)                              ccAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGcc                                                                                       <  2:1202277/139‑1 (MQ=255)                                       ccTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGc                                                                              <  1:957162/139‑1 (MQ=255)                                         tGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCgg                                                                            >  2:671533/1‑139 (MQ=255)                                            gtgtCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGttt                                                                         >  1:99013/1‑139 (MQ=255)                                                 gATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTg                                                                    >  2:1043310/1‑139 (MQ=255)                                                  aTCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCt                                                                                                              <  1:931486/96‑1 (MQ=255)                                                  aTCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCt                                                                                                              >  2:931486/1‑96 (MQ=255)                                                  aTCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                   <  1:1069836/139‑1 (MQ=255)                                                     tAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCAt                                                                >  2:813795/1‑139 (MQ=255)                                                         cGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                            <  2:99013/139‑1 (MQ=255)                                                            ttATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTa                                                         <  1:1669850/139‑1 (MQ=255)                                                            ttATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTa                                                         <  2:654960/139‑1 (MQ=255)                                                              aTGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGcc                                                                                       <  2:190923/107‑1 (MQ=255)                                                              aTGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGcc                                                                                       >  1:190923/1‑107 (MQ=255)                                                                 cGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTa                                                              <  2:753527/129‑1 (MQ=255)                                                                 cGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTa                                                              >  1:753527/1‑129 (MQ=255)                                                                 cGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGaa                                                    >  1:556428/1‑139 (MQ=255)                                                                           tGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACg                                          <  1:813795/139‑1 (MQ=255)                                                                           tGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACg                                          >  1:841904/1‑139 (MQ=255)                                                                           tGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACg                                          <  2:338352/139‑1 (MQ=255)                                                                             cgcTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    >  1:1260872/1‑63 (MQ=255)                                                                             cgcTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGa                                                                                                                    <  2:1260872/63‑1 (MQ=255)                                                                                  aCTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTg                                                      >  1:90724/1‑120 (MQ=255)                                                                                  aCTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTg                                                      <  2:90724/120‑1 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGc                                                                              >  2:486924/1‑86 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGc                                                                              <  1:486924/86‑1 (MQ=255)                                                                                            tcgtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCttttt                         >  2:269860/1‑139 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                  <  1:1354795/48‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGaaa                                                                                                                  >  2:1354795/1‑48 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCta                                                                                                             >  2:1228345/1‑53 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCta                                                                                                             <  1:1228345/53‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                       >  2:261039/1‑91 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                       <  1:261039/91‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTa                                                         >  2:1604305/1‑105 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTa                                                         <  1:1604305/105‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                       <  1:1043310/139‑1 (MQ=255)                                                                                              gtTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTg                       >  1:367198/1‑139 (MQ=255)                                                                                                ttttAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGcc                     <  2:294109/139‑1 (MQ=255)                                                                                                        tCGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTaa             <  2:367198/139‑1 (MQ=255)                                                                                                         cGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                   >  2:369028/1‑84 (MQ=255)                                                                                                         cGGTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGc                                                                   <  1:369028/84‑1 (MQ=255)                                                                                                          ggTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                            >  1:21784/1‑90 (MQ=255)                                                                                                          ggTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACt                                                            <  2:21784/90‑1 (MQ=255)                                                                                                           gTTTTGGTATCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTg          <  2:556428/139‑1 (MQ=255)                                                                                                                   aTCCTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTAtt  <  1:269860/139‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACttt                                                          <  2:581398/80‑1 (MQ=255)                                                                                                                      cTGGCGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACttt                                                          >  1:581398/1‑80 (MQ=255)                                                                                                                          cGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                       >  2:424469/1‑63 (MQ=255)                                                                                                                          cGGATAAGTTCAACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTa                                                                       <  1:424469/63‑1 (MQ=255)                                                                                                                                     |                                                                                                                         CTGCCATTTATGACCATTTACTTGAGTCGCCAGTACAGCCTGAGTGTCGATCTAATCGGTTATGCGATGACAATTGCGCTCACTATTGGCGTCGTTTTTAGCCTCGGTTTTGGTATCCTGGCGGATAAGTTCGACAAGAAACGCTATATGTTACTGGCAATTACCGCCTTCGCCAGCGGTTTTATTGCCATTACTTTAGTGAATAACGTGACGCTGGTTGTGCTCTTTTTTGCCCTCATTAACTGCGCCTATT  >  NC_000913/1621410‑1621662

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