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%	144.0 / NA	45	*378* (TAG→TAA)	ybhS	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/22);  total (23/22)

GCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGC  >  NC_000913/825987‑826249                                                                                                                                     |                                                                                                                                   gcggcgAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATTGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATcc                                                                                                                              >  1:60866/1‑139 (MQ=255)                      aCCTGAATTAGCACGGGTAAAATCAGAACCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                                         >  2:94786/1‑139 (MQ=255)                           aaTTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGt                                                                                                    <  1:889515/139‑1 (MQ=255)                            aTTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTc                                                                                                   >  2:605672/1‑139 (MQ=255)                                     ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGttt                                                                                                              <  2:757648/119‑1 (MQ=255)                                     ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGttt                                                                                                              >  1:757648/1‑119 (MQ=255)                                     ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATa                                                                                          >  1:1638761/1‑139 (MQ=255)                                     ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATa                                                                                          >  1:300974/1‑139 (MQ=255)                                       tAAAATCAGAATTGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATaaa                                                                                        <  2:60866/139‑1 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                          <  2:1638761/139‑1 (MQ=255)                                                     cgcgGGTTCGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGttt                                                                                                              >  1:1539693/1‑103 (MQ=255)                                                     cgcgGGTTCGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGttt                                                                                                              <  2:1539693/103‑1 (MQ=255)                                                                 gTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACa                                                              <  1:94786/139‑1 (MQ=255)                                                                   tCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGc                                                            <  1:605672/139‑1 (MQ=255)                                                                          gcaACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGtt                                                       >  1:1442232/1‑137 (MQ=255)                                                                          gcaACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGtt                                                       <  2:1442232/137‑1 (MQ=255)                                                                              cGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTaccac                                                 <  2:1523068/139‑1 (MQ=255)                                                                                   gCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTaa                                                                                                                                 <  1:1472481/54‑1 (MQ=255)                                                                                   gCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTaa                                                                                                                                 >  2:1472481/1‑54 (MQ=255)                                                                                       ctctTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcaccac                                        <  2:300974/139‑1 (MQ=255)                                                                                                gATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATAtt                               >  1:183835/1‑139 (MQ=255)                                                                                                  ttAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCa                                                                                   >  1:407519/1‑85 (MQ=255)                                                                                                  ttAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCa                                                                                   <  2:407519/85‑1 (MQ=255)                                                                                                    aaCGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCg                           >  2:134068/1‑139 (MQ=255)                                                                                                      cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTc                                                                                                   >  2:339493/1‑65 (MQ=255)                                                                                                      cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTc                                                                                                   <  1:339493/65‑1 (MQ=255)                                                                                                      cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGc                         >  1:871667/1‑139 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGtt                                                                                                               <  2:1479455/48‑1 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGtt                                                                                                               >  1:1479455/1‑48 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAg                                                                                                 <  2:140594/62‑1 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAg                                                                                                 >  1:140594/1‑62 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCg                        <  1:693842/135‑1 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCg                        >  2:693842/1‑135 (MQ=255)                                                                                                                  aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTa                                                     >  2:250510/1‑99 (MQ=255)                                                                                                                  aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTa                                                     <  1:250510/99‑1 (MQ=255)                                                                                                                  aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAT‑ACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCtt            >  1:1682128/1‑139 (MQ=255)                                                                                                                   tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGca                                            <  2:1221243/107‑1 (MQ=255)                                                                                                                   tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGca                                            >  1:1221243/1‑107 (MQ=255)                                                                                                                    gAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGc                         >  2:114932/1‑125 (MQ=255)                                                                                                                    gAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGc                         <  1:114932/125‑1 (MQ=255)                                                                                                                          tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcacca                                         <  1:863132/103‑1 (MQ=255)                                                                                                                          tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcacca                                         >  2:863132/1‑103 (MQ=255)                                                                                                                             tcttctCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGc  <  2:871667/139‑1 (MQ=255)                                                                                                                                 ctcCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGc                                                            <  1:1605112/77‑1 (MQ=255)                                                                                                                                 ctcCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGc                                                            >  2:1605112/1‑77 (MQ=255)                                                                                                                                     |                                                                                                                                   GCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGC  >  NC_000913/825987‑826249

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