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%	148.7 / NA	48	intergenic (+62/‑105)	rihC/dapB	ribonucleoside hydrolase 3/dihydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/24);  total (24/24)

GCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCA  >  NC_000913/28142‑28385                                                                                                                                |                                                                                                                     gCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCAt                                                                                                           >  1:955174/1‑139 (MQ=255)                 tgAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCa                                                                                           >  2:1112593/1‑139 (MQ=255)                       gCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGttt                                                                                     >  2:693600/1‑139 (MQ=255)                           ccAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtgg                                                                                 >  1:785852/1‑139 (MQ=255)                             agcagTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggct                                                                               <  2:955174/139‑1 (MQ=255)                                    ggtggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCACGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGgcggcg                                                                        <  2:785852/139‑1 (MQ=255)                                           gAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGAt                                                                 <  1:1112593/139‑1 (MQ=255)                                              gTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                             <  1:873145/124‑1 (MQ=255)                                              gTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                             >  2:873145/1‑124 (MQ=255)                                              gTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgcg                                                              >  2:1252961/1‑139 (MQ=255)                                                 cTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTg                                                           >  2:297957/1‑139 (MQ=255)                                                         ggCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATccctctct                                                                                                                    <  1:259910/74‑1 (MQ=255)                                                         ggCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATccctctct                                                                                                                    >  2:259910/1‑74 (MQ=255)                                                         ggCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                             <  1:362547/113‑1 (MQ=255)                                                         ggCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                             >  2:362547/1‑113 (MQ=255)                                                           cgtcgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGa                                                 >  2:552723/1‑139 (MQ=255)                                                             tcgtATCCTGTCACATGTTATTGGCATGCAGTCATTCGTTGACTCACGCCTTTCGCTGATATCCCTCTCTTTTTATCTTTAATCTCTAATTATCGGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                               <  2:655185/139‑1 (MQ=255)                                                             tcgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                               <  1:693600/139‑1 (MQ=255)                                                             tcgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                               <  1:297957/139‑1 (MQ=255)                                                                tAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACg                                            <  1:1252961/139‑1 (MQ=255)                                                                          caTGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGtt                                                                                                                 <  2:1329443/60‑1 (MQ=255)                                                                          caTGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGtt                                                                                                                 >  1:1329443/1‑60 (MQ=255)                                                                               tATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTaa                             >  2:1115257/1‑139 (MQ=255)                                                                                aTTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAAc                            >  2:1197454/1‑139 (MQ=255)                                                                                 ttGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACa                           >  1:983204/1‑139 (MQ=255)                                                                                   ggCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACaaa                         >  1:480833/1‑139 (MQ=255)                                                                                    gCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACaaaa                        <  2:480833/139‑1 (MQ=255)                                                                                    gCATGCAGTCATATATCTACTCCTTCCTTTCACTTATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTGGGCTGGCGTTGTAGCGATGAGCTGGTTACTCGGAGAACGGTATAGGCAAATTAACaaaa                        <  2:803990/139‑1 (MQ=255)                                                                                       tGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                               <  2:463358/113‑1 (MQ=255)                                                                                       tGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                               >  1:463358/1‑113 (MQ=255)                                                                                       tGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAagag                     >  1:524637/1‑139 (MQ=255)                                                                                         cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                             >  2:458630/1‑81 (MQ=255)                                                                                         cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                             <  1:458630/81‑1 (MQ=255)                                                                                         cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                   <  1:1115257/139‑1 (MQ=255)                                                                                         cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                   <  1:1197454/139‑1 (MQ=255)                                                                                                 catCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACaa                          <  1:1205154/124‑1 (MQ=255)                                                                                                 catCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACaa                          >  2:1205154/1‑124 (MQ=255)                                                                                                  atCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGc          <  2:524637/139‑1 (MQ=255)                                                                                                      aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                         <  2:607222/88‑1 (MQ=255)                                                                                                      aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                         >  1:607222/1‑88 (MQ=255)                                                                                                      aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCatga      >  1:353983/1‑139 (MQ=255)                                                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGgcggcg                                                                        <  2:1296504/69‑1 (MQ=255)                                                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGgcggcg                                                                        >  1:1296504/1‑69 (MQ=255)                                                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCaa                                  <  2:1213461/107‑1 (MQ=255)                                                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCaa                                  >  1:1213461/1‑107 (MQ=255)                                                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCa  <  2:353983/139‑1 (MQ=255)                                                                                                                 tCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCa                                                                                           <  2:829747/43‑1 (MQ=255)                                                                                                                 tCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCa                                                                                           >  1:829747/1‑43 (MQ=255)                                                                                                                 tCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                   >  1:464212/1‑115 (MQ=255)                                                                                                                 tCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                   <  2:464212/115‑1 (MQ=255)                                                                                                                                |                                                                                                                     GCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCA  >  NC_000913/28142‑28385

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