Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,400,915	A→C	51.7%	A172A (GCT→GCG)	yigB ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,400,915	0	A	C	51.7%	10.3 / 99.9	95	A172A (GCT→GCG)	yigB	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (28/18);  new base C (29/20);  total (57/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCACGTCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAGCGCG  >  minE/2400848‑2400968                                                                    |                                                      ccacGTCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCa                                                       >  1:303112/1‑68 (MQ=255) ccacGTCAGTGGTGAGATCGT‑CCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAg                                                      >  1:2002894/1‑68 (MQ=255)      tCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCa                                                   >  1:1841065/1‑67 (MQ=255)           ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:3054535/69‑1 (MQ=255)           ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:1974259/69‑1 (MQ=255)           ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:2043056/69‑1 (MQ=255)           ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:2592991/69‑1 (MQ=255)              gagaTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACat                                         <  1:1146226/69‑1 (MQ=255)              gagaTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACat                                         <  1:401171/69‑1 (MQ=255)                gaTCGTCCCCTACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1409210/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:691020/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1702160/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:406804/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1659416/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1797347/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1478213/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:191353/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1409519/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2982314/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:34510/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:3404758/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2598814/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1284605/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1067277/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2100533/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2147927/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1142637/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:3271970/60‑4 (MQ=255)                gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2481475/60‑4 (MQ=255)                 aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:1956605/1‑69 (MQ=255)                 aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:3184165/1‑69 (MQ=255)                 aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:959562/1‑69 (MQ=255)                 aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:952083/1‑69 (MQ=255)                 aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:1370992/1‑69 (MQ=255)                 aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACACATc                                      >  1:1379013/1‑69 (MQ=255)                 aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTATCCGCAGCCAAAAAGTACATATc                                      >  1:3316488/1‑69 (MQ=255)                 aTCGACCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:1274250/1‑69 (MQ=255)                   cGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCg                                     >  1:1350946/1‑68 (MQ=255)                         cAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACg                              >  1:512223/1‑69 (MQ=255)                         cAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACg                              >  1:2486762/1‑69 (MQ=255)                               gTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:2484042/69‑1 (MQ=255)                                tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:1460458/68‑1 (MQ=255)                                tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:1883399/68‑1 (MQ=255)                                tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:631514/68‑1 (MQ=255)                                tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:2056456/68‑1 (MQ=255)                                tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:1061653/68‑1 (MQ=255)                                tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:854722/68‑1 (MQ=255)                                  aGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgc                     >  1:3022510/1‑69 (MQ=255)                                  aGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgc                     >  1:1044282/1‑69 (MQ=255)                                    aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTgcgcgc                   >  1:3254618/1‑69 (MQ=255)                                    aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgcgc                   >  1:2166222/1‑69 (MQ=255)                                    aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgcgc                   >  1:800342/1‑69 (MQ=255)                                    aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgcgc                   >  1:1122722/1‑69 (MQ=255)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3236482/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3272100/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:440926/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:734785/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3301241/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3229626/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:711580/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3222384/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:707668/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:683000/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:571602/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1037626/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1989534/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1344377/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1210763/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1523768/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1537530/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1081351/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1692197/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1782095/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:186064/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1997158/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2102633/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2360175/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2375132/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2473070/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1015/3‑50 (MQ=25)                                        caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2992143/3‑50 (MQ=25)                                        aaccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1075678/3‑50 (MQ=25)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCTCTGAACGGTATTGAg                      >  1:3386127/1‑55 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:96022/1‑66 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:938017/1‑66 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:1040414/1‑66 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:1633524/1‑66 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:499362/1‑66 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:1989145/1‑66 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:2684652/1‑66 (MQ=255)                                               cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCACGTGCgg           >  1:647950/1‑66 (MQ=255)                                                     gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:1371458/69‑1 (MQ=255)                                                     gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:1443310/69‑1 (MQ=255)                                                     gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:1227816/69‑1 (MQ=255)                                                     gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:2213096/69‑1 (MQ=255)                                                     gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:3071822/69‑1 (MQ=255)                                                                    |                                                      CCACGTCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAGCGCG  >  minE/2400848‑2400968

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