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

GCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACG  >  NC_000913/825987‑826254                                                                                                                                     |                                                                                                                                        gcggcgAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATcc                                                                                                                                   >  1:860770/1‑139 (MQ=255) gcggcgAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATcc                                                                                                                                   >  2:1510230/1‑139 (MQ=255) gcggcgAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATcc                                                                                                                                   >  2:194423/1‑139 (MQ=255)                  gATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGtttt                                                                                                                  >  2:52748/1‑139 (MQ=255)                      aCCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                                              >  1:1264499/1‑139 (MQ=255)                            aTTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTc                                                                                                        >  1:32056/1‑139 (MQ=255)                             ttAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCa                                                                                                       <  1:230992/139‑1 (MQ=255)                               aGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAgg                                                                                                     <  2:860770/139‑1 (MQ=255)                                 cACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGcc                                                                                                   >  1:453862/1‑139 (MQ=255)                                  aCGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCa                                                                                                       >  2:1645925/1‑134 (MQ=255)                                  aCGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCa                                                                                                       <  1:1645925/134‑1 (MQ=255)                                  aCGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCg                                                                                                  <  2:1155565/139‑1 (MQ=255)                                  aCGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCg                                                                                                  <  1:1520833/139‑1 (MQ=255)                                  aCGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCg                                                                                                  <  1:1293236/139‑1 (MQ=255)                                      gTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATaa                                                                                              <  1:52748/139‑1 (MQ=255)                                            tCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCa                                                                                        <  2:1332822/139‑1 (MQ=255)                                                aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGttt                                                                                                                        <  2:58114/103‑1 (MQ=255)                                                aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGttt                                                                                                                        >  1:58114/1‑103 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                               <  2:1264499/139‑1 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                               >  2:921767/1‑139 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGGTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                               >  1:941985/1‑139 (MQ=255)                                                                      cgcAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACg                                                              >  1:133788/1‑139 (MQ=255)                                                                        cagcaACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGtt                                                            <  1:98311/139‑1 (MQ=255)                                                                               gACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTaccacc                                                     <  1:1228312/139‑1 (MQ=255)                                                                                 cTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAg                                                   <  2:898164/139‑1 (MQ=255)                                                                                 cTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAg                                                   <  1:1078595/139‑1 (MQ=255)                                                                                   gCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGt                                                                                                                     >  2:841431/1‑71 (MQ=255)                                                                                   gCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGt                                                                                                                     <  1:841431/71‑1 (MQ=255)                                                                                   gCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGca                                                 <  2:83305/139‑1 (MQ=255)                                                                                    ctcctctTGGCGTTTAAACGTCCATCACCGATGACACCTAACCTTCCCATTAATCCAGCCGACGTTTGGTCTACAGCCACGTCAGGCCGATAAACTTCACCGCCGAAACGATCAAAAACAGCACGTTTACCACCAGcac                                                <  2:941985/136‑1 (MQ=255)                                                                                    cAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcac                                                <  1:730879/139‑1 (MQ=255)                                                                                      aCTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcacca                                              <  1:1510230/139‑1 (MQ=255)                                                                                                 aTTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTc                                   >  2:1274129/1‑139 (MQ=255)                                                                                                      cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGc                              <  1:1274129/139‑1 (MQ=255)                                                                                                      cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGc                              <  2:32056/139‑1 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGa                         >  2:909499/1‑139 (MQ=255)                                                                                                             aaGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAAc                       <  2:453862/139‑1 (MQ=255)                                                                                                                  aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGtttt                                                                                                                  >  1:775635/1‑43 (MQ=255)                                                                                                                  aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGtttt                                                                                                                  <  2:775635/43‑1 (MQ=255)                                                                                                                  aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCt                  >  1:282306/1‑139 (MQ=255)                                                                                                                  aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCt                  >  2:154111/1‑139 (MQ=255)                                                                                                                   tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGaa                                        <  2:279339/116‑1 (MQ=255)                                                                                                                   tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGaa                                        >  1:279339/1‑116 (MQ=255)                                                                                                                   tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAAc                       >  2:32184/1‑133 (MQ=255)                                                                                                                   tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAAc                       <  1:32184/133‑1 (MQ=255)                                                                                                                   tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCtt                 >  1:1556580/1‑139 (MQ=255)                                                                                                                       aCATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCa             <  1:921767/139‑1 (MQ=255)                                                                                                                           gCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                                              >  2:847178/1‑38 (MQ=255)                                                                                                                           gCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                                              <  1:847178/38‑1 (MQ=255)                                                                                                                           gCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGt         >  1:1324346/1‑139 (MQ=255)                                                                                                                           gCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGt         >  2:1323266/1‑139 (MQ=255)                                                                                                                               ttctCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                                          <  1:1470715/54‑1 (MQ=255)                                                                                                                               ttctCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                                          >  2:1470715/1‑54 (MQ=255)                                                                                                                                 ctcCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCa                                                                <  2:771236/78‑1 (MQ=255)                                                                                                                                 ctcCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCa                                                                >  1:771236/1‑78 (MQ=255)                                                                                                                                 ctcCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGAc   >  2:1609898/1‑139 (MQ=255)                                                                                                                                 ctcCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGAc   >  1:756324/1‑139 (MQ=255)                                                                                                                                  tcCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACg  <  2:665943/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        GCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACG  >  NC_000913/825987‑826254

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