Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	575,416	T→C	N135S (AAT→AGT)	ybjX ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	575,416	0	T	C	100.0%	182.6 / NA	52	N135S (AAT→AGT)	ybjX	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/19);  total (33/19)

GTCAGCTCAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTATTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGGTAACGTA  >  minE/575350‑575479                                                                   |                                                                gTCAGCTCAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCaa                                                             <  1:1359067/71‑1 (MQ=255) gTCAGCTCAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCaa                                                             <  1:373966/71‑1 (MQ=255) gTCAGCTCAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCaa                                                             <  1:2525539/71‑1 (MQ=255) gTCAGCTCAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCa                                                              <  1:110348/70‑1 (MQ=255)        cAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAg                                                      >  1:1278441/1‑71 (MQ=255)        cAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAg                                                      >  1:1990468/1‑71 (MQ=255)         aaGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGc                                                     >  1:815788/1‑71 (MQ=255)           gCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTc                                                               >  1:2408223/1‑59 (MQ=255)           gCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTg                                                   <  1:1325065/71‑1 (MQ=255)             gTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGaa                                                 >  1:899836/1‑71 (MQ=255)             gTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGaa                                                 >  1:2989753/1‑71 (MQ=255)              tGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGAAGCCCGGGGGTACTCAAATATAAGCTGAAt                                                >  1:2409443/1‑71 (MQ=255)               gAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCaaa                                                            >  1:698213/1‑58 (MQ=255)                aaCTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAAttc                                              >  1:2840477/1‑71 (MQ=255)                   tGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAAttctt                                            <  1:2581113/70‑1 (MQ=255)                    gCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTcc                                          >  1:2804996/1‑71 (MQ=255)                    gCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTcc                                          >  1:1803003/1‑71 (MQ=255)                       cGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGGCCGGGGGTACTCAAATATAAGCTGAAttctt                                            >  1:2491520/1‑66 (MQ=255)                       cGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTccgc                                        >  1:625564/1‑70 (MQ=255)                         ccGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAAttctt                                            <  1:2629986/64‑1 (MQ=255)                           gTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGaca                                   <  1:506017/71‑1 (MQ=255)                           gTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGa                                     >  1:1627787/1‑69 (MQ=255)                              tttGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGacata                                 >  1:708622/1‑70 (MQ=255)                                tGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCGGCCGacatacat                              >  1:20689/1‑71 (MQ=255)                                   cTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTcg                           <  1:110226/71‑1 (MQ=255)                                    ttCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTcgc                          >  1:2017588/1‑71 (MQ=255)                                    ttCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTcgc                          >  1:224441/1‑71 (MQ=255)                                     tCCAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTcgcg                         >  1:331157/1‑71 (MQ=255)                                      ccAGCTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTcgcg                         >  1:1153022/1‑70 (MQ=255)                                          cTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTccgccg                                      <  1:839828/53‑1 (MQ=255)                                          cTTCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTa                    <  1:387803/71‑1 (MQ=255)                                            tCGCCAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGaa                                                 >  1:1741702/1‑40 (MQ=255)                                                cAGTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCa               >  1:772174/1‑70 (MQ=255)                                                  gTTGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATaa            >  1:860175/1‑71 (MQ=255)                                                   ttGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTa                    <  1:1520543/62‑1 (MQ=255)                                                   ttGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCAt              <  1:722787/68‑1 (MQ=255)                                                    tGCAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTaa                   >  1:1359299/1‑62 (MQ=255)                                                      cAGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgt        <  1:1310129/71‑1 (MQ=255)                                                       aGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACAtt                             >  1:1528464/1‑49 (MQ=255)                                                       aGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACAtt                             >  1:2309779/1‑49 (MQ=255)                                                       aGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgta       <  1:2497929/71‑1 (MQ=255)                                                       aGCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgta       <  1:948551/71‑1 (MQ=255)                                                        gCCCGGGGGTACTCAAATATAAGCTGATTTCTTccgcc                                       >  1:2601272/1‑38 (MQ=38)                                                        gCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTc                            <  1:184987/49‑1 (MQ=255)                                                        gCCCGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgtaa      >  1:511250/1‑71 (MQ=255)                                                        gCCCGGGGGAACTCAAATATAAGCTGAATTCTTCCGCCGacata                                 >  1:2804441/1‑44 (MQ=255)                                                         cccGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGacataca                               >  1:2609955/1‑45 (MQ=255)                                                         cccGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTa                    >  1:2965086/1‑56 (MQ=255)                                                         cccGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgtaac     <  1:605172/71‑1 (MQ=255)                                                           cGGGGGTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgtaacg    >  1:2308687/1‑70 (MQ=255)                                                             ggggTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgtaacgta  <  1:1625149/70‑1 (MQ=255)                                                              gggTACTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGgtaac     >  1:991695/1‑66 (MQ=255)                                                                   |                                                                GTCAGCTCAAGCGTGAACTGCTCGCCGTTTTTGCCTTCCAGCTTCGCCAGTTGCAGCCCGGGGGTATTCAAATATAAGCTGAATTCTTCCGCCGACATACATTCGCGGAGTAACGCATAATGGTAACGTA  >  minE/575350‑575479

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: ‑