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%	149.0 / 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 (25/23);  total (25/23)

GTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCC  >  NC_000913/28139‑28400                                                                                                                                   |                                                                                                                                    gTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTAt                                                                                                                             >  1:1162290/1‑139 (MQ=255) gTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTAt                                                                                                                             >  1:1322248/1‑139 (MQ=255)    gCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCAt                                                                                                                          >  2:726620/1‑139 (MQ=255)      aTTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTa                                                                                                                        >  1:329010/1‑139 (MQ=255)               ggATGTGAAAGGCTTCCAGTAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAAt                                                                                                               >  1:1536447/1‑139 (MQ=255)                gATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATACAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAAtt                                                                                                              >  2:898235/1‑139 (MQ=255)                    tgAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCa                                                                                                          <  2:1074831/139‑1 (MQ=255)                     gAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAg                                                                                                         >  1:480733/1‑139 (MQ=255)                              ccAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtgg                                                                                                >  1:403411/1‑139 (MQ=255)                               cagcagTGGGTGGCTGAGGTGCGGGCTCTAGCGTCGTAACCTGTCACATGTTATTGGCATGAAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggc                                                                                               <  2:1523574/139‑1 (MQ=255)                                    gtgggtggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTggcg                                                                                          <  2:403411/139‑1 (MQ=255)                                    gtgggtggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTggcg                                                                                          >  2:62846/1‑139 (MQ=255)                                       ggtggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGgcggcg                                                                                       >  1:719065/1‑139 (MQ=255)                                       ggtggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGgcggcg                                                                                       <  1:330242/139‑1 (MQ=255)                                            cTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                                  <  1:867655/139‑1 (MQ=255)                                                ggTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgc                                                                              <  2:1322248/139‑1 (MQ=255)                                                ggTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgc                                                                              <  2:1162290/139‑1 (MQ=255)                                                      ggCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                        <  2:1536447/139‑1 (MQ=255)                                                                tcgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                              <  2:480733/139‑1 (MQ=255)                                                                  gtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                              <  1:794300/137‑1 (MQ=255)                                                                  gtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                              >  2:794300/1‑137 (MQ=255)                                                                  gtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAAc                                                            >  1:228305/1‑139 (MQ=255)                                                                     aCCTGTCACATGTTATTTGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGt                                                         <  2:719065/139‑1 (MQ=255)                                                                     aCCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTACTTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGt                                                         <  1:1051977/139‑1 (MQ=255)                                                                     aCCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGt                                                         <  1:62846/139‑1 (MQ=255)                                                                          tCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGt                                                         >  2:651480/1‑134 (MQ=255)                                                                          tCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGt                                                         <  1:651480/134‑1 (MQ=255)                                                                                  tATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAg                                                                                                         >  2:842353/1‑78 (MQ=255)                                                                                  tATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAg                                                                                                         <  1:842353/78‑1 (MQ=255)                                                                                    ttGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                        >  2:382732/1‑109 (MQ=255)                                                                                    ttGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                        <  1:382732/109‑1 (MQ=255)                                                                                          tGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCaa                                                 <  2:1312599/126‑1 (MQ=255)                                                                                          tGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCaa                                                 >  1:1312599/1‑126 (MQ=255)                                                                                          tGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAagag                                    >  1:874381/1‑139 (MQ=255)                                                                                            cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                                  >  2:653387/1‑139 (MQ=255)                                                                                            cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                                  >  2:851046/1‑139 (MQ=255)                                                                                                    catCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATg                          <  2:874381/139‑1 (MQ=255)                                                                                                       cGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAt                                                                                                                      >  2:86989/1‑44 (MQ=255)                                                                                                       cGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAAt                                                                                                                      <  1:86989/44‑1 (MQ=255)                                                                                                         aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtg                                                                                                 >  1:889208/1‑63 (MQ=255)                                                                                                         aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtg                                                                                                 <  2:889208/63‑1 (MQ=255)                                                                                                         aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgcg                                                                             <  2:41676/83‑1 (MQ=255)                                                                                                         aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgcg                                                                             >  1:41676/1‑83 (MQ=255)                                                                                                                       cTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCcgcg       >  1:1092950/1‑139 (MQ=255)                                                                                                                          atatCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaaa                                                             <  1:894756/82‑1 (MQ=255)                                                                                                                          atatCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaaa                                                             >  2:894756/1‑82 (MQ=255)                                                                                                                          atatCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTg    >  1:390953/1‑139 (MQ=255)                                                                                                                            atCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGcc  <  1:851046/139‑1 (MQ=255)                                                                                                                                ctctctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                              <  2:715264/75‑1 (MQ=255)                                                                                                                                ctctctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                              >  1:715264/1‑75 (MQ=255)                                                                                                                                   |                                                                                                                                    GTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCC  >  NC_000913/28139‑28400

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