Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,803,427	C→T	54.1%	S80S (AGC→AGT)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,803,427	0	C	T	54.1%	9.1 / 84.2	99	S80S (AGC→AGT)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base T (17/36);  total (35/64)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GACGATTTTTACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTG  >  minE/2803359‑2803491                                                                     |                                                                 gACGATTTTTACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGc                                                                  >  1:941319/1‑69 (MQ=255)        tttACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCg                                                           >  1:1145690/1‑69 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:3311690/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2088361/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2081698/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2258977/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:1986282/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2321793/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2325661/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2747154/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:296094/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:1339209/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:1147586/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:754854/69‑1 (MQ=255)               ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:575617/69‑1 (MQ=255)                    cacaCCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTActc                                               >  1:551182/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCtt                                             >  1:1244934/1‑68 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:2706027/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1020280/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:2078331/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1184790/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:614211/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1814971/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:248691/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:2627235/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1300710/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:3052250/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:31330/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1357075/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:401958/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1338827/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:3396803/1‑69 (MQ=255)                       acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:3428960/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:2027805/69‑1 (MQ=255)                           tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:3112402/69‑1 (MQ=255)                           tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:2466577/69‑1 (MQ=255)                           tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:885139/69‑1 (MQ=255)                           tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:2684307/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:3348489/1‑63 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:186613/1‑63 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:3142452/1‑63 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:356790/1‑63 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:657031/1‑63 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCg                                    >  1:1267321/1‑68 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:2263738/1‑69 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:2673569/1‑69 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:1169247/1‑69 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:245036/1‑69 (MQ=255)                                tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:2451892/1‑69 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTTATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:3269905/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2553773/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:337460/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:132900/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:3421255/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:315000/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2363066/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1375952/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1423375/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2934170/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:650432/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2864139/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2857353/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:279841/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2771554/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1533163/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1566619/68‑1 (MQ=255)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1264739/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)                                           aTGGCGCGTTTGCAGGAAAGCGGTAGTCCAATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1426088/68‑1 (MQ=255)                                            tGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1280100/67‑1 (MQ=255)                                            tGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTgg                                <  1:2631724/60‑1 (MQ=255)                                                    ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:76406/69‑1 (MQ=255)                                                    ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1012141/69‑1 (MQ=255)                                                    ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:3190555/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:1372530/68‑1 (MQ=255)                                                     tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:2211413/68‑1 (MQ=255)                                                     tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:2684050/68‑1 (MQ=255)                                                     tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:3052348/68‑1 (MQ=255)                                                     tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1324618/68‑1 (MQ=255)                                                        aGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt            <  1:2779217/68‑1 (MQ=255)                                                        aGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt            <  1:2600957/68‑1 (MQ=255)                                                         ggAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt            <  1:2605038/67‑1 (MQ=255)                                                                ggTAGGCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1972487/57‑1 (MQ=255)                                                                 gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:37620/69‑1 (MQ=255)                                                                 gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:2689008/69‑1 (MQ=255)                                                                 gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:1609581/69‑1 (MQ=255)                                                                 gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:2524284/69‑1 (MQ=255)                                                                 gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:3357789/69‑1 (MQ=255)                                                                     |                                                                 GACGATTTTTACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTG  >  minE/2803359‑2803491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑