Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	28,269	C→T	intergenic (+62/‑105)	rihC → / → dapB	ribonucleoside hydrolase 3/dihydrodipicolinate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	28,269	0	C	T	100.0%	103.1 / NA	34	intergenic (+62/‑105)	rihC/dapB	ribonucleoside hydrolase 3/dihydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/17);  total (17/17)

GGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGC  >  NC_000913/28138‑28399                                                                                                                                    |                                                                                                                                   ggTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTa                                                                                                                             >  2:775992/1‑139 (MQ=255) ggTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTa                                                                                                                             >  1:252198/1‑139 (MQ=255)     gCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCAt                                                                                                                         >  1:230007/1‑139 (MQ=255)       aTTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTa                                                                                                                       >  1:676378/1‑139 (MQ=255)         tGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAt                                                                                                                     <  2:751408/139‑1 (MQ=255)          ggATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAtt                                                                                                                    >  1:525331/1‑139 (MQ=255)                     tgAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCa                                                                                                         >  2:215206/1‑139 (MQ=255)                                   cagTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                                           <  1:775992/139‑1 (MQ=255)                                   cagTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCAGTTATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCACCGTTTTtggctgg                                                                                           <  1:781965/139‑1 (MQ=255)                                           ggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAg                                                                                   <  1:568734/139‑1 (MQ=255)                                             cTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                                 >  2:252181/1‑139 (MQ=255)                                             cTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                                 <  2:230007/139‑1 (MQ=255)                                               gAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGAt                                                                               <  1:254966/139‑1 (MQ=255)                                                aGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                                           <  1:38749/126‑1 (MQ=255)                                                aGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                                           >  2:38749/1‑126 (MQ=255)                                                    gCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTTGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCCGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCt                                                                          <  2:483319/139‑1 (MQ=255)                                                       ggCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                       <  2:676378/139‑1 (MQ=255)                                                            tGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTActc                                                                  >  2:323906/1‑139 (MQ=255)                                                                ttacgacgcTGTCACATGTTATTGGCATGGAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTATCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaa                                                              >  1:467591/9‑139 (MQ=255)                                                                  cgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaaa                                                            <  1:252181/139‑1 (MQ=255)                                                                           tCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATg                                                   <  2:467591/139‑1 (MQ=255)                                                                                            gCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAgaga                                  <  1:701000/139‑1 (MQ=255)                                                                                             cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAt                                                                                                                     <  1:290873/55‑1 (MQ=255)                                                                                             cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAt                                                                                                                     >  2:290873/1‑55 (MQ=255)                                                                                             cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                                 >  1:513091/1‑139 (MQ=255)                                                                                             cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                                 >  2:697194/1‑139 (MQ=255)                                                                                             cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                                 <  2:252198/139‑1 (MQ=255)                                                                                             aaGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTATAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                                 <  2:525331/138‑1 (MQ=255)                                                                                                          aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCatga                    >  1:140174/1‑139 (MQ=255)                                                                                                          aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGC‑ATTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCatgat                   >  2:190989/1‑139 (MQ=255)                                                                                                              aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                       >  2:693039/1‑84 (MQ=255)                                                                                                              aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                       <  1:693039/84‑1 (MQ=255)                                                                                                              aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCa                >  2:119427/1‑139 (MQ=255)                                                                                                                            tatCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGc  <  2:140174/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   GGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGC  >  NC_000913/28138‑28399

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