Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,429,412	T→C	N59N (AAT→AAC)	yrdA →	protein YrdA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,429,412	0	T	C	100.0%	105.4 / NA	31	N59N (AAT→AAC)	yrdA	protein YrdA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/12);  total (19/12)

GGCTGATGATGTGGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATT  >  NC_000913/3429334‑3429477                                                                               |                                                                  ggCTGATGATGTGGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAacat                                                                 >  2:586892/1‑81 (MQ=255)    tgatgaTGTGGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATc                                                                >  1:586890/1‑79 (MQ=255)      atgatgTGGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAgg                                                            <  1:586843/81‑1 (MQ=255)         atgTGGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATg                                                         >  2:586875/1‑81 (MQ=255)         atgTGGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATg                                                         >  1:586871/1‑81 (MQ=255)           gtgGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATgg                                                        >  1:586884/1‑80 (MQ=255)            tgGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGc                                                       <  1:586844/80‑1 (MQ=255)               ggATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTa                                                   <  1:586845/81‑1 (MQ=255)                   cTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGtt                                               >  2:586886/1‑81 (MQ=255)                   cTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGtt                                               >  2:586874/1‑81 (MQ=255)                      gCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGt                                                <  1:586846/77‑1 (MQ=255)                           tCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCa                                       >  2:586882/1‑81 (MQ=255)                             gTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACt                                     >  2:586879/1‑81 (MQ=255)                             gTTGTGATTCGTGGAGATGTACATTATGTACAGATCGAAGCACGCACCAACATCCAGGATGACAGTATGTTGCATGTCa                                       <  2:586846/79‑1 (MQ=255)                                 tgATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCAt                                  >  2:586899/1‑80 (MQ=255)                                 tgATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATa                                 <  1:586848/81‑1 (MQ=255)                                      cGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATcc                            >  2:586894/1‑81 (MQ=255)                                        tGGAGATGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTc                          <  2:586848/81‑1 (MQ=255)                                          gagaTGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCg                         >  2:586895/1‑80 (MQ=255)                                           agaTGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTa                       <  1:586850/81‑1 (MQ=255)                                             aTGTACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTAc                      >  1:586898/1‑80 (MQ=255)                                                tACATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAAc                   <  2:586850/80‑1 (MQ=255)                                                 aCATTATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAAcc                  >  2:586884/1‑80 (MQ=255)                                                    ttATGTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGa              <  1:586852/81‑1 (MQ=255)                                                        gTACAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATgg           >  1:586909/1‑80 (MQ=255)                                                           cAGATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAAc       >  2:586893/1‑81 (MQ=255)                                                             gATCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAAcc      >  1:586899/1‑80 (MQ=255)                                                             gATCGGAGCACGCACCAACATCCAGGATGACAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAAccc     <  1:586853/81‑1 (MQ=255)                                                              aTCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCa    >  1:586901/1‑81 (MQ=255)                                                               tCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCAt   <  1:586854/81‑1 (MQ=255)                                                               tCGGAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCAt   >  2:586904/1‑81 (MQ=255)                                                                 ggAGCACGCACCAACATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCagt  >  2:586907/1‑78 (MQ=255)                                                                               |                                                                  GGCTGATGATGTGGGGATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATT  >  NC_000913/3429334‑3429477

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