Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,899,378	T→C	T6A (ACC→GCC)	yieK ←	putative glucosamine‑6‑phosphate deaminase YieK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,899,378	0	T	C	100.0%	109.8 / NA	32	T6A (ACC→GCC)	yieK	putative glucosamine‑6‑phosphate deaminase YieK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/13);  total (19/13)

GTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGA  >  NC_000913/3899301‑3899454                                                                              |                                                                             gTTAACACGACGCGTCTTAGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAAt                                                                           >  1:706793/1‑81 (MQ=255)    aaCACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGAt                                                                        >  1:706766/1‑81 (MQ=255)      cacGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTa                                                                      <  2:706725/81‑1 (MQ=255)          acgCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAAttt                                                                  >  1:706727/1‑81 (MQ=255)              gTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATc                                                              <  2:706727/81‑1 (MQ=255)              gTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATc                                                              >  1:706765/1‑81 (MQ=255)                cTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCt                                                             >  2:706769/1‑80 (MQ=255)                   cGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCtt                                                            <  2:706729/78‑1 (MQ=255)                   cGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTa                                                           <  2:706728/79‑1 (MQ=255)                    gACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATg                                                        <  1:706731/81‑1 (MQ=255)                    gACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATg                                                        >  2:706761/1‑81 (MQ=255)                    gACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTACCTGGTAATCTTCGGCAATGATTAATTTCATCTTAAt                                                         <  1:706732/80‑1 (MQ=255)                       atatAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCt                                                     <  1:706733/81‑1 (MQ=255)                         ataACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCtta                                                   >  1:706763/1‑81 (MQ=255)                               aaGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCttattttat                                             <  1:706734/81‑1 (MQ=255)                                 gCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTtattttattt                                           <  2:706734/81‑1 (MQ=255)                                             ccGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTg                               <  2:706735/81‑1 (MQ=255)                                               gcgACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTa                             >  1:706780/1‑81 (MQ=255)                                                   cacGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACaaa                         <  1:706737/81‑1 (MQ=255)                                                      ggCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGc                      >  2:706763/1‑81 (MQ=255)                                                       gCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCt                     >  1:706779/1‑81 (MQ=255)                                                          cATTTCCTGGTAATCTTCAGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCa                  >  2:706777/1‑81 (MQ=255)                                                           aTTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGc                   <  2:706737/79‑1 (MQ=255)                                                            tttCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCa                  >  1:706786/1‑79 (MQ=255)                                                             ttCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAg               >  1:706794/1‑81 (MQ=255)                                                                   gTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGt         >  1:706781/1‑81 (MQ=255)                                                                      aTCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGc      >  1:706783/1‑81 (MQ=255)                                                                        cTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGc      >  2:706766/1‑79 (MQ=255)                                                                         ttCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGc      >  1:706769/1‑78 (MQ=255)                                                                           cGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGcc     <  1:706739/77‑1 (MQ=255)                                                                           cGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGa  >  2:706772/1‑80 (MQ=255)                                                                           cGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGa  >  2:706773/1‑80 (MQ=255)                                                                              |                                                                             GTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGA  >  NC_000913/3899301‑3899454

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