Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,947,335	A→G	21.1%	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	21.1%	198.9 / 50.5	96	L158P (CTG→CCG)	dnaT	DNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (29/46);  new base G (8/12);  total (37/59)
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.99e-01

ACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGT  >  minE/2947270‑2947401                                                                  |                                                                   aCCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTccg                                                               >  1:1324783/1‑71 (MQ=255)      ctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGtt                                                          >  1:1361176/1‑71 (MQ=255)      ctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGt                                                                  >  1:2225063/1‑63 (MQ=255)      ctGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTc                                                                 >  1:204228/1‑64 (MQ=255)         aaTCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGc                                                        >  1:25493/1‑70 (MQ=255)                gtggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACg                                                 <  1:4089650/70‑1 (MQ=255)                gtggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACg                                                 <  1:2911353/70‑1 (MQ=255)                gtggAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACg                                                 <  1:4085005/70‑1 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCg                                             <  1:3484343/71‑1 (MQ=255)                   gAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCg                                             <  1:2472059/71‑1 (MQ=255)                    aaTTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCTGTCCGCCGTTGCTGGCACGACCGa                                            <  1:352130/71‑1 (MQ=255)                    aaTTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGa                                            <  1:4003542/71‑1 (MQ=255)                    aaTTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGa                                            <  1:2090146/71‑1 (MQ=255)                    aaTTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGa                                            <  1:1221745/71‑1 (MQ=255)                    aaTTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGa                                            <  1:3325448/71‑1 (MQ=255)                    aaTTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGa                                            <  1:3895453/71‑1 (MQ=255)                     aTTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            >  1:3020577/1‑70 (MQ=255)                        tGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTg                                        >  1:1607478/1‑71 (MQ=255)                        tGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTg                                        >  1:3725634/1‑71 (MQ=255)                              tCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:2791620/59‑1 (MQ=255)                              tCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:35618/59‑1 (MQ=255)                              tCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:2256399/59‑1 (MQ=255)                              tCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:1574332/59‑1 (MQ=255)                              tCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:708792/59‑1 (MQ=255)                              tCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:1684319/59‑1 (MQ=255)                                atgTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgc                                <  1:2221708/69‑1 (MQ=255)                                aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgc                                >  1:1769864/1‑71 (MQ=255)                                aGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgc                                >  1:1950453/1‑71 (MQ=255)                                 ggTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGa                                            <  1:2094236/58‑1 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGc                                                        <  1:2450629/45‑1 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTc                                                                 >  1:1501724/1‑36 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTc                                                                 >  1:3803706/1‑36 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTc                                                                 >  1:762214/1‑36 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:1124813/55‑1 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:1377951/55‑1 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:3659447/55‑1 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:1916802/55‑1 (MQ=255)                                  gTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:2832298/55‑1 (MQ=255)                                   ttCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGAcc                                              <  1:3332479/54‑1 (MQ=255)                                    tCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgcgcgc                            >  1:1524957/1‑71 (MQ=255)                                    tCGCTGACCGTATTCACATCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTgcgcgc                            >  1:311235/1‑71 (MQ=255)                                        tGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGt                        >  1:21792/1‑71 (MQ=255)                                           ccGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGtttt                     <  1:3467423/71‑1 (MQ=255)                                           ccGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGtttt                     <  1:1091479/71‑1 (MQ=255)                                           ccGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGtttt                     <  1:1760985/71‑1 (MQ=255)                                           cagTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGtttt                     <  1:1994468/69‑1 (MQ=255)                                             gTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttg                   <  1:818178/71‑1 (MQ=255)                                             gTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttg                   <  1:1813791/71‑1 (MQ=255)                                              tATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:3838195/45‑1 (MQ=255)                                              tATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:3876111/45‑1 (MQ=255)                                              tATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:1835968/45‑1 (MQ=255)                                              tATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:60049/45‑1 (MQ=255)                                              tATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:2747235/45‑1 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCTCCGTTGCTGGCACGACCGa                                            <  1:3415809/43‑1 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:3936601/43‑1 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:2914663/43‑1 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:2913711/43‑1 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:2413352/43‑1 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:275145/1‑45 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:3584105/1‑45 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:3567422/1‑45 (MQ=255)                                                ttCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttg                   <  1:4132374/68‑1 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGa                                            <  1:1368100/41‑1 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:2727565/1‑43 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:3587733/1‑43 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:3613810/1‑43 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:2137713/1‑43 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGAtt                                          >  1:2787214/1‑43 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCGAGTTTttgt                  >  1:443325/1‑67 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc               <  1:2575466/70‑1 (MQ=255)                                                  cacaTCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc               <  1:1161766/70‑1 (MQ=255)                                                      tctcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc               <  1:3350372/66‑1 (MQ=255)                                                      tctcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc               <  1:433365/66‑1 (MQ=255)                                                      tctcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGc               <  1:977629/66‑1 (MQ=255)                                                       ctcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTg          >  1:3933973/1‑70 (MQ=255)                                                        tcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCa             >  1:2883927/1‑66 (MQ=255)                                                        tcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAcgg          >  1:989778/1‑67 (MQ=255)                                                        tcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAcgg          >  1:303287/1‑67 (MQ=255)                                                        tcGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca        >  1:639289/1‑71 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcaca      >  1:2925978/1‑71 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcaca      >  1:493206/1‑71 (MQ=255)                                                          gTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcaca      >  1:192402/1‑71 (MQ=255)                                                          gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTTca        <  1:3905319/69‑1 (MQ=255)                                                          gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcaca      >  1:3840426/1‑71 (MQ=255)                                                          gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcaca      >  1:1062178/1‑71 (MQ=255)                                                          gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca        <  1:3808178/69‑1 (MQ=255)                                                          gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca        <  1:3747789/69‑1 (MQ=255)                                                          gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca        <  1:1515674/69‑1 (MQ=255)                                                          gTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGca        <  1:1722784/69‑1 (MQ=255)                                                          gTTTCGACGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGcaca      >  1:528806/1‑71 (MQ=255)                                                             tCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATgg   >  1:2012254/1‑71 (MQ=255)                                                               ggCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGcc                           <  1:1015762/45‑1 (MQ=255)                                                               ggCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGcc                           <  1:2291671/45‑1 (MQ=255)                                                               ggCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGt  <  1:296254/70‑1 (MQ=255)                                                               ggCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGt  <  1:3579810/70‑1 (MQ=255)                                                               ggCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGt  <  1:3860020/70‑1 (MQ=255)                                                                  |                                                                   ACCCTCTGAATCCTGGTGGAATTTGGCTGTCAGGTTCGCTGACCGTATTCACATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCACTGCACATGGT  >  minE/2947270‑2947401

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