Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	406,585	A→G	D8D (GAT→GAC)	ybeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	406,585	0	A	G	100.0%	104.9 / NA	29	D8D (GAT→GAC)	ybeB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/17);  total (12/17)

CGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTC  >  minE/406530‑406650                                                        |                                                                  cGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTa                                                    <  1:2160688/71‑1 (MQ=255) cGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTa                                                    <  1:1482042/71‑1 (MQ=255)      aaGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAg                                                           <  1:1196793/59‑1 (MQ=255)      aaGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAg                                                           <  1:1734368/59‑1 (MQ=255)       aGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTGACCCTg                                               >  1:928467/1‑70 (MQ=255)           gatgatGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGt                                          >  1:1082839/1‑71 (MQ=255)            atgatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTg                                                              <  1:1433140/50‑1 (MQ=255)            atgatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTg                                                              <  1:2732011/50‑1 (MQ=255)            atgatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTg                                                              <  1:1900454/50‑1 (MQ=255)            atgatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGt                                          <  1:3196803/70‑1 (MQ=255)             tgatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAg                                                           <  1:1475486/52‑1 (MQ=255)                     tGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTg                                >  1:2043485/1‑71 (MQ=255)                       aCCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTggg                              <  1:2132027/71‑1 (MQ=255)                       aCCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTgg                               <  1:261339/70‑1 (MQ=255)                           ttGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGaa                          >  1:2577896/1‑71 (MQ=255)                             gAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAAt                         >  1:865320/1‑70 (MQ=255)                                 tcatcaATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaaa                     <  1:1671771/70‑1 (MQ=255)                                    tcaATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGa                 <  1:2158589/71‑1 (MQ=255)                                     caATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGa                 >  1:1948165/1‑70 (MQ=255)                                      aaTTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATa               <  1:2262328/71‑1 (MQ=255)                                           tGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATaaaaa           <  1:1239365/70‑1 (MQ=255)                                            gTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTcccc                                   >  1:2526121/1‑45 (MQ=255)                                            gTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTcccc                                   >  1:1451405/1‑45 (MQ=255)                                            gTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATg         >  1:76870/1‑71 (MQ=255)                                               gATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGt       >  1:662771/1‑70 (MQ=255)                                                  aaCAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGt   >  1:2577953/1‑71 (MQ=255)                                                  aaCAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGt   >  1:745790/1‑71 (MQ=255)                                                   aCAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGt   <  1:3098895/70‑1 (MQ=255)                                                   aCAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTc  <  1:3103587/71‑1 (MQ=255)                                                     aaagTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTggg                              <  1:161043/41‑1 (MQ=255)                                                     aaagTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTggg                              <  1:1429699/41‑1 (MQ=255)                                                        |                                                                  CGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTC  >  minE/406530‑406650

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