Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,803,445	T→C	45.6%	T86T (ACT→ACC)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,803,445	0	T	C	45.6%	10.4 / 76.1	90	T86T (ACT→ACC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (19/30);  new base C (17/24);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTGGTTTTGCTTATCTGGCAGA  >  minE/2803378‑2803510                                                                    |                                                                  cacaCCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTActc                                                                  >  1:551182/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCtt                                                                >  1:1244934/1‑68 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:1300710/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:1814971/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:2078331/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:248691/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:2627235/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:2706027/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:1357075/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:1338827/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:3052250/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:31330/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:3396803/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:3428960/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:401958/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:1184790/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:614211/1‑69 (MQ=255)    acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                                               >  1:1020280/1‑69 (MQ=255)     ccGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTg                                                               >  1:3034850/1‑68 (MQ=255)     ccGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTg                                                               >  1:2065100/1‑68 (MQ=255)       gTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                                            >  1:960513/1‑69 (MQ=255)        tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                                           <  1:3112402/69‑1 (MQ=255)        tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                                           <  1:2684307/69‑1 (MQ=255)        tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                                           <  1:885139/69‑1 (MQ=255)        tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                                           <  1:2466577/69‑1 (MQ=255)        tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                                           <  1:2027805/69‑1 (MQ=255)        tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTAGCGGa                                                           <  1:2254475/69‑1 (MQ=255)            atCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATc                                                        <  1:2884277/68‑1 (MQ=255)            atCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCg                                                       >  1:3251129/1‑69 (MQ=255)             tCTTTACTTAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                                      >  1:390423/1‑69 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                                            >  1:3142452/1‑63 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                                            >  1:3348489/1‑63 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                                            >  1:356790/1‑63 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                                            >  1:186613/1‑63 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                                            >  1:657031/1‑63 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCg                                                       >  1:1267321/1‑68 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                                      >  1:245036/1‑69 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                                      >  1:2451892/1‑69 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                                      >  1:2263738/1‑69 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                                      >  1:2673569/1‑69 (MQ=255)             tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                                      >  1:1169247/1‑69 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTTATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:3269905/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1264739/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2857353/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:279841/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2864139/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:132900/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2771554/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:315000/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:337460/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2934170/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2363066/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1566619/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1533163/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:650432/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1215531/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1251465/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2633081/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1423375/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2553773/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:3421255/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1375952/68‑1 (MQ=255)                        aTGGCGCGTTTGCAGGAAAGCGGTAGTCCAATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1426088/68‑1 (MQ=255)                         tGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:1280100/67‑1 (MQ=255)                         tGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTgg                                                   <  1:2631724/60‑1 (MQ=255)                                 ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:3190555/69‑1 (MQ=255)                                 ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:1012141/69‑1 (MQ=255)                                 ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:76406/69‑1 (MQ=255)                                 ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:1633636/69‑1 (MQ=255)                                 ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:148451/69‑1 (MQ=255)                                  tGCAGGAAAGCGGTAGCCCTATCGATGTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:2731301/68‑1 (MQ=255)                                  tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:1324618/68‑1 (MQ=255)                                  tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:2211413/68‑1 (MQ=255)                                  tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:1372530/68‑1 (MQ=255)                                  tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:3052348/68‑1 (MQ=255)                                  tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:2684050/68‑1 (MQ=255)                                     aGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt                               <  1:2600957/68‑1 (MQ=255)                                     aGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt                               <  1:2779217/68‑1 (MQ=255)                                      ggAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt                               <  1:2605038/67‑1 (MQ=255)                                             ggTAGGCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc                                  <  1:1972487/57‑1 (MQ=255)                                              gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg                     <  1:2689008/69‑1 (MQ=255)                                              gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg                     <  1:3357789/69‑1 (MQ=255)                                              gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg                     <  1:2524284/69‑1 (MQ=255)                                              gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg                     <  1:37620/69‑1 (MQ=255)                                              gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg                     <  1:1609581/69‑1 (MQ=255)                                                    cTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                                            <  1:2316012/40‑1 (MQ=255)                                                     tATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCa                                      >  1:1488169/1‑45 (MQ=255)                                                     tATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTGGTTTTg               >  1:2226948/1‑68 (MQ=255)                                                         gATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTGGTTTTGCtt            <  1:169927/67‑1 (MQ=255)                                                         gATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTGGTTTTGCTTa           <  1:2810452/68‑1 (MQ=255)                                                                  aCTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTGGTTTTGCTTATCTGGCaga  >  1:2871210/1‑68 (MQ=255)                                                                    |                                                                  CACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTGGTTTTGCTTATCTGGCAGA  >  minE/2803378‑2803510

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