Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,232,360	C→G	9.7%	M132I (ATG→ATC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,232,360	0	C	G	9.7%	210.7 / 11.2	72	M132I (ATG→ATC)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base G (3/4);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

ATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCA  >  minE/1232292‑1232429                                                                     |                                                                      aTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAt                                                                     <  1:21634/71‑1 (MQ=255) aTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAt                                                                     <  1:3233465/71‑1 (MQ=255)    ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                                  <  1:1742057/71‑1 (MQ=255)    ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                                  <  1:3652393/71‑1 (MQ=255)      gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                >  1:242790/1‑71 (MQ=255)      gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAAAGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                >  1:2515542/1‑71 (MQ=255)              gTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAtt                                                                    >  1:1028786/1‑59 (MQ=255)               tCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAtt                                                                    >  1:3405419/1‑58 (MQ=255)                 gggAAAGAAAATCTTCTGGCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                                     <  1:2173542/71‑1 (MQ=255)                 gggAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                                     <  1:1189974/71‑1 (MQ=255)                   gaaagaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaaaa                                                   <  1:1435691/71‑1 (MQ=255)                      agaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                >  1:1766811/1‑55 (MQ=255)                       gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                               <  1:1310724/71‑1 (MQ=255)                       gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                               <  1:2323041/71‑1 (MQ=255)                       gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                               >  1:546977/1‑71 (MQ=255)                       gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                               >  1:1637556/1‑71 (MQ=255)                        aaaaTCTTCTGCCAGCGACATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                              >  1:2319817/1‑71 (MQ=255)                         aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                              <  1:1746360/70‑1 (MQ=255)                         aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                              >  1:2356153/1‑70 (MQ=255)                         aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                             <  1:1673387/71‑1 (MQ=255)                         aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                             <  1:14212/71‑1 (MQ=255)                           aTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAtt                                                                    >  1:2485635/1‑46 (MQ=255)                            tcttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCa                                                                      >  1:765386/1‑43 (MQ=255)                            tcttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCa                                                                      >  1:2741840/1‑43 (MQ=255)                            tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                              >  1:3135235/1‑51 (MQ=21)                            tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                              >  1:1842334/1‑51 (MQ=21)                            tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                              >  1:2025436/1‑51 (MQ=21)                             cttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTAt                                           >  1:107123/1‑69 (MQ=255)                               tctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                       >  1:1840121/1‑71 (MQ=255)                                ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                          <  1:2257400/67‑1 (MQ=255)                                ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                          <  1:43657/67‑1 (MQ=255)                                ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                          <  1:1772999/67‑1 (MQ=255)                                    cAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgc                                  <  1:1553584/71‑1 (MQ=255)                                      gcgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                               <  1:3221466/56‑1 (MQ=255)                                       cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                               >  1:329111/1‑71 (MQ=255)                                       cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                               >  1:355919/1‑71 (MQ=255)                                        gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                              <  1:2452894/71‑1 (MQ=255)                                        gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                              <  1:3243234/71‑1 (MQ=255)                                         ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgcgc                                <  1:2305181/68‑1 (MQ=255)                                         ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                             >  1:2349938/1‑71 (MQ=255)                                          cATAATGGTGTTATCGGGCGTAATCGCCATTCCGCACATCGAATAAAACTCACTATGAGCTGCGCGCTTAt                            <  1:3102388/71‑1 (MQ=255)                                           aTAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTAtt                           >  1:664892/1‑71 (MQ=255)                                            tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTa                                            <  1:957888/53‑1 (MQ=255)                                            tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                          >  1:1891808/1‑71 (MQ=255)                                            tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                          >  1:3282933/1‑71 (MQ=255)                                            tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                          >  1:1004755/1‑71 (MQ=255)                                            tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                          >  1:1857135/1‑71 (MQ=255)                                            tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                          >  1:3287211/1‑71 (MQ=255)                                             aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGGGCGCTTATTgg                         >  1:2246104/1‑71 (MQ=255)                                             aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                         >  1:3619948/1‑71 (MQ=255)                                             aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                         >  1:3656216/1‑71 (MQ=255)                                             aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                         >  1:3205421/1‑71 (MQ=255)                                             aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                         >  1:472771/1‑71 (MQ=255)                                             aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                         >  1:2117673/1‑71 (MQ=255)                                                      aTCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGa                   >  1:1101784/1‑68 (MQ=255)                                                       tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa               <  1:2859813/71‑1 (MQ=255)                                                       tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa               <  1:1338074/71‑1 (MQ=255)                                                        cGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa               <  1:470352/70‑1 (MQ=255)                                                         gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                       <  1:1042718/45‑1 (MQ=18)                                                         gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                       <  1:3324596/45‑1 (MQ=18)                                                         gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                       <  1:1601209/45‑1 (MQ=18)                                                         gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                       <  1:113402/45‑1 (MQ=18)                                                            cGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                         <  1:1999930/40‑1 (MQ=255)                                                             gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc         <  1:2189975/71‑1 (MQ=255)                                                             gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc         >  1:1735470/1‑71 (MQ=255)                                                              tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                       >  1:1805924/1‑40 (MQ=255)                                                              tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc        >  1:238452/1‑71 (MQ=255)                                                              tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc        <  1:2043933/71‑1 (MQ=255)                                                                 tCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                         >  1:2978945/1‑35 (MQ=255)                                                                   gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGc                 <  1:1708764/57‑1 (MQ=255)                                                                   gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGc                 <  1:3155017/57‑1 (MQ=255)                                                                   gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTc   <  1:2489568/71‑1 (MQ=255)                                                                    ccATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTc   >  1:1117393/1‑70 (MQ=255)                                                                    ccATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCa  >  1:3390240/1‑71 (MQ=255)                                                                     |                                                                      ATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCA  >  minE/1232292‑1232429

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