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

GTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCG  >  NC_000913/28174‑28404                                                                                                |                                                                                                                                        gtgggtggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTggcg                                                                                              >  1:70114/1‑139 (MQ=255)   gggtggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTggcggc                                                                                            >  2:769442/1‑139 (MQ=255)       ggCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAg                                                                                        >  1:851966/1‑139 (MQ=255)           gAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGAt                                                                                    <  1:800752/139‑1 (MQ=255)              gTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgcg                                                                                 >  1:61657/1‑139 (MQ=255)                   ggCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                            <  2:61657/139‑1 (MQ=255)                           cgtcgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGtt                                                                                                                                    >  1:752946/1‑75 (MQ=255)                           cgtcgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCGCTCTGtt                                                                                                                                    <  2:752946/75‑1 (MQ=255)                              cgtAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaaa                                                                 <  2:486877/139‑1 (MQ=255)                                  aCCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGt                                                             <  1:437622/139‑1 (MQ=255)                                  aCCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGt                                                             <  2:70114/139‑1 (MQ=255)                                                 ttGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACa                                              <  1:126119/139‑1 (MQ=255)                                                     cATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                                      <  1:378207/95‑1 (MQ=255)                                                     cATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                                      >  2:378207/1‑95 (MQ=255)                                                       tGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAagag                                        >  1:250209/1‑139 (MQ=255)                                                        gCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                                      <  2:924529/92‑1 (MQ=255)                                                        gCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                                      >  1:924529/1‑92 (MQ=255)                                                        gCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAgaga                                       >  1:657801/1‑139 (MQ=255)                                                         cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                                                <  2:341642/81‑1 (MQ=255)                                                         cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTtggctgg                                                                                                >  1:341642/1‑81 (MQ=255)                                                         cAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                                      <  2:683093/139‑1 (MQ=255)                                                                    cGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCa                            >  2:373990/1‑138 (MQ=255)                                                                    cGACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCa                            <  1:373990/138‑1 (MQ=255)                                                                     gACTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCatg                          <  2:851966/139‑1 (MQ=255)                                                                      aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgc                                                                                  <  2:442728/82‑1 (MQ=255)                                                                      aCTCATGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATgc                                                                                  >  1:442728/1‑82 (MQ=255)                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                            <  1:761472/84‑1 (MQ=255)                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGt                                                                            >  2:761472/1‑84 (MQ=255)                                                                          aTGCCTTTCACTGATATCCCTCTCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCa                     >  2:27663/1‑139 (MQ=255)                                                                                          tccctctcTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCa     >  2:744545/1‑139 (MQ=255)                                                                                             ctctctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATcg  >  2:193265/1‑139 (MQ=255)                                                                                                |                                                                                                                                        GTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCG  >  NC_000913/28174‑28404

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