Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,537,453	C→A	Q38K (CAG→AAG)	yfeK →	protein YfeK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,537,453	0	C	A	100.0%	91.5 / NA	26	Q38K (CAG→AAG)	yfeK	protein YfeK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/10);  total (16/10)

AATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTG  >  NC_000913/2537377‑2537527                                                                             |                                                                           aaTGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGa                                                                         >  1:450986/1‑80 (MQ=255)    gACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGaaa                                                                       >  1:450982/1‑79 (MQ=255)           cccGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGAttt                                                               >  1:450978/1‑80 (MQ=255)                ttACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATAtt                                                         >  1:450974/1‑81 (MQ=255)                 tACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATAtt                                                         <  2:450948/80‑1 (MQ=255)                      gAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTgcg                                                   <  1:450949/81‑1 (MQ=255)                      gAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTgcg                                                   >  1:450980/1‑81 (MQ=255)                         gTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCaa                                                <  2:450950/81‑1 (MQ=255)                           tAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACg                                              >  1:450996/1‑81 (MQ=255)                           tAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACg                                              >  1:450979/1‑81 (MQ=255)                              cTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGtg                                           <  1:450951/81‑1 (MQ=255)                              cTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGt                                            >  2:450985/1‑80 (MQ=255)                                 cccATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGtgat                                         <  1:450952/80‑1 (MQ=255)                                        agaagCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATAc                                 >  1:450989/1‑81 (MQ=255)                                         gaagCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATAcc                                <  2:450953/81‑1 (MQ=255)                                          aagCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCt                               >  1:450976/1‑81 (MQ=255)                                          aagCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCt                               >  2:450992/1‑81 (MQ=255)                                           agCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTg                              >  1:450984/1‑81 (MQ=255)                                            gCCCGCATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGc                             >  1:451000/1‑81 (MQ=255)                                                 cATCAACGCCATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGa                        <  1:450954/81‑1 (MQ=255)                                                         ccATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGttt                <  1:450955/81‑1 (MQ=255)                                                          cATGCTGGAGGGATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTtc               >  1:450993/1‑81 (MQ=255)                                                                    ggATTAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGt     <  1:450956/81‑1 (MQ=255)                                                                       ttAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTg  <  2:450957/81‑1 (MQ=255)                                                                        tAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTg  >  1:450990/1‑80 (MQ=255)                                                                        tAGCAAAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTg  >  2:450991/1‑80 (MQ=255)                                                                             |                                                                           AATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTG  >  NC_000913/2537377‑2537527

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