Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,232,357	C→G	9.8%	A133A (GCG→GCC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,232,357	0	C	G	9.8%	210.6 / 11.6	71	A133A (GCG→GCC)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base G (3/4);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

AACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCC  >  minE/1232287‑1232423                                                                       |                                                                   aaCGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATc                                                                    >  1:1901630/1‑71 (MQ=255)    gAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGc                                                                  >  1:822838/1‑70 (MQ=255)     aaTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATc                                                                    >  1:2973247/1‑67 (MQ=255)      aTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAt                                                               <  1:21634/71‑1 (MQ=255)      aTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAt                                                               <  1:3233465/71‑1 (MQ=255)         ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                            <  1:3652393/71‑1 (MQ=255)         ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                            <  1:1742057/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)                            gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGc                                                                  <  1:2475121/46‑1 (MQ=255)                            gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         <  1:1310724/71‑1 (MQ=255)                            gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         >  1:1637556/1‑71 (MQ=255)                            gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         <  1:2323041/71‑1 (MQ=255)                            gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         >  1:546977/1‑71 (MQ=255)                             aaaaTCTTCTGCCAGCGACATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                        >  1:2319817/1‑71 (MQ=255)                              aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                        >  1:2356153/1‑70 (MQ=255)                              aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                        <  1:1746360/70‑1 (MQ=255)                              aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                       <  1:14212/71‑1 (MQ=255)                              aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                       <  1:1673387/71‑1 (MQ=255)                                aTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAtt                                                              >  1:2485635/1‑46 (MQ=255)                                 tcttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCa                                                                >  1:2741840/1‑43 (MQ=255)                                 tcttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCa                                                                >  1:765386/1‑43 (MQ=255)                                 tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                        >  1:1842334/1‑51 (MQ=21)                                 tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                        >  1:3135235/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:355919/1‑71 (MQ=255)                                            cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                         >  1:329111/1‑71 (MQ=255)                                             gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                        <  1:3243234/71‑1 (MQ=255)                                             gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                        <  1:2452894/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:1857135/1‑71 (MQ=255)                                                 tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:3287211/1‑71 (MQ=255)                                                 tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:1891808/1‑71 (MQ=255)                                                 tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:3282933/1‑71 (MQ=255)                                                 tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:1004755/1‑71 (MQ=255)                                                  aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGGGCGCTTATTgg                   >  1:2246104/1‑71 (MQ=255)                                                  aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:3205421/1‑71 (MQ=255)                                                  aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:472771/1‑71 (MQ=255)                                                  aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:3656216/1‑71 (MQ=255)                                                  aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:2117673/1‑71 (MQ=255)                                                  aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:3619948/1‑71 (MQ=255)                                                           aTCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGa             >  1:1101784/1‑68 (MQ=255)                                                            tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa         <  1:1338074/71‑1 (MQ=255)                                                            tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa         <  1:2859813/71‑1 (MQ=255)                                                             cGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa         <  1:470352/70‑1 (MQ=255)                                                              gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:113402/45‑1 (MQ=18)                                                              gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:1042718/45‑1 (MQ=18)                                                              gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:3324596/45‑1 (MQ=18)                                                              gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:1601209/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:2043933/71‑1 (MQ=255)                                                                   tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc  >  1:238452/1‑71 (MQ=255)                                                                      tCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                   >  1:2978945/1‑35 (MQ=255)                                                                       |                                                                   AACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCC  >  minE/1232287‑1232423

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