Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	62,471	T→A	5.2%	T330S (ACT→TCT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	62,471	0	T	A	5.2%	239.9 / 3.9	77	T330S (ACT→TCT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (48/25);  new base A (2/2);  total (50/27)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGA  >  minE/62405‑62539                                                                   |                                                                     aaaaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGt                                                                      >  1:3990966/1‑67 (MQ=255) aaaaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGt                                                                      >  1:4814498/1‑67 (MQ=255) aaaaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGt                                                                      >  1:1309574/1‑67 (MQ=255) aaaaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGt                                                                      >  1:1522583/1‑67 (MQ=255) aaaaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGt                                                                      >  1:5227100/1‑67 (MQ=255) aaaaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGt                                                                      >  1:6889826/1‑67 (MQ=255)   aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCa                                                                 >  1:5404202/1‑70 (MQ=255)   aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAg                                                                >  1:2615766/1‑71 (MQ=255)   aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAg                                                                >  1:6787259/1‑71 (MQ=255)   aaTGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAg                                                                >  1:2338859/1‑71 (MQ=255)     tGGCATGATTTCATTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:5620236/71‑1 (MQ=255)     tGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:2180690/71‑1 (MQ=255)     tGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:6758609/71‑1 (MQ=255)     tGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                                              <  1:694694/71‑1 (MQ=255)          tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         <  1:252925/71‑1 (MQ=255)          tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         <  1:4466893/71‑1 (MQ=255)          tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         <  1:848712/71‑1 (MQ=255)          tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         >  1:1031841/1‑71 (MQ=255)          tGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGt                                                         >  1:1933420/1‑71 (MQ=255)           gATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGtt                                                        >  1:5307856/1‑71 (MQ=255)                 aGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCt                                                   <  1:6394136/70‑1 (MQ=255)                     ttGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              >  1:2461105/1‑71 (MQ=255)                      tgtgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTAt                                                     >  1:1799094/1‑63 (MQ=255)                      tgtgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              <  1:4014366/70‑1 (MQ=255)                      tgtgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              >  1:808554/1‑70 (MQ=255)                        tgGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGa                                            <  1:2167280/70‑1 (MQ=255)                         ggAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGa                                            >  1:2489917/1‑69 (MQ=255)                          gAAGATGCGCGTGACCAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGaaca                                         >  1:241191/1‑69 (MQ=255)                             gATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGaaa                                          <  1:5216757/67‑1 (MQ=255)                             gATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggc                                      >  1:6362388/1‑71 (MQ=255)                             gATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggc                                      >  1:1253820/1‑71 (MQ=255)                               tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTAt                                                     >  1:5266919/1‑54 (MQ=255)                               tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTAt                                                     >  1:1082083/1‑54 (MQ=255)                               tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:1019963/1‑71 (MQ=255)                               tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:5958180/1‑71 (MQ=255)                               tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:6150460/1‑71 (MQ=255)                               tCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcggcgg                                    >  1:3210811/1‑71 (MQ=255)                                cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCCGCGCAGGTTTATCTTCGGCGAAACGGCGGt                                   >  1:4736683/1‑71 (MQ=255)                                cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGt                                   >  1:932236/1‑71 (MQ=255)                                cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGt                                   >  1:607650/1‑71 (MQ=255)                                cTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCCGt                                   >  1:6365210/1‑71 (MQ=255)                                   cGTGACCAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGAACATGCGGTcga                                <  1:1547846/71‑3 (MQ=255)                                       aCCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCg                                             >  1:3229676/1‑54 (MQ=255)                                       aCCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCg                                             >  1:4784393/1‑54 (MQ=255)                                       aCCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCg                                             >  1:4776020/1‑54 (MQ=255)                                       aCCAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGAACATGCGGTCGAtctc                            >  1:2422970/1‑71 (MQ=255)                                        ccAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGa                                            <  1:2914662/54‑1 (MQ=255)                                        ccAGCTCGCAATCCTTGAGGTTGTACCAGGCCAGCGCCGGCTTGTCCTCGGCGAACATGCGGTCGATCTCg                           <  1:2096503/71‑1 (MQ=255)                                          aGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGc                                              >  1:2484461/1‑50 (MQ=255)                                           gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATc                                                    <  1:191269/43‑1 (MQ=255)                                           gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTc                                  >  1:3597422/1‑61 (MQ=255)                                           gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTc                                  >  1:6141302/1‑61 (MQ=255)                                           gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTc                                  >  1:6314886/1‑61 (MQ=255)                                           gCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTc                         >  1:5993921/1‑70 (MQ=255)                                            cTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGTCGGTCAATTTCGTCCa                       <  1:2241917/71‑1 (MQ=255)                                            cTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCTGTCAATTTCGTCCa                       <  1:2330194/71‑1 (MQ=255)                                            cTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCa                       <  1:3561457/71‑1 (MQ=255)                                              cGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGaaa                                          >  1:4517680/1‑50 (MQ=255)                                              cGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCaa                                >  1:5567953/1‑60 (MQ=255)                                              cGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCaa                                >  1:6584173/1‑60 (MQ=255)                                               gCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCAt                      <  1:5604015/69‑1 (MQ=255)                                                cAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCg                   <  1:23250/71‑1 (MQ=255)                                                  aTCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcg                                        >  1:3756049/1‑48 (MQ=255)                                                  aTCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTcc                        >  1:601340/1‑64 (MQ=255)                                                  aTCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTcc                        >  1:5704933/1‑64 (MQ=255)                                                   tCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATc                >  1:5338803/1‑71 (MQ=255)                                                     tttCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTc                         <  1:6052000/60‑1 (MQ=255)                                                       tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAAcg                                        >  1:2513290/1‑43 (MQ=255)                                                       tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCAt                      >  1:4784659/1‑61 (MQ=255)                                                       tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCAt                      >  1:6651069/1‑61 (MQ=255)                                                       tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATcc               <  1:5823581/68‑1 (MQ=255)                                                       tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCa             <  1:2528479/70‑1 (MQ=255)                                                       tCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCa             <  1:2566889/70‑1 (MQ=255)                                                          ggTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACggg         <  1:198289/71‑1 (MQ=255)                                                              ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTa      <  1:2125439/70‑1 (MQ=255)                                                              ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTa      >  1:178466/1‑70 (MQ=255)                                                              ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTa      >  1:3036434/1‑70 (MQ=255)                                                              ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTAt     >  1:1854288/1‑71 (MQ=255)                                                              ataAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTAt     >  1:2905641/1‑71 (MQ=255)                                                               taAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATc    <  1:5432507/71‑1 (MQ=255)                                                                aaGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATc    >  1:2200664/1‑70 (MQ=255)                                                                aaGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATc    >  1:1007033/1‑70 (MQ=255)                                                                 aGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGa  >  1:6806464/1‑71 (MQ=255)                                                                   |                                                                     AAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGA  >  minE/62405‑62539

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