Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	857,377	T→C	intergenic (‑8/+4)	trpB ← / ← trpC	tryptophan synthase, beta subunit/fused indole‑3‑glycerolphosphate synthetase and N‑(5‑phosphoribosyl)anthranilate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	857,377	0	T	C	100.0%	253.7 / NA	70	intergenic (‑8/+4)	trpB/trpC	tryptophan synthase, beta subunit/fused indole‑3‑glycerolphosphate synthetase and N‑(5‑phosphoribosyl)anthranilate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/40);  total (30/40)

GATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGC  >  minE/857310‑857444                                                                    |                                                                    gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGttcctctcc                                                                  >  1:3240084/1‑71 (MQ=255) gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGttcctctcc                                                                  >  1:3681742/1‑71 (MQ=255) gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGttcctctcc                                                                  >  1:2384344/1‑71 (MQ=255) gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGttcctctcc                                                                  >  1:650389/1‑71 (MQ=255) gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGttcctctcc                                                                  >  1:1844314/1‑71 (MQ=255) gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGttcctctc                                                                   >  1:138809/1‑70 (MQ=255) gATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGttcctctc                                                                   >  1:719029/1‑70 (MQ=255)         gCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTaa                                                              <  1:1344393/67‑1 (MQ=255)              tACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATgcgcg                                                      >  1:3335235/1‑70 (MQ=255)              tACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATgcgc                                                       <  1:1840027/69‑1 (MQ=255)               aCATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATgcgcg                                                      >  1:2906811/1‑69 (MQ=255)                cATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCa                                                    <  1:1677797/70‑1 (MQ=255)                 aTGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTaa                                                              <  1:2347979/59‑1 (MQ=255)                   gccgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTa                                                               <  1:1866001/56‑1 (MQ=255)                   gccgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTa                                                               <  1:1618661/56‑1 (MQ=255)                   gccgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATgc                                                         >  1:1065051/1‑62 (MQ=255)                    ccgccAAACTCACCTAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGt                                                <  1:943535/70‑1 (MQ=255)                    ccgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCa                                                    <  1:4104221/66‑1 (MQ=255)                    ccgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGt                                                <  1:1652190/70‑1 (MQ=255)                    ccgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGt                                                <  1:3292710/70‑1 (MQ=255)                     cgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAat                                                             >  1:2668107/1‑56 (MQ=255)                     cgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAat                                                             >  1:3059501/1‑56 (MQ=255)                     cgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAat                                                             <  1:238427/56‑1 (MQ=255)                     cgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGc                                                  >  1:769112/1‑67 (MQ=255)                     cgccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCt                                              <  1:904278/71‑1 (MQ=255)                      gccAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGc                                                  <  1:2045381/66‑1 (MQ=255)                        cAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAata                                                            >  1:2588876/1‑54 (MQ=255)                              cACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAAcc                                      >  1:1033354/1‑70 (MQ=255)                               aCCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCTCAGCGt                                                >  1:3206648/1‑59 (MQ=255)                               aCCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCa                                                    <  1:3616591/55‑1 (MQ=255)                               aCCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAAcc                                      <  1:2171893/69‑1 (MQ=255)                               aCCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCg                                     <  1:4076824/70‑1 (MQ=255)                               aCCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCg                                     <  1:3560020/70‑1 (MQ=255)                               aCCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCg                                     <  1:3107566/70‑1 (MQ=255)                               aCCAAAACAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGa                                    <  1:3225483/71‑1 (MQ=255)                                ccAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTTGAAAACCGAg                                   <  1:183393/71‑1 (MQ=255)                                ccAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAg                                   <  1:2905923/71‑1 (MQ=255)                                  aaaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGc                                 >  1:651485/1‑71 (MQ=255)                                   aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATgc                                                         <  1:174329/46‑1 (MQ=255)                                   aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCa                                                    >  1:111613/1‑51 (MQ=255)                                   aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAg                                                   >  1:192260/1‑52 (MQ=255)                                   aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGcc                                <  1:1869486/71‑1 (MQ=255)                                   aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGc                                 >  1:2414053/1‑70 (MQ=255)                                   aaaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGc                                 >  1:3519084/1‑70 (MQ=255)                                    aaTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGa                                    >  1:3025732/1‑66 (MQ=255)                                     aTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGa                                    >  1:3625814/1‑65 (MQ=255)                                     aTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGa                                    >  1:2854783/1‑65 (MQ=255)                                     aTAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCaa                              <  1:2911588/71‑1 (MQ=255)                                      tAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAata                                                            >  1:4071018/1‑40 (MQ=255)                                      tAGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAata                                                            >  1:2947187/1‑40 (MQ=255)                                       aGGGGTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGa                                    >  1:859996/1‑63 (MQ=255)                                         gggTTAATTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCa                                                    <  1:1484593/45‑1 (MQ=255)                                          ggTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGAc                         <  1:140980/71‑1 (MQ=255)                                           gTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCg                                     <  1:1247663/58‑1 (MQ=255)                                           gTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCg                                     <  1:1067251/58‑1 (MQ=255)                                           gTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGAc                         >  1:2869190/1‑70 (MQ=255)                                           gTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGAc                         <  1:177335/70‑1 (MQ=255)                                           gTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACg                        >  1:3921184/1‑71 (MQ=255)                                           gTTAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACg                        >  1:2137039/1‑71 (MQ=255)                                            ttAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGAc                         <  1:1358497/69‑1 (MQ=255)                                            ttAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGt                       >  1:978633/1‑71 (MQ=255)                                            ttAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGt                       <  1:3578575/71‑1 (MQ=255)                                            ttAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGt                       <  1:2501023/71‑1 (MQ=255)                                            ttAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGt                       >  1:3417270/1‑71 (MQ=255)                                             tAAGTAATGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTg                      >  1:570217/1‑71 (MQ=255)                                                  aaTGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAg                                   <  1:3295784/53‑1 (MQ=255)                                                  aaTGTTGTCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCt                 >  1:402874/1‑71 (MQ=255)                                                    tgttgtCATTGTTCGTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCttt               >  1:3866811/1‑71 (MQ=255)                                                    tgttgtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCttt               <  1:2738642/71‑1 (MQ=255)                                                      ttgtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTgg                                            <  1:2175179/40‑1 (MQ=255)                                                      ttgtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTg              <  1:1812151/70‑1 (MQ=255)                                                      ttgtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTg              <  1:1807375/70‑1 (MQ=255)                                                        gtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGAc                         >  1:1850755/1‑57 (MQ=255)                                                        gtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCttt               <  1:3300237/67‑1 (MQ=255)                                                        gtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGAt            <  1:1530235/70‑1 (MQ=255)                                                        gtCATTGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGAt            <  1:3251349/70‑1 (MQ=255)                                                             tGTTCCTCTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCgg      >  1:3081331/1‑71 (MQ=255)                                                                 cctctccttAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGc  >  1:1170559/1‑71 (MQ=255)                                                                    |                                                                    GATTTGTGGCACGTACATGCCGCCAAACTCACCAAAATAGGGGTTAAGTAATGTTGTCATTGTTCCTTTCCTTAATATGCGCGCAGCGTCTGGAAAACCGAGGCCAAAAGACGTGCGTCTTTGATGCCCGGTTGC  >  minE/857310‑857444

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