Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,432,788	C→T	15.1%	G22E (GGG→GAG)	trxA ←	thioredoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,432,788	0	C	T	15.1%	211.0 / 25.1	86	G22E (GGG→GAG)	trxA	thioredoxin 1
Reads supporting (aligned to +/- strand):  ref base C (34/39);  new base T (7/6);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

CATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCTCATATATAA  >  minE/2432720‑2432858                                                                     |                                                                       cATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcccc                                                                       <  1:3183243/70‑1 (MQ=255)     cAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCg                                                                  >  1:2286162/1‑71 (MQ=255)     cAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCg                                                                  >  1:1320161/1‑71 (MQ=255)      aGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGc                                                                 >  1:3518101/1‑71 (MQ=255)       gAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCGCt                                                                >  1:56577/1‑71 (MQ=255)       gAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCt                                                                <  1:3627464/71‑1 (MQ=255)         aTCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCGCttt                                                              <  1:110019/71‑1 (MQ=255)         aTCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCttt                                                              <  1:2813693/71‑1 (MQ=255)         aTCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCttt                                                              <  1:164716/71‑1 (MQ=255)          tCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCttt                                                              <  1:2277486/70‑1 (MQ=255)          tCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTg                                                             <  1:1113035/71‑1 (MQ=255)          tCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTg                                                             <  1:3224031/71‑1 (MQ=255)            ggggCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAg                                                           <  1:3934709/71‑1 (MQ=255)                   tCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATcc                                                    <  1:1541362/71‑1 (MQ=255)                    cATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATc                                                     <  1:3405650/69‑1 (MQ=255)                    cATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCg                                                   >  1:3667789/1‑71 (MQ=255)                     aTTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt                                                  <  1:1679883/71‑1 (MQ=255)                     aTTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt                                                  <  1:415565/71‑1 (MQ=255)                     aTTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCg                                                   <  1:451684/70‑1 (MQ=255)                      ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgtg                                                 >  1:352178/1‑71 (MQ=255)                      ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt                                                  <  1:1234217/70‑1 (MQ=255)                       tttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTc                                                                    >  1:2194115/1‑51 (MQ=255)                        ttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAg                                                           >  1:219135/1‑59 (MQ=255)                        ttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgtgt                                                <  1:2913881/70‑1 (MQ=255)                              ggACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACt                                         <  1:1439229/71‑1 (MQ=255)                              ggACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACt                                         >  1:2396577/1‑71 (MQ=255)                              ggACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACt                                         >  1:307726/1‑71 (MQ=255)                              ggACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACt                                         >  1:3621382/1‑71 (MQ=255)                               gACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACTg                                        <  1:3122807/71‑1 (MQ=255)                               gACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACt                                         >  1:2180197/1‑70 (MQ=255)                                aCCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTg                                        <  1:2716693/70‑1 (MQ=255)                                 ccGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTgtc                                      <  1:3947535/71‑1 (MQ=255)                                  cGCACCACTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCGCTTTGAg                                                           >  1:1201484/1‑49 (MQ=255)                                  cGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTg                                        <  1:1791920/68‑1 (MQ=255)                                   gCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTgtcgt                                    <  1:2504594/71‑1 (MQ=255)                                   gCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTgtcgt                                    <  1:1886595/71‑1 (MQ=255)                                    caccacTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACTgtc                                      >  1:3088676/1‑68 (MQ=255)                                    caccacTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATACGTGTCAAAAc                                          >  1:591322/1‑64 (MQ=255)                                      ccacTCTGCCCAGAAATCGACGAGGATCGCCTCGt                                                                     >  1:1959428/1‑35 (MQ=255)                                      ccacTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCgtca                                  >  1:2339349/1‑70 (MQ=255)                                       cacTCTGCCCAGAAATCGACGAGGATCGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCgtcagt                                <  1:630275/71‑1 (MQ=255)                                       cacTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCgtcagt                                <  1:2193303/71‑1 (MQ=255)                                        acTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCgtcagt                                >  1:34964/1‑70 (MQ=255)                                         ctctGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa                                                            >  1:3177113/1‑41 (MQ=255)                                         ctctGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACt                                         <  1:66172/60‑1 (MQ=255)                                         ctctGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTgtc                                      <  1:1524369/63‑1 (MQ=255)                                         ctctGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCgtcagt                                >  1:2773855/1‑69 (MQ=255)                                           ctGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa                                                            >  1:147780/1‑39 (MQ=255)                                           ctGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa                                                            >  1:3858753/1‑39 (MQ=255)                                            tGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAAc                                          >  1:566481/1‑56 (MQ=255)                                            tGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGtcagtcag                             <  1:2632998/69‑1 (MQ=255)                                              cccAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTg                          <  1:1490082/70‑1 (MQ=255)                                               ccAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaa                        <  1:2030196/71‑1 (MQ=255)                                               ccAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaa                        <  1:2037407/71‑1 (MQ=255)                                                cAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACAt                                                      >  1:2067227/1‑40 (MQ=255)                                                cAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt                                                  <  1:2649358/44‑1 (MQ=255)                                                cAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaat                       <  1:2327585/71‑1 (MQ=255)                                                cAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaat                       <  1:3883965/71‑1 (MQ=255)                                                 aGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGtcagtcag                             >  1:1269233/1‑64 (MQ=255)                                                 aGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaat                       >  1:3028670/1‑70 (MQ=255)                                                  gAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaataa                     <  1:657029/71‑1 (MQ=255)                                                  gAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaataa                     >  1:757873/1‑71 (MQ=255)                                                    aaTCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCgtcagtca                              <  1:3823592/60‑1 (MQ=255)                                                     aTCGACGAGGATCGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAAtt                   <  1:3717133/70‑1 (MQ=255)                                                     aTCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAgg                            >  1:567532/1‑61 (MQ=255)                                                     aTCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAgg                            >  1:949154/1‑61 (MQ=255)                                                     aTCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAAttt                  <  1:2122435/71‑1 (MQ=255)                                                     aTCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAAttt                  <  1:2225040/71‑1 (MQ=255)                                                      tCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAgg                            >  1:505626/1‑60 (MQ=255)                                                       cgacgaGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGa                         >  1:3716076/1‑62 (MQ=255)                                                       cgacgaGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGa                         >  1:1673178/1‑62 (MQ=255)                                                        gacgaGGATCGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTAt               <  1:909530/71‑1 (MQ=255)                                                        gacgaGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTgtcgtc                                   <  1:2339297/51‑1 (MQ=255)                                                        gacgaGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTAt               >  1:1567172/1‑71 (MQ=255)                                                           gaGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaataat                    <  1:3473829/63‑1 (MQ=255)                                                           gaGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGc            >  1:1663207/1‑71 (MQ=255)                                                           gaGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGc            >  1:3293197/1‑71 (MQ=255)                                                            aGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTgtcgt                                    <  1:418614/46‑1 (MQ=255)                                                            aGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGc            >  1:1668386/1‑70 (MQ=255)                                                            aGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCt           <  1:3514493/71‑1 (MQ=255)                                                             ggATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCgtcagtc                               <  1:2077625/50‑1 (MQ=255)                                                             ggATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCt           >  1:552406/1‑70 (MQ=255)                                                             ggATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCt           >  1:2706784/1‑70 (MQ=255)                                                              gATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCTc          <  1:2036471/70‑1 (MQ=255)                                                                 cGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTAt               >  1:73609/1‑62 (MQ=255)                                                                 cGCCTCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTAt               >  1:943009/1‑62 (MQ=255)                                                                  gCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGaataa                     >  1:3283269/1‑55 (MQ=255)                                                                   ccccGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCTCatata     >  1:965828/1‑70 (MQ=255)                                                                     ccGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCTCATATATaa  >  1:3335747/1‑71 (MQ=255)                                                                     |                                                                       CATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGTGTCAAAACTGTCGTCAGTCAGGTGAATAATTTTATCGCTCATATATAA  >  minE/2432720‑2432858

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