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%	163.5 / NA	47	S218P (TCC→CCC)	ydcT	putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/27);  total (20/27)

AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGA  >  NC_000913/1513337‑1513593                                                                                                                                    |                                                                                                                              aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        <  1:1064651/139‑1 (MQ=255) aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        <  2:305223/139‑1 (MQ=255)   cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                      <  1:146406/139‑1 (MQ=255)   cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                      <  2:282230/139‑1 (MQ=255)   cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                      <  2:152528/139‑1 (MQ=255)   cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                      <  1:562822/139‑1 (MQ=255)   cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                      <  1:594789/139‑1 (MQ=255)   cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                      <  1:614981/139‑1 (MQ=255)     ggAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                    <  1:761914/139‑1 (MQ=255)      gAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        >  1:139443/1‑134 (MQ=255)      gAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        <  2:139443/134‑1 (MQ=255)       aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        >  2:416769/1‑133 (MQ=255)       aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        <  1:416769/133‑1 (MQ=255)              aaaaCTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGACGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATg                                                                                                           <  2:308817/139‑1 (MQ=255)                             tctcGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        <  2:731163/111‑1 (MQ=255)                             tctaGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        >  1:731163/1‑111 (MQ=255)                                    aTCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTttgttg                                                                                     >  2:896332/1‑139 (MQ=255)                                       aCTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCg                                                                                  >  2:1003179/1‑139 (MQ=255)                                             aTCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCg                                                                            >  2:953203/1‑139 (MQ=255)                                                 tCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTgg                                                                        <  1:547703/139‑1 (MQ=255)                                                   gTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTa                                                                      >  2:178113/1‑139 (MQ=255)                                                     tACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACa                                                                    <  1:953203/139‑1 (MQ=255)                                                           cGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAAt                                                              >  2:130799/1‑139 (MQ=255)                                                            gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        <  2:441793/80‑1 (MQ=255)                                                            gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        >  1:441793/1‑80 (MQ=255)                                                            gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATg                                                             >  1:720499/1‑139 (MQ=255)                                                                    cGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtc                                                                                                                   <  1:1068488/77‑1 (MQ=255)                                                                    cGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtc                                                                                                                   >  2:1068488/1‑77 (MQ=255)                                                                      aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGt                                                                                           <  2:422935/99‑1 (MQ=255)                                                                      aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGt                                                                                           >  1:422935/1‑99 (MQ=255)                                                                      aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATgg                                                   >  1:554468/1‑139 (MQ=255)                                                                                   gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATcccc                                                                                                                            <  1:814609/53‑1 (MQ=255)                                                                                   gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATcccc                                                                                                                            >  2:814609/1‑53 (MQ=255)                                                                                   gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        >  1:63982/1‑57 (MQ=255)                                                                                   gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                        <  2:63982/57‑1 (MQ=255)                                                                                        aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACt                                 <  1:130799/139‑1 (MQ=255)                                                                                          cGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACttt                               <  2:554468/139‑1 (MQ=255)                                                                                                      ttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTgg                                                                        >  1:880263/1‑86 (MQ=255)                                                                                                      ttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTgg                                                                        <  2:880263/86‑1 (MQ=255)                                                                                                                gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGc                          <  1:958601/122‑1 (MQ=255)                                                                                                                gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGc                          >  2:958601/1‑122 (MQ=255)                                                                                                                gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGc         >  2:898980/1‑139 (MQ=255)                                                                                                                  cGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGccc       >  2:991280/1‑139 (MQ=255)                                                                                                                   gCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCt      <  1:896332/139‑1 (MQ=255)                                                                                                                   gCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCt      <  1:1003179/139‑1 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCt      <  1:98606/128‑1 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCt      >  2:98606/1‑128 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGa  >  1:971667/1‑132 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGa  <  2:971667/132‑1 (MQ=255)                                                                                                                                    |                                                                                                                              AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGA  >  NC_000913/1513337‑1513593

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