Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,252,298	C→T	intergenic (+15/‑208)	yjbI → / → ubiC	pentapeptide repeats‑containing protein YjbI/chorismate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,252,298	0	C	T	100.0%	111.6 / NA	31	intergenic (+15/‑208)	yjbI/ubiC	pentapeptide repeats‑containing protein YjbI/chorismate lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/17);  total (14/17)

CGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTC  >  NC_000913/4252222‑4252375                                                                             |                                                                              cgcTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCg                                                                           <  2:801977/81‑1 (MQ=255)    ttttCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTc                                                                        >  1:802005/1‑81 (MQ=255)            aTTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATc                                                                >  2:802014/1‑81 (MQ=255)              ttttACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCt                                                              >  2:802015/1‑81 (MQ=255)               tttACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCt                                                              <  2:801978/80‑1 (MQ=255)                ttACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTg                                                             >  1:802021/1‑80 (MQ=255)                ttACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTg                                                             <  2:801979/80‑1 (MQ=255)                  aCTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGa                                                           <  2:801980/80‑1 (MQ=255)                   cTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACg                                                         >  1:802014/1‑81 (MQ=255)                       gAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGa                                                     <  2:801981/81‑1 (MQ=255)                           tGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGcc                                                 >  2:802018/1‑81 (MQ=255)                              cGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCt                                              <  1:801983/81‑1 (MQ=255)                                aGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTAc                                            <  1:801984/81‑1 (MQ=255)                                aGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTAc                                            >  2:802012/1‑81 (MQ=255)                                aGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTAc                                            >  1:802023/1‑81 (MQ=255)                                 gCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACg                                           >  1:802020/1‑81 (MQ=255)                                  cTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACgg                                          <  2:801984/81‑1 (MQ=255)                                    cATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGtt                                        >  1:802012/1‑81 (MQ=255)                                         tataATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCt                                   <  2:801985/81‑1 (MQ=255)                                         tataATGAATAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCt                                   >  2:802016/1‑81 (MQ=255)                                                  tAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCAc                          <  2:801986/81‑1 (MQ=255)                                                  tAAGGGTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCAc                          <  2:801987/81‑1 (MQ=255)                                                      ggTGTTTAAGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGc                      <  1:801989/81‑1 (MQ=255)                                                             aaGTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGt                >  2:802043/1‑80 (MQ=255)                                                               gTAAAGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTa               >  2:802023/1‑79 (MQ=255)                                                                 aaaGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGc           >  1:802032/1‑81 (MQ=255)                                                                  aaGGAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGc           <  1:801990/80‑1 (MQ=255)                                                                    ggAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTa        <  2:801990/81‑1 (MQ=255)                                                                     gAAAACATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTaa       <  2:801991/81‑1 (MQ=255)                                                                         aCATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGtt   <  1:801993/81‑1 (MQ=255)                                                                          cATTACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTc  <  2:801993/81‑1 (MQ=255)                                                                             |                                                                              CGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTC  >  NC_000913/4252222‑4252375

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