Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,526,945	C→T	62.7%	E13K (GAG→AAG)	ttk ←	division inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,526,945	0	C	T	62.7%	46.4 / 61.6	59	E13K (GAG→AAG)	ttk	division inhibitor
Reads supporting (aligned to +/- strand):  ref base C (15/7);  new base T (13/24);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GCCGTCGTGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATAAGCACTT  >  minE/2526878‑2527011                                                                    |                                                                   gCCGTCGTGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTtgc                                                                  <  1:3195920/70‑1 (MQ=255) gCCGTCGTGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTcgc                                                                  <  1:887072/70‑1 (MQ=255) gCCGTCGTGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTcgc                                                                  <  1:1160287/70‑1 (MQ=255)    gtcgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACg                                                              >  1:957764/1‑71 (MQ=255)    gtcgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACg                                                              >  1:2654009/1‑71 (MQ=255)     tcgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACgg                                                             >  1:1092342/1‑71 (MQ=255)     tcgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACgg                                                             >  1:3430017/1‑71 (MQ=255)     tcgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACgg                                                             >  1:2799708/1‑71 (MQ=255)     tcgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACgg                                                             >  1:3357788/1‑71 (MQ=255)      cgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGt                                                            >  1:3476055/1‑71 (MQ=255)      cgtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGt                                                            >  1:2960472/1‑71 (MQ=255)       gtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGtt                                                           >  1:824244/1‑71 (MQ=255)       gtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGt                                                            <  1:3439298/70‑1 (MQ=255)       gtGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGt                                                            <  1:941603/70‑1 (MQ=255)          aTACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTcgc                                                                  <  1:2004091/61‑1 (MQ=255)          aTACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTcgc                                                                  <  1:163281/61‑1 (MQ=255)          aTACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCt                                                        >  1:2211026/1‑71 (MQ=255)             cGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCttt                                                      <  1:1557918/70‑1 (MQ=255)             cGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCttt                                                      <  1:1781739/70‑1 (MQ=255)               ttGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCg                                                   <  1:1412364/71‑1 (MQ=255)                   cTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAgt                                               <  1:117116/71‑1 (MQ=255)                   cTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAgt                                               >  1:1916627/1‑71 (MQ=255)                   cTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAgt                                               <  1:2074228/71‑1 (MQ=255)                       cATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAgtttgt                                           <  1:2614003/71‑1 (MQ=255)                       cATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAgtttg                                            <  1:2649733/70‑1 (MQ=255)                        aTCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCa                                                 >  1:2078079/1‑64 (MQ=255)                         tCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAgtttgttt                                         <  1:904327/71‑1 (MQ=255)                         tCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAgtttgttt                                         <  1:1188105/71‑1 (MQ=255)                         tCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAgtttgttt                                         <  1:1413189/71‑1 (MQ=255)                           gCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGtttt                                        <  1:582727/70‑1 (MQ=255)                            cTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTc                                      >  1:147066/1‑71 (MQ=255)                            cTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTc                                      >  1:573738/1‑71 (MQ=255)                            cTGGATTCCAGCATCAGAGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTc                                      <  1:2004282/71‑1 (MQ=255)                             tGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTtgc                                                                  <  1:1745243/42‑1 (MQ=255)                             tGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTc                                      <  1:1051137/70‑1 (MQ=255)                               gATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGc                                   <  1:2293376/71‑1 (MQ=255)                                 ttCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCATTTCGCAGTTTGTTTTTCTGCCa                                 >  1:1197096/1‑71 (MQ=255)                                   ccAGCATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTc                                                    >  1:3659298/1‑50 (MQ=255)                                   ccAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTcc                                                         >  1:2408431/1‑45 (MQ=255)                                     aGCATCAGCGCCAGAGCCTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGt                              >  1:1790786/1‑70 (MQ=255)                                       cATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCtt                                                       <  1:3447283/43‑1 (MQ=255)                                       cATCAGCGCCAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTa                            >  1:780626/1‑70 (MQ=255)                                               ccAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAgt                                               <  1:3087082/43‑1 (MQ=255)                                                cAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATAcg                   <  1:1686168/70‑2 (MQ=255)                                                cAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATAccc                  <  1:2607271/71‑1 (MQ=255)                                                cAGAGACTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATAccc                  <  1:1094081/71‑1 (MQ=255)                                                  gagaCTGAAGTATTTCCTTGCGACGGTTCCTTTTCGc                                                  >  1:497423/1‑37 (MQ=255)                                                  gagaCTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGGCATGTTACAAAATAccc                  <  1:3072212/69‑1 (MQ=255)                                                  gagaCTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCt                >  1:1682802/1‑71 (MQ=255)                                                  gagaCTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCt                >  1:2034655/1‑71 (MQ=255)                                                  gagaCTGAAGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGACATGTTACAAAATAccc                  >  1:145442/1‑69 (MQ=255)                                                   agaCTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAGGTTGttttt                                       <  1:2731384/47‑1 (MQ=255)                                                    gaCTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTg               >  1:317030/1‑70 (MQ=255)                                                      cTGAAGTATTTCCTCGCGGCGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGaaa            <  1:3591955/71‑1 (MQ=255)                                                         aaGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAAt          <  1:1554937/70‑1 (MQ=255)                                                         aaGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATa         <  1:781019/71‑1 (MQ=255)                                                         aaGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATa         <  1:125599/71‑1 (MQ=255)                                                         aaGTATTTCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAAt          >  1:907378/1‑70 (MQ=255)                                                          aGTATTTCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATaa        <  1:2869149/71‑1 (MQ=255)                                                            tATTTCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACaa                         <  1:394307/52‑1 (MQ=255)                                                            tATTTCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACaa                         <  1:3722103/52‑1 (MQ=255)                                                              tttCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATAAGCAc    >  1:1153279/1‑71 (MQ=255)                                                                tCCTTGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATAAGCACt   <  1:275252/70‑1 (MQ=255)                                                                tCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATAAGCACtt  >  1:3680184/1‑71 (MQ=255)                                                                tCCTCGCGACGGTTCCTTTTCGCAGTTTGTATTTCTGCCATGTTACAAAATACACCTGAAAATAAGCACt   >  1:661454/1‑70 (MQ=255)                                                                    |                                                                   GCCGTCGTGATACGTTGGCTTCCATCGCTGGATTCCAGCATCAGCGCCAGAGACTGAAGTATTTCCTCGCGACGGTTCCTTTTCGCAGTTTGTTTTTCTGCCATGTTACAAAATACCCCTGAAAATAAGCACTT  >  minE/2526878‑2527011

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