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

AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAAC  >  NC_000913/1513337‑1513600                                                                                                                                    |                                                                                                                                     aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               <  2:493515/139‑1 (MQ=255) aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               <  2:786128/139‑1 (MQ=255)   cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                             <  1:127629/139‑1 (MQ=255)    tGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGa                                                                                                                            <  2:1419229/139‑1 (MQ=255)       aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               <  1:825470/133‑1 (MQ=255)       aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                               >  2:825470/1‑133 (MQ=255)        aCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtctc                                                                                                                        >  2:407268/1‑139 (MQ=255)           gAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCtat                                                                                                                     >  2:931928/1‑139 (MQ=255)                                      cACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGcc                                                                                          <  1:402291/139‑1 (MQ=255)                                       aCTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCg                                                                                         >  1:975348/1‑139 (MQ=255)                                         ttttATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                       <  1:611582/139‑1 (MQ=255)                                         ttttATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                       <  2:1444830/139‑1 (MQ=255)                                         ttttATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                       <  2:975348/139‑1 (MQ=255)                                             aTCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtct                                                                                                                         >  2:742921/1‑101 (MQ=255)                                             aTCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtct                                                                                                                         <  1:742921/101‑1 (MQ=255)                                                 tCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTgg                                                                               >  1:1038583/1‑139 (MQ=255)                                                  cGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGt                                                                              <  2:80117/139‑1 (MQ=255)                                                            gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGtt                                                                                 <  1:392838/126‑1 (MQ=255)                                                            gATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGtt                                                                                 >  2:392838/1‑126 (MQ=255)                                                                 gggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGttttt                                                               >  2:151584/1‑139 (MQ=255)                                                                    cGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                            <  2:1038583/139‑1 (MQ=255)                                                                      aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                           <  2:1202696/74‑1 (MQ=255)                                                                      aaGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                           >  1:1202696/1‑74 (MQ=255)                                                                       aGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                             <  2:1401760/71‑1 (MQ=255)                                                                       aGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                             >  1:1401760/1‑71 (MQ=255)                                                                            tATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGAt                                                    <  1:407268/139‑1 (MQ=255)                                                                                   gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCagag                                             >  2:157445/1‑139 (MQ=255)                                                                                      cGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGaaa                                          >  2:501916/1‑139 (MQ=255)                                                                                            tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                           <  2:1377622/52‑1 (MQ=255)                                                                                            tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                           >  1:1377622/1‑52 (MQ=255)                                                                                            tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                     >  2:199895/1‑90 (MQ=255)                                                                                            tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                     <  1:199895/90‑1 (MQ=255)                                                                                                                gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGc                >  1:303725/1‑139 (MQ=255)                                                                                                                     aTTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCtac           >  1:433121/1‑137 (MQ=255)                                                                                                                          gCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTCTGATGGa                                                         <  2:1260211/88‑1 (MQ=255)                                                                                                                          gCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTCTGATGGa                                                         >  1:1260211/1‑88 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCGCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                            <  2:901041/81‑1 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                            >  1:901041/1‑81 (MQ=255)                                                                                                                              gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAAc  >  2:1301092/1‑139 (MQ=255)                                                                                                                               tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGAcc       <  2:1162309/133‑1 (MQ=255)                                                                                                                               tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGAcc       >  1:1162309/1‑133 (MQ=255)                                                                                                                                    |                                                                                                                                     AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAAC  >  NC_000913/1513337‑1513600

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