Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,513,468	T→C	S218P (TCC→CCC)	ydcT →	putative ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,513,468	0	T	C	100.0%	251.0 / NA	69	S218P (TCC→CCC)	ydcT	putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/33);  total (36/33)

TGCAGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACC  >  NC_000913/1513334‑1513595                                                                                                                                       |                                                                                                                                tGCAGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCcccc                                                                                                                             >  1:1237430/1‑139 (MQ=255)    aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  2:711719/139‑1 (MQ=255)      cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                        <  1:1589446/139‑1 (MQ=255)        ggAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGaa                                                                                                                      <  2:142028/139‑2 (MQ=255)          aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  1:1228331/1‑133 (MQ=255)          aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  2:1228331/133‑1 (MQ=255)              gAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCtat                                                                                                                >  1:823489/1‑139 (MQ=255)               aaaaaaCTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  2:1503930/128‑1 (MQ=255)               aaaaaaCTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  1:1503930/1‑128 (MQ=255)                    aCTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATgcgc                                                                                                          >  1:39722/1‑139 (MQ=255)                          aCAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCCCTATATGCGCCCgcgc                                                                                                    <  2:275391/139‑1 (MQ=255)                                tctcGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCg                                                                                              >  2:828034/1‑139 (MQ=255)                                  tcGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCa                                                                                                   >  2:493298/1‑132 (MQ=255)                                  tcGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCa                                                                                                   <  1:493298/132‑1 (MQ=255)                                        tCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGc                                                                                      >  1:1704752/1‑139 (MQ=255)                                                            cACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGa                                                                  <  2:1704752/139‑1 (MQ=255)                                                               gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATcccc                                                                                                                              <  2:999216/76‑1 (MQ=255)                                                               gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATcccc                                                                                                                              >  1:999216/1‑76 (MQ=255)                                                               gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATg                                                               >  1:988793/1‑139 (MQ=255)                                                                    gggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  1:164955/1‑75 (MQ=255)                                                                    gggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  2:164955/75‑1 (MQ=255)                                                                    gggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGttttt                                                          >  2:377573/1‑139 (MQ=255)                                                                     ggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                        >  2:964623/1‑76 (MQ=255)                                                                     ggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                        <  1:964623/76‑1 (MQ=255)                                                                     ggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                  >  1:1467323/1‑114 (MQ=255)                                                                     ggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                  <  2:1467323/114‑1 (MQ=255)                                                                     ggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCg                                                                                    >  2:41899/1‑112 (MQ=255)                                                                     ggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCg                                                                                    <  1:41899/112‑1 (MQ=255)                                                                       cGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                       <  2:823489/139‑1 (MQ=255)                                                                         aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGATGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  2:1228668/1‑70 (MQ=255)                                                                         aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  1:1228668/70‑1 (MQ=255)                                                                         aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGCCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                <  2:1577996/112‑1 (MQ=255)                                                                         aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGCCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                >  1:1577996/1‑112 (MQ=255)                                                                           gCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  1:125315/68‑1 (MQ=255)                                                                           gCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  2:125315/1‑68 (MQ=255)                                                                            cGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACt                                                  >  1:37738/1‑139 (MQ=255)                                                                               tATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGAt                                               >  1:1085596/1‑139 (MQ=255)                                                                                   gATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCa                                           <  2:37738/139‑1 (MQ=255)                                                                                      gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCagag                                        >  2:457133/1‑139 (MQ=255)                                                                                       tCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAgaga                                       <  1:828034/139‑1 (MQ=255)                                                                                           aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACt                                   >  2:326891/1‑139 (MQ=255)                                                                                           aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACt                                   >  1:1274918/1‑139 (MQ=255)                                                                                              gtgtGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTtg                                <  2:1274918/139‑1 (MQ=255)                                                                                               tgtgGCGGCTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTtgt                               >  1:1272000/1‑139 (MQ=255)                                                                                                 tgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  2:1123747/1‑46 (MQ=255)                                                                                                 tgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  1:1123747/46‑1 (MQ=255)                                                                                                  ggCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGc                            <  1:457133/139‑1 (MQ=255)                                                                                                   gCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGcc                                                                                     <  1:1578685/81‑1 (MQ=255)                                                                                                   gCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGcc                                                                                     >  2:1578685/1‑81 (MQ=255)                                                                                                   gCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCa                           <  2:1237430/139‑1 (MQ=255)                                                                                                      gTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  1:1058901/1‑41 (MQ=255)                                                                                                      gTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  2:1058901/41‑1 (MQ=255)                                                                                                       ttttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGt                                                              <  1:956312/100‑1 (MQ=255)                                                                                                       ttttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGt                                                              >  2:956312/1‑100 (MQ=255)                                                                                                        tttCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          >  2:45275/1‑39 (MQ=255)                                                                                                        tttCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                          <  1:45275/39‑1 (MQ=255)                                                                                                         ttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACgg                     >  1:1153874/1‑139 (MQ=255)                                                                                                          tCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACggg                    <  2:1272000/139‑1 (MQ=255)                                                                                                                atGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCtt              >  2:531080/1‑139 (MQ=255)                                                                                                                   gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGc           >  1:821278/1‑139 (MQ=255)                                                                                                                    aCGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGt                                                              <  1:464965/87‑1 (MQ=255)                                                                                                                    aCGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGt                                                              >  2:464965/1‑87 (MQ=255)                                                                                                                      gCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCt        <  2:39722/139‑1 (MQ=255)                                                                                                                           gAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGAc   <  1:326891/139‑1 (MQ=255)                                                                                                                            aGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGAcc  >  1:998565/1‑139 (MQ=255)                                                                                                                               aGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCg                                                                   <  2:1195526/71‑1 (MQ=255)                                                                                                                               aGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCg                                                                   >  1:1195526/1‑71 (MQ=255)                                                                                                                                  tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCg                                                                                    <  2:1156304/51‑1 (MQ=255)                                                                                                                                  tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCg                                                                                    >  1:1156304/1‑51 (MQ=255)                                                                                                                                  tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGa    <  1:1098276/131‑1 (MQ=255)                                                                                                                                  tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGa    >  2:1098276/1‑131 (MQ=255)                                                                                                                                       |                                                                                                                                TGCAGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACC  >  NC_000913/1513334‑1513595

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: ‑