Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	826,119	C→T	*378* (TAG→TAA)	ybhS ←	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	826,119	0	C	T	100.0%	141.0 / NA	44	*378* (TAG→TAA)	ybhS	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/21);  total (23/21)

CCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGA  >  NC_000913/826009‑826231                                                                                                               |                                                                                                                 ccTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                       <  2:1416737/138‑1 (MQ=255) ccTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                       >  1:1416737/1‑138 (MQ=255)    gAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGAGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACg                                                                                   >  1:1382895/1‑139 (MQ=255)     aaTTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGAc                                                                                                     >  2:1336727/1‑120 (MQ=255)     aaTTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGAc                                                                                                     <  1:1336727/120‑1 (MQ=255)          gCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAgg                                                                              >  1:1039884/1‑138 (MQ=255)          gCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAgg                                                                              <  2:1039884/138‑1 (MQ=255)          gCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGc                                                                             >  1:232402/1‑139 (MQ=255)               ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATa                                                                        >  2:1252889/1‑139 (MQ=255)                gTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATaa                                                                       >  2:323856/1‑139 (MQ=255)                gTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATaa                                                                       >  2:489546/1‑139 (MQ=255)                    aaTCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                   <  1:323856/139‑1 (MQ=255)                    aaTCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                   <  1:489546/139‑1 (MQ=255)                    aaTCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                   <  1:1252889/139‑1 (MQ=255)                          aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAg                                                                                        <  2:795352/112‑1 (MQ=255)                          aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAg                                                                                        >  1:795352/1‑112 (MQ=255)                          aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAgg                                                                              <  2:1108022/122‑1 (MQ=255)                          aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAgg                                                                              >  1:1108022/1‑122 (MQ=255)                          aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCAccgc                                                             >  2:1374149/1‑139 (MQ=255)                               cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                        <  2:47396/139‑1 (MQ=255)                               cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                        <  2:232402/139‑1 (MQ=255)                               cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                        >  2:1016643/1‑139 (MQ=255)                               cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                        >  1:1312181/1‑139 (MQ=255)                                              cgcgCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTg                                                                                                >  2:1457879/1‑84 (MQ=255)                                              cgcgCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTg                                                                                                <  1:1457879/84‑1 (MQ=255)                                              cgcgCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCa                                         <  2:1312181/139‑1 (MQ=255)                                                        cGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTaccac                               <  1:1374149/139‑1 (MQ=255)                                                               aaCTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGc                                                                             <  1:47817/86‑1 (MQ=255)                                                               aaCTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGc                                                                             >  2:47817/1‑86 (MQ=255)                                                                           aTTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTc            >  1:483943/1‑139 (MQ=255)                                                                             tAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCa                                         >  2:824151/1‑108 (MQ=255)                                                                             tAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCa                                         <  1:824151/108‑1 (MQ=255)                                                                                cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGAc                                                                                                     <  2:1467564/45‑1 (MQ=255)                                                                                cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGAc                                                                                                     >  1:1467564/1‑45 (MQ=255)                                                                                         gCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                   <  1:1492971/70‑1 (MQ=255)                                                                                         gCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                   >  2:1492971/1‑70 (MQ=255)                                                                                         gCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGa  >  2:749150/1‑135 (MQ=255)                                                                                         gCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGa  <  1:749150/135‑1 (MQ=255)                                                                                              gAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCaaa                                                <  2:54258/84‑1 (MQ=255)                                                                                              gAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCaaa                                                >  1:54258/1‑84 (MQ=255)                                                                                                    tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                        >  1:526798/1‑70 (MQ=255)                                                                                                    tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                        <  2:526798/70‑1 (MQ=255)                                                                                                    tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcaccac                      <  2:350561/104‑1 (MQ=255)                                                                                                    tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcaccac                      >  1:350561/1‑104 (MQ=255)                                                                                                               |                                                                                                                 CCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGA  >  NC_000913/826009‑826231

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