Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,947,335	A→G	53.7%	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	53.7%	18.3 / 92.4	67	L158P (CTG→CCG)	dnaT	DNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (14/17);  new base G (7/29);  total (21/46)
Fisher's exact test for biased strand distribution p-value = 3.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

ACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTGA  >  minE/2947270‑2947403                                                                  |                                                                     aCCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccg                                                                 >  1:1653739/1‑71 (MQ=255)  cccTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgc                                                                <  1:39198/71‑1 (MQ=255)   ccTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgcc                                                               <  1:1980798/71‑1 (MQ=255)   ccTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgcc                                                               <  1:2714324/71‑1 (MQ=255)   ccTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgc                                                                <  1:1596367/70‑1 (MQ=255)    ctctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgccg                                                              >  1:3204615/1‑71 (MQ=255)     tctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGt                                                             >  1:2419061/1‑71 (MQ=255)               ggtggAATTTTGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccgcc                                                               <  1:1654094/59‑1 (MQ=255)                 tggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGt                                                             <  1:1475704/59‑1 (MQ=255)                 tggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGt                                                             <  1:2222525/59‑1 (MQ=255)                 tggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGt                                                             <  1:1061829/59‑1 (MQ=255)                 tggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGt                                                             <  1:1021886/59‑1 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGtt                                                            >  1:1159475/1‑58 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGtt                                                            >  1:3344200/1‑58 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGtt                                                            >  1:9276/1‑58 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCg                                               <  1:2593761/71‑1 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCg                                               <  1:2482518/71‑1 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCg                                               <  1:692491/71‑1 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCg                                               <  1:1545995/71‑1 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCg                                               <  1:909779/71‑1 (MQ=255)                          gCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCa                                        <  1:2413617/71‑1 (MQ=255)                           cTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTccg                                                                 <  1:3413191/45‑1 (MQ=255)                           cTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTccg                                                                 <  1:3379662/45‑1 (MQ=255)                           cTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTccg                                                                 <  1:3681731/45‑1 (MQ=255)                           cTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTccg                                                                 <  1:2837370/45‑1 (MQ=255)                           cTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTccg                                                                 <  1:3876710/45‑1 (MQ=255)                           cTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTccg                                                                 <  1:2031517/45‑1 (MQ=255)                           cTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAg                                       <  1:1248662/71‑1 (MQ=255)                                aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTg                                   <  1:835481/70‑1 (MQ=255)                                aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTg                                   <  1:1323737/70‑1 (MQ=255)                                aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTg                                   <  1:635008/70‑1 (MQ=255)                                aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTg                                   <  1:1936727/70‑1 (MQ=255)                                   ttCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgcgc                                >  1:1157564/1‑70 (MQ=255)                                    tCGCTGCCCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgcgcg                               >  1:634741/1‑70 (MQ=255)                                     cGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgcgcgc                              >  1:3478658/1‑70 (MQ=255)                                     cGCTGACCGTATTCACATCTCGTTTCGGCAGTccgcc                                                               <  1:2254773/37‑1 (MQ=255)                                         gACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGTCACGACCGATTTGCAGGCTGCGCGCCAGtt                         >  1:850372/1‑71 (MQ=255)                                         gACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGtt                         >  1:1500685/1‑71 (MQ=255)                                             gTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttg                     <  1:1451751/71‑1 (MQ=255)                                             gTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttg                     <  1:3317947/71‑1 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt                   <  1:880879/70‑1 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc                 <  1:1665621/70‑1 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc                 <  1:3634072/70‑1 (MQ=255)                                                    caTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAgg                                      >  1:622040/1‑47 (MQ=255)                                                    caTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGttt                        >  1:977122/1‑61 (MQ=255)                                                    caTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGttt                        >  1:1182496/1‑61 (MQ=255)                                                     aTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGttt                        <  1:4260380/60‑1 (MQ=255)                                                     aTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGttt                        <  1:1640253/60‑1 (MQ=255)                                                       ctcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg            >  1:2689839/1‑70 (MQ=255)                                                       ctcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg            >  1:2629945/1‑70 (MQ=255)                                                       ctcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg            >  1:643685/1‑70 (MQ=255)                                                        tcGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca          >  1:2436366/1‑71 (MQ=255)                                                         cGTTTCGGCAGTCCGCCGTTGCTGGCAGGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca          >  1:2084265/1‑70 (MQ=255)                                                         cGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca          >  1:1508093/1‑70 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg            <  1:1110531/67‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg            <  1:3971871/67‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg            <  1:3979365/67‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg            <  1:794692/67‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcaca        <  1:2386090/71‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcac         <  1:1023830/70‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcac         <  1:2038525/70‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcac         <  1:3771930/70‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca          <  1:1417486/69‑1 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCACGCTGCGCGCCAGTTTTTGTTGCCACTGcac         <  1:2706163/70‑1 (MQ=255)                                                             tCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATg      >  1:1311230/1‑70 (MQ=255)                                                              cggcggTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCa               >  1:1732475/1‑60 (MQ=255)                                                              cggcggTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCa               >  1:839413/1‑60 (MQ=255)                                                                gCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTGa  <  1:35136/71‑1 (MQ=255)                                                                gCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTGa  <  1:1056137/71‑1 (MQ=255)                                                                  |                                                                     ACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGTGA  >  minE/2947270‑2947403

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