Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,602,506	T→A	D93E (GAT→GAA)	yjjU →	patatin family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,602,506	0	T	A	100.0%	111.5 / NA	29	D93E (GAT→GAA)	yjjU	patatin family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (17/12);  total (17/12)

CCTCTCGGCGTTTATCTGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGATCCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGCCACTGCAAGCCAGATGCCGTT  >  NZ_CP009273/4602425‑4602583                                                                                  |                                                                              ccTCTCGGCGTTTATCTGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTgcg                                                                       >  2:241248/1‑90 (MQ=255)   tctcGGCGTTTATCTGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCt                                                                     >  1:570196/1‑90 (MQ=255)          gTTTATCTGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCg                                                              >  1:367525/1‑90 (MQ=255)           tttATCTGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGt                                                             <  1:269818/90‑1 (MQ=255)           tttATCTGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGt                                                             <  2:502461/90‑1 (MQ=255)                 tGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTggagga                                                       <  1:122811/90‑1 (MQ=255)                 tGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTggagga                                                       <  2:326260/90‑1 (MQ=255)                  gCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGaa                                                      >  2:284959/1‑90 (MQ=255)                         gCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTAt                                               >  1:115822/1‑90 (MQ=255)                            cGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGa                                            >  2:227321/1‑90 (MQ=255)                                 aCGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCg                                       >  2:200643/1‑90 (MQ=255)                                   gcgcgcAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGAc                                     <  1:192663/90‑1 (MQ=255)                                   gcgcgcAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGAc                                     >  1:598093/1‑90 (MQ=255)                                       gcAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGActggc                                 >  2:438403/1‑90 (MQ=255)                                        cAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGActggct                                <  2:135299/90‑1 (MQ=255)                                                tcatGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAgg                        >  2:66526/1‑90 (MQ=255)                                                  atGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGcc                      <  1:554740/90‑1 (MQ=255)                                                  atGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGcc                      >  2:255323/1‑90 (MQ=255)                                                  atGCGCTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGcc                      <  2:34647/90‑1 (MQ=255)                                                       cTATACCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGCCACTGc                 <  1:126606/90‑1 (MQ=255)                                                         ataCCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGActggct                                >  1:556881/1‑73 (MQ=255)                                                         ataCCACAAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGActggct                                <  2:556881/73‑1 (MQ=255)                                                              acaAAACGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGCCACTGCAAGCCAg          >  2:202308/1‑90 (MQ=255)                                                                aaaaCGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGCCACTGCAAGCCAGAt        <  1:445460/90‑1 (MQ=255)                                                                 aaaCGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGGGGAGGCCACTGCAAGCCAGATg       >  1:407063/1‑90 (MQ=255)                                                                   aCGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGCCACTGCAAGCCAGATGcc     >  1:60572/1‑90 (MQ=255)                                                                   aCGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGGGTAGGCCCCTGCAAGCCAGATGcc     >  1:394911/1‑90 (MQ=255)                                                                   aCGCGAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGGGGAGGCCACTGCAAGCCAGATGcc     >  1:280320/1‑90 (MQ=255)                                                                      cgAATTTTTCGAACCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGCCACTGCAAGCCAGATGCCGtt  <  2:570196/90‑1 (MQ=255)                                                                                  |                                                                              CCTCTCGGCGTTTATCTGCAATCAGCCCGGTTACGCGCGCAAAGTCATCATGCGCTATACCACAAAACGCGAATTTTTCGATCCATTGCGCTTTGTCCGTGGAGGAAATCTTATCGATCTCGACTGGCTGGTGGAGGCCACTGCAAGCCAGATGCCGTT  >  NZ_CP009273/4602425‑4602583

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