Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,303,538	T→C	60.9%	F84L (TTT→CTT)	sbcB →	exonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,303,538	0	T	C	60.9%	35.7 / 49.0	46	F84L (TTT→CTT)	sbcB	exonuclease I
Reads supporting (aligned to +/- strand):  ref base T (12/6);  new base C (14/14);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCG  >  minE/1303483‑1303605                                                        |                                                                    aattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTAtt                                                       <  1:1414972/70‑1 (MQ=255) aattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTgccg                                                               >  1:4151483/1‑62 (MQ=255)    taCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACt                                                   >  1:980360/1‑71 (MQ=255)    taCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACt                                                   >  1:19058/1‑71 (MQ=255)          tATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCttt                                             >  1:3296131/1‑71 (MQ=255)           aTTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCttt                                             >  1:1730645/1‑70 (MQ=255)            ttACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCtttt                                            <  1:3171513/70‑1 (MQ=255)            ttACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCtt                                              >  1:1777728/1‑68 (MQ=255)              aCCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtac                                         >  1:329966/1‑71 (MQ=255)              aCCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTta                                          >  1:1914073/1‑70 (MQ=255)               ccccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtacc                                        <  1:1160859/71‑1 (MQ=255)               ccccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtacc                                        <  1:3368709/71‑1 (MQ=255)                ccccCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtacc                                        <  1:3992270/70‑1 (MQ=255)                cccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtaccg                                       <  1:525072/71‑1 (MQ=255)                cccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtaccg                                       <  1:2924260/71‑1 (MQ=255)                        aaGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAg                               <  1:4186490/71‑1 (MQ=255)                        aaGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAg                               <  1:4042553/71‑1 (MQ=255)                         aGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAgc                              >  1:4115928/1‑70 (MQ=255)                         aGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGa                              >  1:3775071/1‑71 (MQ=255)                         aGCACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGa                              >  1:2630547/1‑71 (MQ=255)                           cACGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtaccgtaccg                                  <  1:3085844/65‑1 (MQ=255)                            aCGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCt                           >  1:2847244/1‑71 (MQ=255)                            aCGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCt                           >  1:2575934/1‑71 (MQ=255)                             cGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtaccg                                       <  1:1004509/58‑1 (MQ=255)                             cGGGCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtaccg                                       <  1:2367759/58‑1 (MQ=255)                              gggCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAg                               <  1:1065354/65‑1 (MQ=255)                              gggCGAAAGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTtaccg                                       <  1:4181241/57‑1 (MQ=255)                                     aGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtaccgtaccg                                  >  1:1256187/1‑55 (MQ=255)                                     aGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTc                     >  1:2362448/1‑68 (MQ=255)                                     aGGAGAAAACGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTc                     >  1:4213329/1‑68 (MQ=255)                                      ggAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtac                                         <  1:1340836/47‑1 (MQ=255)                                       gagaAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGc                >  1:3703283/1‑71 (MQ=255)                                          aaaaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTa              >  1:1027841/1‑70 (MQ=255)                                          aaaaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTa              >  1:3942129/1‑70 (MQ=255)                                          aaaaCGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCt               <  1:1141074/69‑1 (MQ=255)                                            aaCGAAGCCGCGCTTGCCGCCCGTATTCACTCGCtttt                                            >  1:4157838/1‑38 (MQ=255)                                            aaCGAAGCCGCGCTTGCCGCCCGTATTCACTCGCtttt                                            >  1:4251513/1‑38 (MQ=255)                                            aaCGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTacaa           >  1:4207284/1‑71 (MQ=255)                                             aCGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTacaac          <  1:3852533/71‑1 (MQ=255)                                             aCGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTacaac          <  1:3409413/71‑1 (MQ=255)                                             aCGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTacaac          <  1:2289509/71‑1 (MQ=255)                                              cGAAGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTacaaca         <  1:380405/71‑1 (MQ=255)                                                 aGCCGCGCTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAAtg      >  1:1337108/1‑71 (MQ=255)                                                     gcgTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCg  >  1:3604116/1‑71 (MQ=255)                                                     gcgTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCg  >  1:484250/1‑71 (MQ=255)                                                     gcgTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCg  >  1:1874417/1‑71 (MQ=255)                                                        |                                                                    AATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCG  >  minE/1303483‑1303605

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