Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	61,929	G→A	F510F (TTC→TTT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	61,929	0	G	A	90.7%	122.1 / 4.6	43	F510F (TTC→TTT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (17/22);  total (19/24)
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.27e-01

CGCTTTGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGC  >  minE/61863‑61986                                                                   |                                                          cGCTTTGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGc                                                       <  1:182680/71‑1 (MQ=255) cGCTTTGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGAGCGAAAAAGc                                                       <  1:1782563/71‑1 (MQ=255)  gCTTTGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCg                                                      >  1:914863/1‑71 (MQ=255)   cTTTGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCgg                                                     >  1:1411395/1‑71 (MQ=255)   cTTTGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCgg                                                     <  1:2393528/71‑1 (MQ=255)    tttGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCgg                                                     >  1:1208858/1‑70 (MQ=255)      tGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCa                                                   <  1:181501/70‑1 (MQ=255)      tGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAg                                                  >  1:946603/1‑71 (MQ=255)      tGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAg                                                  >  1:903361/1‑71 (MQ=255)      tGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAg                                                  >  1:1375911/1‑71 (MQ=255)        gTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAgg                                                 <  1:1386748/70‑1 (MQ=255)        gTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGc                                                >  1:1697201/1‑71 (MQ=255)         tttGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCg                                               <  1:1877267/71‑1 (MQ=255)               gCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCggtggtt                                         >  1:834859/1‑70 (MQ=255)                cATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGc                                        >  1:1020083/1‑71 (MQ=255)                cATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGc                                        >  1:961323/1‑71 (MQ=255)                 aTGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGcc                                       <  1:2063428/71‑1 (MQ=255)                 aTGATCTGATGACCACGCATGGTGATAGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGcc                                       <  1:68231/71‑1 (MQ=255)                  tGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCgaaaaa                                                         <  1:1109204/52‑1 (MQ=255)                    aTCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCgg                                                     <  1:1925140/54‑1 (MQ=255)                     tCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCgg                                                     <  1:1911078/53‑1 (MQ=255)                     tCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCggtgg                                           <  1:2109757/63‑1 (MQ=255)                       tgatgaCCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGgtggtg                                         >  1:2141632/1‑63 (MQ=255)                         atgaCCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACg                                <  1:1836008/70‑1 (MQ=255)                         atgaCCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGc                               <  1:719802/71‑1 (MQ=255)                             ccACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCg                                               >  1:2114427/1‑51 (MQ=255)                              cACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATaa                          <  1:1835772/71‑1 (MQ=255)                              cACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATaa                          <  1:750797/71‑1 (MQ=255)                                cGCATGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATaa                          >  1:1803621/1‑69 (MQ=255)                                   aTGGTGATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGc                     <  1:667012/71‑1 (MQ=255)                                        gATCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCggtg                                            >  1:377302/1‑43 (MQ=255)                                         aTCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat               <  1:1023813/71‑1 (MQ=255)                                          tCGACGATGCCAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGc                               <  1:1966775/54‑1 (MQ=255)                                            gacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatgat            >  1:749627/1‑71 (MQ=255)                                            gacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatga             >  1:1862627/1‑70 (MQ=255)                                                   ccAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt      <  1:2393370/70‑1 (MQ=255)                                                   ccAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt      <  1:2363956/70‑1 (MQ=255)                                                   ccAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTAATGATGAtttt      <  1:2211882/70‑1 (MQ=255)                                                    cAGCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGc                               <  1:1251797/44‑1 (MQ=255)                                                      gCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAgc  >  1:1159282/1‑71 (MQ=255)                                                      gCCGCGGATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg   >  1:1613618/1‑70 (MQ=255)                                                           ggATCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATa                           <  1:2342744/41‑1 (MQ=255)                                                              tCGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCa    >  1:770677/1‑61 (MQ=255)                                                               cGAAAAAGCGGCAGGCGGTGGTGCCGAGCACGCCATaaa                         <  1:1113091/39‑1 (MQ=255)                                                                gaagaaGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCAtt                  <  1:919380/45‑1 (MQ=255)                                                                gaaaaaGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCAt                   <  1:451785/44‑1 (MQ=255)                                                                   |                                                          CGCTTTGGTTTGCCGCATGATCTGATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGC  >  minE/61863‑61986

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