Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,666,716	A→G	12.8%	V14A (GTT→GCT)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,666,716	0	A	G	12.8%	221.6 / 21.0	86	V14A (GTT→GCT)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base G (6/5);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

GTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGG  >  minE/1666647‑1666779                                                                      |                                                                gTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                                                                 <  1:3284539/70‑1 (MQ=255) gTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                                                                 <  1:3288980/70‑1 (MQ=255)  ttttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcg                                                               >  1:2798672/1‑71 (MQ=255)    ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtc                                                             >  1:1791440/1‑71 (MQ=255)    ttCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtc                                                             >  1:6743/1‑71 (MQ=255)      cGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcgt                                                           >  1:230087/1‑71 (MQ=255)          aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCaaa                                                       >  1:2868694/1‑71 (MQ=255)          aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                                                       >  1:2046044/1‑71 (MQ=255)          aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                                                       >  1:1394026/1‑71 (MQ=255)          aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                                                       >  1:3052031/1‑71 (MQ=255)          aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                        >  1:1756643/1‑70 (MQ=255)          aGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                        <  1:3094256/70‑1 (MQ=255)             tGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtc                                                             >  1:1659791/1‑62 (MQ=255)             tGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                    >  1:1273855/1‑71 (MQ=255)             tGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                    >  1:1242140/1‑71 (MQ=255)             tGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                    >  1:3438673/1‑71 (MQ=255)              gggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGt                                                    >  1:2829962/1‑70 (MQ=255)               ggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGt                                                    >  1:3436748/1‑69 (MQ=255)               ggTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCa                                                  <  1:3020097/71‑1 (MQ=255)                gTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCag                                                 <  1:2891029/71‑1 (MQ=255)                gTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCag                                                 <  1:875610/71‑1 (MQ=255)                   cTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGc                                              >  1:3554117/1‑71 (MQ=255)                   cTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGc                                              >  1:75410/1‑71 (MQ=255)                      aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGcttc                                           >  1:2406564/1‑71 (MQ=255)                      aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAAc                                                      >  1:149477/1‑60 (MQ=255)                      aGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCaga                                                >  1:1746754/1‑66 (MQ=255)                      aGGTCAGCAGTATTAGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGctt                                            <  1:1515430/70‑1 (MQ=255)                           agcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTa                                       >  1:3519533/1‑70 (MQ=255)                           agcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAg                                      >  1:534596/1‑71 (MQ=255)                           agcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcg                                                            <  1:2571234/49‑1 (MQ=255)                            gcagTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGc                                     <  1:2638736/71‑1 (MQ=255)                               gTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                   >  1:2994311/1‑70 (MQ=255)                               gTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                   >  1:2310868/1‑70 (MQ=255)                               gTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGa                                   >  1:1735611/1‑70 (MQ=255)                                tATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATa                                 >  1:2958759/1‑71 (MQ=255)                                  ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCttt                                        <  1:754534/62‑1 (MQ=255)                                  ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAg                                      <  1:1378269/64‑1 (MQ=255)                                  ttCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                               >  1:1224883/1‑71 (MQ=255)                                    cGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGct                                             >  1:2270758/1‑55 (MQ=255)                                    cGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                               >  1:3802386/1‑69 (MQ=255)                                    cGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                               >  1:2453133/1‑69 (MQ=255)                                     ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                            <  1:3097739/71‑1 (MQ=255)                                     ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                            <  1:3249499/71‑1 (MQ=255)                                     ggCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAg                            <  1:1882657/71‑1 (MQ=255)                                      gCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTa                                       <  1:3307880/59‑1 (MQ=255)                                      gCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGc                           <  1:285192/71‑1 (MQ=255)                                         gAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCaa                                                        >  1:2457951/1‑39 (MQ=255)                                         gAACGGTAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATg                        <  1:247278/71‑1 (MQ=255)                                         gAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAAc                                                      >  1:1191941/1‑41 (MQ=255)                                          aaCGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGc                                     >  1:1007246/1‑57 (MQ=255)                                          aaCGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATgg                       >  1:3273720/1‑71 (MQ=255)                                           aCGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                                                        <  1:1128277/37‑1 (MQ=255)                                            cGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATggg                      >  1:3426135/1‑70 (MQ=255)                                              gTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCaa                   <  1:2370512/71‑1 (MQ=255)                                               tAGAGTGAGCAGGAACGAGGAGAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAat                  <  1:2112447/71‑1 (MQ=255)                                                 gagTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCaga                                                >  1:3721634/1‑39 (MQ=255)                                                 gagTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAatat                >  1:2810280/1‑71 (MQ=255)                                                  agTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAat                  <  1:1214278/68‑1 (MQ=255)                                                    tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATAc                                <  1:3757001/52‑1 (MQ=255)                                                    tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATAc                                <  1:700385/52‑1 (MQ=255)                                                    tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATAc                                <  1:2205072/52‑1 (MQ=255)                                                    tGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtctc             <  1:30148/71‑1 (MQ=255)                                                      aGCAGGAACGAGGAGTACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtctc             <  1:942095/69‑1 (MQ=255)                                                      aGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa           <  1:1351616/71‑1 (MQ=255)                                                      aGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa           <  1:314118/71‑1 (MQ=255)                                                      aGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa           >  1:3659705/1‑71 (MQ=255)                                                       gCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa           <  1:3579649/70‑1 (MQ=255)                                                       gCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGa           <  1:2295946/70‑1 (MQ=255)                                                        cAGGAACGAGGGGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtc               >  1:2169834/1‑65 (MQ=255)                                                        cAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATAtc               >  1:1022467/1‑65 (MQ=255)                                                         aGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGc                                     <  1:3020500/42‑1 (MQ=255)                                                         aGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCg                                    >  1:2953220/1‑43 (MQ=255)                                                         aGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCa        <  1:1099526/71‑1 (MQ=255)                                                          ggAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCa                          >  1:3460017/1‑52 (MQ=255)                                                           gAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACg                               <  1:2532104/46‑1 (MQ=255)                                                           gAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCaga      <  1:1104507/71‑1 (MQ=255)                                                            aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag     <  1:2958528/71‑1 (MQ=255)                                                            aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag     <  1:3040165/71‑1 (MQ=255)                                                            aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag     <  1:2134692/71‑1 (MQ=255)                                                            aaCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCagag     <  1:1434370/71‑1 (MQ=255)                                                             aCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTa                             >  1:3712982/1‑46 (MQ=255)                                                             aCGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTa                             >  1:2041612/1‑46 (MQ=255)                                                              cGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCGTGGGCAATATCTCGACCAGAGTg   <  1:3061678/71‑1 (MQ=255)                                                               gaggagAGCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTgg  <  1:2524806/71‑1 (MQ=255)                                                               gaggagAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTgg  >  1:590842/1‑71 (MQ=255)                                                                  gagaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCaa                   <  1:579746/51‑1 (MQ=255)                                                                  gagaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCaa                   <  1:1524914/51‑1 (MQ=255)                                                                   agaACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGAc          >  1:1858928/1‑59 (MQ=255)                                                                     aaCGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTgg  <  1:1044575/65‑1 (MQ=255)                                                                      |                                                                GTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGG  >  minE/1666647‑1666779

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