Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	734,426	A→G	8.1%	K16K (AAA→AAG)	tmk →	thymidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	734,426	0	A	G	8.1%	170.8 / 2.6	84	K16K (AAA→AAG)	tmk	thymidylate kinase
Reads supporting (aligned to +/- strand):  ref base A (34/43);  new base G (5/2);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGTGCTTAAGGAAAAAAATGCGCAGTAAGTATATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCAC  >  minE/734362‑734491                                                                 |                                                                  aGTGCTTAAGGAAAAAAATGCGCAGTAAGTATATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAAc                                                                 <  1:345382/67‑1 (MQ=255) aGTGCTTAAGGAAAAAAATGCGCAGTAAGTATATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAAc                                                                 <  1:1984089/67‑1 (MQ=255) aGTGCTTAAGGAAAAAAATGCGCAGTAAGTATATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAAc                                                                 <  1:2334191/67‑1 (MQ=255)           gAAAAAAATGCGCAGTAAGTATATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTa                                                      <  1:2745066/68‑1 (MQ=255)           gAAAAAAATGCGCAGTAAGTATATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTa                                                      <  1:2760077/68‑1 (MQ=255)                            agtaTATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGc                                    >  1:3254342/1‑69 (MQ=255)                            agtaTATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGc                                    >  1:839756/1‑69 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGc                                    >  1:3348561/1‑67 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGc                                    >  1:2735405/1‑67 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGc                                    >  1:2533096/1‑67 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGc                                    >  1:2791930/1‑67 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGc                                    >  1:885220/1‑67 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:334260/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:3271403/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:2390225/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:2179417/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1971648/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1953114/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:177870/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:2668157/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:3030614/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:3408281/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1698840/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1697224/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1536154/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:152758/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1474351/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1289713/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:125999/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1187566/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1147805/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:1143725/69‑2 (MQ=255)                              tataTCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:2903561/69‑2 (MQ=255)                               atatCGTCATTGGGGGGCTGGAAGGCGCAGGCAAACCTACCGCGCGTAATGTGGTGGTTGAGACGCTCg                                 >  1:2117978/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAGACTACCGCGCGTAATGTGGTGGTTGAGACGCTCg                                 >  1:2000900/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAGACTACCGCGCGTAATGTGGTGGTTGAGACGCTCg                                 >  1:1110047/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAGACTACCGCGCGTAATGTGGTGGTTGAGACGCTCg                                 >  1:1773273/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAGACTACCGCGCGTAATGTGGTGGTTGAGACGCTCg                                 >  1:2223715/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAGACTACCGCGCGTAATGAGGTGGTTGAGACGCTCg                                 >  1:238225/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCg                                 >  1:2122250/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCg                                 >  1:865047/1‑69 (MQ=255)                               atatCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:243621/68‑2 (MQ=255)                                  tCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACcg                                                            >  1:2956420/1‑39 (MQ=255)                                  tCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGc                              >  1:2581064/1‑69 (MQ=255)                                  tCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGc                              >  1:2823492/1‑69 (MQ=255)                                  tCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGc                              >  1:2290977/1‑69 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCaa                            >  1:3370313/1‑69 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:988106/68‑1 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:564830/68‑1 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:626151/68‑1 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:3077465/68‑1 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:2978805/68‑1 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:2782734/68‑1 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:2113940/68‑1 (MQ=255)                                    gTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:1505334/68‑1 (MQ=255)                                         tGAGGTGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAACGTGGTGGTTGAGACGCtt                                  <  1:39863/58‑2 (MQ=255)                                            ggggCTGGAAGGCGCAGGCAAGACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTa                     <  1:1703028/68‑1 (MQ=255)                                            ggggCTGGAAGGCGCAGGCAAGACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTa                     <  1:110085/68‑1 (MQ=255)                                               gCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATcc                  >  1:1660386/1‑68 (MQ=255)                                                cTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCcgc                <  1:742608/69‑1 (MQ=255)                                                cTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCcgc                <  1:3385461/69‑1 (MQ=255)                                                cTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCcgc                <  1:738439/69‑1 (MQ=255)                                                cTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCcgc                <  1:2245976/69‑1 (MQ=255)                                                 tGGTAGGCGCAGGCAAAACTAC‑‑CGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             <  1:396894/53‑1 (MQ=255)                                                       gcgcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATgg         >  1:3312647/1‑69 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:2890473/1‑47 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:1301242/1‑47 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:1209454/1‑47 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:1123355/1‑47 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:956023/1‑47 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:1000807/1‑47 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:2110073/1‑47 (MQ=255)                                                         gcAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCa                             >  1:2233809/1‑47 (MQ=255)                                                           aGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACTTGGgttt     <  1:2644822/69‑1 (MQ=255)                                                           aGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGtttt     <  1:33011/69‑1 (MQ=255)                                                            ggCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTggg                       >  1:2984586/1‑50 (MQ=255)                                                            ggCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTc    >  1:1734813/1‑69 (MQ=255)                                                            ggCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTc    >  1:3354245/1‑69 (MQ=255)                                                            ggCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTc    >  1:1749016/1‑69 (MQ=255)                                                            ggCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTc    >  1:2112896/1‑69 (MQ=255)                                                              cAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCAc  >  1:2086378/1‑69 (MQ=255)                                                              cAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCAc  >  1:3270908/1‑69 (MQ=255)                                                              cAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCAc  >  1:899404/1‑69 (MQ=255)                                                              cAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCAc  >  1:2820702/1‑69 (MQ=255)                                                                 |                                                                  AGTGCTTAAGGAAAAAAATGCGCAGTAAGTATATCGTCATTGAGGGGCTGGAAGGCGCAGGCAAAACTACCGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCAC  >  minE/734362‑734491

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