Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,371,634	C→T	12.8%	E151E (GAG→GAA)	fre ←	flavin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,371,634	0	C	T	12.8%	209.7 / 15.8	78	E151E (GAG→GAA)	fre	flavin reductase
Reads supporting (aligned to +/- strand):  ref base C (41/27);  new base T (5/5);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGTTAC  >  minE/2371565‑2371704                                                                      |                                                                       gCCCGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCt                                                                       <  1:1768755/71‑1 (MQ=255) gCCCGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCt                                                                       <  1:2350966/71‑1 (MQ=255)   ccGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCg                                                                     <  1:2861332/71‑1 (MQ=255)   ccGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCg                                                                     <  1:1971796/71‑1 (MQ=255)     gCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCg                                                                     >  1:1846764/1‑69 (MQ=255)     gCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCa                                                                   <  1:2981758/71‑1 (MQ=255)         ccGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCagag                                                                <  1:2081176/70‑1 (MQ=255)         ccGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCa                                                                   >  1:304392/1‑67 (MQ=255)         ccGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCa                                                                   >  1:1686592/1‑67 (MQ=255)            gttgttCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCa                                                            <  1:3839467/71‑1 (MQ=255)             ttgttCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCAt                                                           <  1:853198/71‑1 (MQ=255)              tgttCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAgaga                                                               >  1:606244/1‑66 (MQ=255)                ttCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGAt                                                              >  1:479455/1‑65 (MQ=255)                 tCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCa                                                                   >  1:3846268/1‑59 (MQ=255)                 tCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATa                                                          >  1:2821495/1‑68 (MQ=255)                   aaCCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATAc                                                         <  1:3923741/67‑1 (MQ=255)                   aaCCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGa                                                      <  1:787465/70‑1 (MQ=255)                     ccaccGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGAtgc                                                   <  1:383486/71‑1 (MQ=255)                     ccaccGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGAtgc                                                   >  1:2797055/1‑71 (MQ=255)                     ccaccGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGAtgc                                                   <  1:1983342/71‑1 (MQ=255)                           gCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGctct                                             <  1:1654927/71‑1 (MQ=255)                           gCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGctct                                             <  1:309184/71‑1 (MQ=255)                            caccacTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCtt                                            >  1:337927/1‑71 (MQ=255)                            caccacTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCtt                                            <  1:2383512/71‑1 (MQ=255)                             accacTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTc                                           <  1:1252092/71‑1 (MQ=255)                             accacTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCtt                                            >  1:1237556/1‑70 (MQ=255)                               cacTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCa                                          >  1:1271079/1‑70 (MQ=255)                                 cTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgc                                       >  1:2470555/1‑71 (MQ=255)                                 cTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgc                                       >  1:2324621/1‑71 (MQ=255)                                 cTTGCAGACCAGGATGCTTCAACGAAAGCGCCGCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgc                                       >  1:2926495/1‑71 (MQ=255)                                  ttGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCAcgc                                       <  1:3751968/70‑1 (MQ=255)                                  ttGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgcc                                      >  1:1552928/1‑71 (MQ=255)                                   tGCAGACCAGGATGCTTCAACGAAAGTGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgccc                                     >  1:3563579/1‑71 (MQ=255)                                   tGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgccc                                     >  1:2261366/1‑71 (MQ=255)                                    gCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGAtgctgc                                                >  1:368196/1‑59 (MQ=255)                                    gCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGAtg                                                    >  1:3911913/1‑55 (MQ=255)                                    gCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGctct                                             >  1:2206245/1‑62 (MQ=255)                                    gCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCAcgcc                                      <  1:1276584/69‑1 (MQ=255)                                     cAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCAcgcccgc                                   >  1:568914/1‑71 (MQ=255)                                     cAGACCAGGATGCTTCAACGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCAcgcccgc                                   >  1:2900552/1‑71 (MQ=255)                                     cAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCa                                          >  1:3692273/1‑64 (MQ=255)                                       gACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACgcccgccc                                 >  1:790300/1‑71 (MQ=255)                                         ccAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACgcccgccc                                 >  1:819602/1‑69 (MQ=255)                                            ggATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGt                             >  1:1873094/1‑70 (MQ=255)                                            ggATCCTTCAACGAAAGCGCCTCAACTTCGCAGAGATc                                                             >  1:281299/1‑38 (MQ=25)                                              aTGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTaaa                          >  1:801068/1‑71 (MQ=255)                                               tGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATAc                                                         <  1:1214756/39‑1 (MQ=255)                                               tGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAAt                         <  1:3397703/71‑1 (MQ=255)                                                  ttCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGctct                                             <  1:2881020/48‑1 (MQ=255)                                                  ttCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCtt                                            <  1:1701744/49‑1 (MQ=255)                                                    cAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGa                    >  1:1167203/1‑71 (MQ=255)                                                    cAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATCCTGCTCTTCACGCCCGCCCCAGTAAATGGTGa                    >  1:3561391/1‑71 (MQ=255)                                                     aaCGAAAGCGCCTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCa                                          >  1:1475082/1‑48 (MQ=255)                                                     aaCGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATgctgct                                               <  1:431439/43‑1 (MQ=255)                                                     aaCGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATgctgct                                               <  1:3427229/43‑1 (MQ=255)                                                      aCGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTc                                           >  1:1475830/1‑46 (MQ=255)                                                      aCGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGata                  <  1:449322/71‑1 (MQ=255)                                                        gAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTa                            <  1:2185467/59‑1 (MQ=255)                                                         aaaGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCa                                          >  1:2236215/1‑44 (MQ=255)                                                          aaGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACgcccgccc                                 >  1:591700/1‑52 (MQ=255)                                                           aGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGata                  >  1:1459212/1‑66 (MQ=255)                                                            gcgcCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTg                     <  1:2106395/62‑1 (MQ=255)                                                             cgcCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGt                      <  1:3267886/60‑1 (MQ=255)                                                               ccTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCa                               >  1:802245/1‑49 (MQ=255)                                                               ccTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGttt         >  1:2557319/1‑71 (MQ=255)                                                               acTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCa               >  1:2052049/2‑65 (MQ=255)                                                                cTCAAGTTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGata                  <  1:1897461/61‑1 (MQ=255)                                                                cTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGtt          <  1:1270201/69‑1 (MQ=255)                                                                cTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGtt          <  1:1680704/69‑1 (MQ=255)                                                                cTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTg        >  1:1939014/1‑71 (MQ=255)                                                                 tCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTgg       >  1:695699/1‑71 (MQ=255)                                                                  cAAGCTCGCAGAGATCATAGAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTggg      >  1:2431088/1‑71 (MQ=255)                                                                   caGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGt     >  1:3728545/2‑71 (MQ=255)                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGt                      >  1:2179763/1‑54 (MQ=255)                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGat                   >  1:3815075/1‑57 (MQ=255)                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGat                   >  1:3371182/1‑57 (MQ=255)                                                                   aaGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTgg       <  1:1282961/69‑1 (MQ=255)                                                                    aGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGt     >  1:1857383/1‑70 (MQ=255)                                                                      cTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGTTAc  <  1:3427/71‑1 (MQ=255)                                                                      cTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGTTAc  <  1:1106691/71‑1 (MQ=255)                                                                      |                                                                       GCCCGCTTCCGGTTGTTCAACCACCGGCACCACTTGCAGACCAGGATGCTTCAACGAAAGCGCCTCAAGCTCGCAGAGATCATACAGATGCTGCTCTTCACGCCCGCCCCAGTAAATGGTGATATCACGGTTTGGGTTAC  >  minE/2371565‑2371704

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