Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,947,335	A→G	17.8%	L158P (CTG→CCG)	dnaT ←	DNA biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,947,335	0	A	G	17.8%	97.8 / 17.1	45	L158P (CTG→CCG)	dnaT	DNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base G (5/3);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TCGTTACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTGA  >  minE/2947265‑2947403                                                                       |                                                                     tCGTTACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCa                                                                      <  1:1159655/71‑1 (MQ=255) tCGTTACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCa                                                                      <  1:1162222/71‑1 (MQ=255)   gTTACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGt                                                                    >  1:319894/1‑71 (MQ=255)     tACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTcc                                                                  >  1:833659/1‑71 (MQ=255)       cccTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgc                                                                <  1:711768/71‑1 (MQ=255)         ctctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTccgc                                                                <  1:396922/69‑1 (MQ=255)         ctctGAATCCTGGT‑GAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGt                                                             <  1:1210468/71‑1 (MQ=255)          tctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGt                                                             >  1:2149135/1‑71 (MQ=255)          tctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGACGt                                                             >  1:503893/1‑71 (MQ=255)           ctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGtt                                                            >  1:834102/1‑71 (MQ=255)            tGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGt                                                             <  1:82822/69‑1 (MQ=255)               aTCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTg                                                        >  1:1124544/1‑71 (MQ=255)                  cTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGc                                                      <  1:487878/70‑1 (MQ=255)                  cTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGc                                                      <  1:2156600/70‑1 (MQ=255)                      tggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgc                                                                <  1:628444/56‑1 (MQ=255)                           tttGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGc                                                      <  1:1792964/61‑1 (MQ=255)                           tttGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGc                                                      <  1:2088874/61‑1 (MQ=255)                           tttGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                            >  1:2137647/1‑71 (MQ=255)                            ttGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGa                                                  >  1:517948/1‑64 (MQ=255)                             tGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAttt                                           <  1:502602/70‑1 (MQ=255)                               gCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTg                                          >  1:1558993/1‑69 (MQ=255)                                     aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgc                                  >  1:1029701/1‑71 (MQ=255)                                     aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTg                                   >  1:1082721/1‑70 (MQ=255)                                                ccGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGtttt                       <  1:777347/71‑1 (MQ=255)                                                ccGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGtttt                       <  1:1438285/71‑1 (MQ=255)                                                   tATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgt                    >  1:2053155/1‑71 (MQ=255)                                                   tATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgt                    <  1:1878312/71‑1 (MQ=255)                                                    aTTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt                   >  1:550666/1‑71 (MQ=255)                                                     ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGc                                                      <  1:1594428/35‑1 (MQ=255)                                                     ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGcc                             >  1:404758/1‑60 (MQ=255)                                                     ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGttttt                      <  1:67254/67‑1 (MQ=255)                                                     ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGttttt                      <  1:1438683/67‑1 (MQ=255)                                                      tCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGttttt                      <  1:1366691/66‑1 (MQ=255)                                                       cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGTCCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc                 >  1:1985721/1‑70 (MQ=255)                                                       cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAttt                                           >  1:762168/1‑44 (MQ=255)                                                        acaTCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGcc                <  1:1786375/70‑1 (MQ=255)                                                        acaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgt                    >  1:368216/1‑66 (MQ=255)                                                             tcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca          >  1:134559/1‑71 (MQ=255)                                                              cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCa               >  1:1064797/1‑65 (MQ=255)                                                              cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc              >  1:532574/1‑66 (MQ=255)                                                              cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc              >  1:1355496/1‑66 (MQ=255)                                                               gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca          <  1:901316/69‑1 (MQ=255)                                                               gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca          <  1:566563/69‑1 (MQ=255)                                                               gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTT‑TTGCCACTGcaca        <  1:1018201/70‑1 (MQ=255)                                                                     gCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTg   >  1:1993668/1‑70 (MQ=255)                                                                      cAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCa               >  1:1851624/1‑57 (MQ=255)                                                                      cAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTGa  <  1:1016655/70‑1 (MQ=255)                                                                       |                                                                     TCGTTACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTGA  >  minE/2947265‑2947403

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