Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,704,553	G→A	9.8%	R255H (CGC→CAC)	damX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,704,553	0	G	A	9.8%	262.8 / 12.9	91	R255H (CGC→CAC)	damX	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/43);  new base A (6/3);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGTTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGAAGCCG  >  minE/2704483‑2704623                                                                       |                                                                       ccTGTTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTAcgcg                                                                        >  1:3130811/1‑71 (MQ=255)  cTGTTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACgcgc                                                                       >  1:1384545/1‑71 (MQ=255)   tGTTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGcc                                                                      >  1:1629920/1‑71 (MQ=255)   tGTTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGcc                                                                      >  1:2013202/1‑71 (MQ=255)     ttCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGcc                                                                      <  1:1148317/69‑1 (MQ=255)        gcaatggcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCt                                                                 >  1:1309809/1‑71 (MQ=255)         caatggcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTc                                                                <  1:3549222/71‑1 (MQ=255)           atggcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGt                                                              >  1:3675136/1‑71 (MQ=255)           atggcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGt                                                              <  1:589670/71‑1 (MQ=255)           atggcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGt                                                              <  1:1046289/71‑1 (MQ=255)              gcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGt                                                              >  1:630315/1‑68 (MQ=255)              gcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcag                                                           >  1:3602974/1‑71 (MQ=255)               caatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCg                                                               <  1:2126152/66‑1 (MQ=255)                aatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcagca                                                         <  1:3897415/71‑1 (MQ=255)                  tgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcagcag                                                        >  1:3898570/1‑70 (MQ=255)                    cATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGc                                                      <  1:1776595/70‑1 (MQ=255)                     aTCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTc                                                                >  1:2502585/1‑59 (MQ=255)                     aTCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTc                                                                >  1:2253924/1‑59 (MQ=255)                     aTCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCgg                                                    <  1:1005726/71‑1 (MQ=255)                     aTCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCgg                                                    <  1:155314/71‑1 (MQ=255)                     aTCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCgg                                                    <  1:1284300/71‑1 (MQ=255)                      tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcagca                                                         <  1:1936874/65‑1 (MQ=255)                      tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGt                                                   >  1:2300817/1‑71 (MQ=255)                      tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGt                                                   <  1:2888753/71‑1 (MQ=255)                       cgcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcag                                                           >  1:1295142/1‑62 (MQ=255)                       cgcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGt                                                   >  1:713874/1‑70 (MQ=255)                        gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGa                                                 >  1:1690083/1‑71 (MQ=255)                           tGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                              <  1:2719396/71‑1 (MQ=255)                           tGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                              <  1:2724033/71‑1 (MQ=255)                            gACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGAtt                                               >  1:2615797/1‑69 (MQ=255)                            gACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                              >  1:1921845/1‑70 (MQ=255)                            gACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                              >  1:2418540/1‑70 (MQ=255)                             acacGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGaa                                            >  1:560660/1‑71 (MQ=255)                             acacGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGa                                             <  1:887609/70‑1 (MQ=255)                             acacGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGa                                             <  1:2121115/70‑1 (MQ=255)                               acGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcagca                                                         >  1:1580840/1‑56 (MQ=255)                               acGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                              <  1:3574860/67‑1 (MQ=255)                               acGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAAcc                                          <  1:2532842/71‑1 (MQ=255)                                 ggCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGt                                                              <  1:3452132/49‑1 (MQ=255)                                 ggCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCg                                         >  1:1332362/1‑70 (MQ=255)                                  gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCgg                                                    >  1:2789624/1‑58 (MQ=255)                                  gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaa                                       <  1:3228063/71‑1 (MQ=255)                                  gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaa                                       <  1:1235529/71‑1 (MQ=255)                                  gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGa                                        >  1:364601/1‑70 (MQ=255)                                   cGAAAACGCAAACCGCTGAACTTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaa                                       >  1:2224644/1‑70 (MQ=255)                                   cGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                     >  1:1478773/1‑40 (MQ=255)                                    gAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaa                                     <  1:2508859/71‑1 (MQ=255)                                    gAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcagcag                                                        >  1:3738163/1‑52 (MQ=255)                                    gAAAACGCAAACCGCTGAACCTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaa                                     <  1:243376/71‑1 (MQ=255)                                     aaaaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaa                                     >  1:2635660/1‑70 (MQ=255)                                      aaaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                                   <  1:3341301/71‑1 (MQ=255)                                      aaaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                                   <  1:610725/71‑1 (MQ=255)                                      aaaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                                   <  1:2106460/71‑1 (MQ=255)                                       aaCGCAAACCGCTGAACGTCCGTCCACTACGCGcc                                                                      >  1:1694391/1‑35 (MQ=255)                                       aaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcagcag                                                        >  1:8804/1‑49 (MQ=255)                                       aaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcagcag                                                        >  1:1695466/1‑49 (MQ=255)                                        aCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                                   >  1:1685742/1‑69 (MQ=255)                                        aCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAAcc                                 <  1:2514803/71‑1 (MQ=255)                                        aCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAAcc                                 <  1:1749328/71‑1 (MQ=255)                                        aCGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAAc                                  >  1:426135/1‑70 (MQ=255)                                           cAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgca                              <  1:3267144/71‑1 (MQ=255)                                            aaaCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCg                                <  1:2187038/68‑1 (MQ=255)                                            aaaCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCg                                <  1:1771405/68‑1 (MQ=255)                                              aCCGCTTAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaag                            >  1:2626826/1‑70 (MQ=255)                                              aCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGa                                                 >  1:368525/1‑49 (MQ=255)                                              aCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaagc                           <  1:203429/71‑1 (MQ=255)                                              aCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaagc                           <  1:343795/71‑1 (MQ=255)                                              aCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaag                            >  1:1155436/1‑70 (MQ=255)                                              aCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaag                            >  1:3084596/1‑70 (MQ=255)                                              aCCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaag                            >  1:3728185/1‑70 (MQ=255)                                              aCCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaag                            >  1:1034951/1‑70 (MQ=255)                                               ccGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaagc                           >  1:1858440/1‑70 (MQ=255)                                                cGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaagcaa                         <  1:2927211/71‑1 (MQ=255)                                                  cTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGa                                             <  1:2523193/49‑1 (MQ=255)                                                    gAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaa                                       <  1:572561/53‑1 (MQ=255)                                                     aaCGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCg                      <  1:1405770/69‑1 (MQ=255)                                                     aaCGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCg                      <  1:3047069/69‑1 (MQ=255)                                                       cgtccgtcCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGa                   >  1:1193962/1‑70 (MQ=255)                                                        gtccgtccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGaaa                 <  1:3558224/71‑1 (MQ=255)                                                        gtccgtccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGaaa                 <  1:2348908/71‑1 (MQ=255)                                                        gtccgtcaACTACGCGCCCAGCTCGTCAGCAGGCGGGGATTgc                                             >  1:3775631/1‑42 (MQ=255)                                                         tccgtccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGaaa                 <  1:3708627/70‑1 (MQ=255)                                                          ccgtccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAAc               >  1:433055/1‑71 (MQ=255)                                                           cgtccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACg              >  1:1783924/1‑71 (MQ=255)                                                           cgtccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACg              >  1:1833063/1‑71 (MQ=255)                                                             tccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACg              <  1:1559000/69‑1 (MQ=255)                                                               cACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAg           <  1:21073/70‑1 (MQ=255)                                                               cACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGa                                        >  1:3830892/1‑41 (MQ=255)                                                                aCTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaa                                     <  1:3367853/43‑1 (MQ=255)                                                                 cTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAAc               <  1:2353186/64‑1 (MQ=255)                                                                 cTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGcc         <  1:2480337/70‑1 (MQ=255)                                                                  tACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                                   >  1:1211020/1‑43 (MQ=255)                                                                   aCGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCg        >  1:1874270/1‑69 (MQ=255)                                                                    cgcgCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGa       <  1:1323475/69‑1 (MQ=255)                                                                    cgcaCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                                   >  1:3783461/1‑41 (MQ=255)                                                                    cgcaCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                                   >  1:143974/1‑41 (MQ=255)                                                                      cgcCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGAAGcc   <  1:1848840/71‑1 (MQ=255)                                                                       gcCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGc          >  1:2918226/1‑63 (MQ=255)                                                                       gcCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGAAGCCg  >  1:1220404/1‑71 (MQ=255)                                                                       |                                                                       CCTGTTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGAAGCCG  >  minE/2704483‑2704623

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