Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,680,754	T→C	Y1407C (TAC→TGC)	yfhM ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,680,754	0	T	C	100.0%	191.5 / NA	55	Y1407C (TAC→TGC)	yfhM	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/27);  total (28/27)

CGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTG  >  minE/1680686‑1680820                                                                     |                                                                   cgTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCAc                                                                  >  1:1755461/1‑71 (MQ=255)     aCTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGc                                                               <  1:316139/70‑1 (MQ=255)     aCTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGc                                                               <  1:510374/70‑1 (MQ=255)           cTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCgg                                                        <  1:169175/71‑1 (MQ=255)           cTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCgg                                                        <  1:81088/71‑1 (MQ=255)           cTGCGTCGATAG‑CAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCgg                                                        <  1:1615861/70‑1 (MQ=255)             gCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTTCt                                                              <  1:2335372/63‑1 (MQ=255)              cGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCggg                                                       >  1:2626144/1‑69 (MQ=255)                  gATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCa                                                  <  1:167114/70‑1 (MQ=255)                  gATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAg                                                 >  1:2706730/1‑71 (MQ=255)                   aTAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGt                                                <  1:1348099/71‑1 (MQ=255)                    tAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGtt                                               >  1:228304/1‑71 (MQ=255)                    tAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGtt                                               <  1:1427573/71‑1 (MQ=255)                       ccAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAt                                            >  1:2895643/1‑71 (MQ=255)                        cAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGttt                                              >  1:404828/1‑68 (MQ=255)                        cAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtt                                           <  1:1242606/71‑1 (MQ=255)                        cAGCGTCCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtt                                           <  1:2687354/71‑1 (MQ=255)                         aGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCAc                                                                  >  1:2330434/1‑47 (MQ=255)                           cGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTttc                                        >  1:3014331/1‑71 (MQ=255)                             ttCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCAc                                                                  >  1:1170278/1‑43 (MQ=255)                             ttCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCAc                                                                  >  1:1000669/1‑43 (MQ=255)                             ttCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCa                                                  <  1:1877388/59‑1 (MQ=255)                             ttCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtt                                           <  1:1955337/66‑1 (MQ=255)                               cACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTttcttc                                     >  1:992296/1‑70 (MQ=255)                               cACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTttcttc                                     >  1:1223497/1‑70 (MQ=255)                                  cAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGc                                 <  1:2867609/71‑1 (MQ=255)                                   aaaCGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTttct                                       <  1:501350/64‑1 (MQ=255)                                    aaCGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGc                                 >  1:664069/1‑69 (MQ=255)                                      cGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAgg                             <  1:1848582/71‑1 (MQ=255)                                         ctgctgCGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGga                            >  1:259229/1‑69 (MQ=255)                                           gctgcGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGagag                         >  1:1264614/1‑70 (MQ=255)                                            ctgcGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTttctt                                      <  1:1400986/56‑1 (MQ=255)                                            ctgcGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTttctt                                      <  1:1450640/56‑1 (MQ=255)                                             tgcGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCAt                      >  1:992295/1‑71 (MQ=255)                                               cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCAtt                     <  1:1313381/70‑1 (MQ=255)                                               cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCAtt                     >  1:2004886/1‑70 (MQ=255)                                               cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCAtt                     <  1:1147007/70‑1 (MQ=255)                                               cGAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGAAAGGAGAGCAtt                     >  1:2605932/1‑70 (MQ=255)                                                gAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTAtttt                                         >  1:1516968/1‑49 (MQ=255)                                                gAAAGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTaa                   <  1:2160648/71‑1 (MQ=255)                                                 aaaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAc                  >  1:148533/1‑71 (MQ=255)                                                  aaGTGTGTTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAAcg                 >  1:2313167/1‑71 (MQ=255)                                                   aGTGTGTTCAGCAAAGTGCACTGCTCATCGTGTAGCAGTTTATTTTCTTCCAGCaa                               >  1:1150422/1‑56 (MQ=255)                                                    gtgtgtTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCAtt                     >  1:2158135/1‑65 (MQ=255)                                                     tgtgtTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGt              >  1:855848/1‑71 (MQ=255)                                                      gtgtTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGtt             >  1:2922347/1‑71 (MQ=255)                                                       tgtTCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGtt             <  1:1651121/70‑1 (MQ=255)                                                          tCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTttcttaa                                    >  1:151562/1‑42 (MQ=39)                                                          tCAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTcg         >  1:70655/1‑71 (MQ=255)                                                           cAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTa                    <  1:2524980/59‑1 (MQ=255)                                                           cAGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTcgc        >  1:1111586/1‑71 (MQ=255)                                                            aGCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTcgcg       <  1:1068288/71‑1 (MQ=255)                                                             gCAAAGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCgcgc      <  1:1709852/71‑1 (MQ=255)                                                                 aGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAg                              <  1:2950104/43‑1 (MQ=255)                                                                 aGTGCACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTg  <  1:2460826/71‑1 (MQ=255)                                                                     |                                                                   CGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTG  >  minE/1680686‑1680820

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