Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,666,180	T→A	9.1%	I193F (ATC→TTC)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,666,180	0	T	A	9.1%	178.8 / 9.4	66	I193F (ATC→TTC)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base A (2/4);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCAC  >  minE/1666110‑1666246                                                                       |                                                                   ccTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGAt                                                                    >  1:1307125/1‑71 (MQ=255) ccTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGAt                                                                    <  1:3281615/71‑1 (MQ=255)   ttGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGAt                                                                    >  1:2778599/1‑69 (MQ=255)    tGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGAt                                                                    >  1:3440373/1‑68 (MQ=255)      cTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAgg                                                                <  1:507401/70‑1 (MQ=255)       tCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTg                                                              >  1:2588593/1‑71 (MQ=255)         gTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTgg                                                             >  1:1754340/1‑70 (MQ=255)          tCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAg                                                                 >  1:697054/1‑65 (MQ=255)          tCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGaa                                                           <  1:1868675/71‑1 (MQ=255)           cTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAAt                                                          >  1:896378/1‑71 (MQ=255)              tACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTc                                                        <  1:996921/70‑1 (MQ=255)              tACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCAGATCAGGTGGAATTc                                                        <  1:3057426/70‑1 (MQ=255)                 aGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtca                                                    <  1:3742939/71‑1 (MQ=255)                   gCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtcatc                                                  <  1:3244218/71‑1 (MQ=255)                    cGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtcatca                                                 <  1:3278491/71‑1 (MQ=255)                     gTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCaa                                                >  1:1624327/1‑71 (MQ=255)                        cGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtcatc                                                  >  1:2381107/1‑66 (MQ=255)                        cGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAacca                                             >  1:739063/1‑71 (MQ=255)                        cGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAacc                                              <  1:3012158/70‑1 (MQ=255)                        cGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAac                                               >  1:1036908/1‑69 (MQ=255)                         ggTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAaccac                                            >  1:3439994/1‑71 (MQ=255)                          gTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGAACAGGTGGaa                                                           <  1:483521/55‑1 (MQ=255)                          gTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGAACAGGTGGaa                                                           <  1:3494879/55‑1 (MQ=255)                          gTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGAACAGGTGGaa                                                           <  1:2459660/55‑1 (MQ=255)                           tttACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTg                                                              <  1:3399490/51‑1 (MQ=255)                             tACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAAccacca                                          >  1:600583/1‑69 (MQ=255)                             tACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAg                                         <  1:1089457/70‑1 (MQ=255)                               cGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAgcgc                                      <  1:15440/71‑1 (MQ=255)                               cGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAgcg                                       <  1:2199647/70‑1 (MQ=255)                               cGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAgcg                                       <  1:3457410/70‑1 (MQ=255)                                gTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCt                                     <  1:2224588/71‑1 (MQ=255)                                   ccGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCtttt                                  >  1:2376289/1‑71 (MQ=255)                                   ccGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCttt                                   <  1:1141018/70‑1 (MQ=255)                                    cGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGTATTCGTCATCAACCACCAGCGCttttt                                 <  1:2413306/71‑1 (MQ=255)                                    cGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCtttt                                  <  1:291876/70‑1 (MQ=255)                                     gCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCtt                                    <  1:1947877/67‑1 (MQ=255)                                      cTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTc                                >  1:1329292/1‑70 (MQ=255)                                      cTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTATTc                                >  1:2904351/1‑70 (MQ=255)                                       tttCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCt                                     <  1:152952/64‑1 (MQ=255)                                         tCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCg                                                       <  1:797148/44‑1 (MQ=255)                                             gAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGtcat                                                   <  1:2400291/44‑1 (MQ=255)                                              aaGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGttt                       >  1:2776017/1‑71 (MQ=255)                                              aaGTCTTTCACGGTGATCATGCCGAACAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGttt                       >  1:2450036/1‑71 (MQ=255)                                               aGTCTTTCACGGTGATCATGCTGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGtttt                      <  1:273310/71‑1 (MQ=255)                                                gTCTTTCACGGTGATCATGCCGATCAGGTGGAATTc                                                        >  1:2993011/1‑36 (MQ=255)                                                   tttCACGGTGATCATGCCGATCAGGTGGAATTCGtcatc                                                  <  1:2358107/39‑1 (MQ=255)                                                   tttCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGt                  >  1:2659289/1‑71 (MQ=255)                                                    ttCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCttt                                   <  1:796128/53‑1 (MQ=255)                                                    ttCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCttt                                   <  1:1019100/53‑1 (MQ=255)                                                    ttCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCttt                                   <  1:263323/53‑1 (MQ=255)                                                     tCACGGTGATCATGCCGATCAGGTGGAATTCGtcat                                                   >  1:135148/1‑36 (MQ=255)                                                     tCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAac                                               >  1:1118916/1‑40 (MQ=255)                                                     tCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCa                               <  1:931553/56‑1 (MQ=255)                                                     tCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTg                 >  1:2207649/1‑70 (MQ=255)                                                     tCACGGTGATCATGCCGAACAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTg                 >  1:1328797/1‑70 (MQ=255)                                                       aCGGTGATCATGCCGATCAGGTGGAATTCGTCATCaa                                                >  1:1594029/1‑37 (MQ=255)                                                       aCGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAt              >  1:3253902/1‑71 (MQ=255)                                                       aCGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAt              <  1:3125968/71‑1 (MQ=255)                                                         ggTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAttt            <  1:1884560/71‑1 (MQ=255)                                                         ggTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAttt            <  1:3680269/71‑1 (MQ=255)                                                         ggTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAttt            <  1:722884/71‑1 (MQ=255)                                                         ggTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAttt            <  1:2352960/71‑1 (MQ=255)                                                         ggTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAttt            <  1:2273874/71‑1 (MQ=255)                                                          gTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAtttt           >  1:3021120/1‑71 (MQ=255)                                                          gTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAtttt           >  1:343870/1‑71 (MQ=255)                                                             aTCATGCCGATCAGGTGGATTTCGTCATCAACCACCAgc                                        >  1:3863264/1‑39 (MQ=255)                                                                 tGCCGAACAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGc    <  1:2495703/71‑1 (MQ=255)                                                                  gCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTc                    >  1:3201433/1‑54 (MQ=255)                                                                   ccGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGcac  >  1:2170687/1‑71 (MQ=255)                                                                     gATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCAtttt           >  1:148397/1‑60 (MQ=255)                                                                       |                                                                   CCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGCCAGCAC  >  minE/1666110‑1666246

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