Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	854,906	A→G	20.0%	H28H (CAT→CAC)	yciG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	854,906	0	A	G	20.0%	105.1 / 22.2	55	H28H (CAT→CAC)	yciG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base G (4/7);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCAC  >  minE/854836‑854973                                                                       |                                                                    cGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                                     <  1:2394448/71‑1 (MQ=255)   cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATg                                                                   >  1:1814136/1‑71 (MQ=255)      tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCt                                                                 >  1:599170/1‑70 (MQ=255)          taccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            >  1:5201179/2‑70 (MQ=255)          gACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                                     <  1:262660/62‑1 (MQ=255)           accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTgtc                                                              <  1:2902977/68‑3 (MQ=255)           accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTccac                                                           >  1:2807254/1‑69 (MQ=255)           accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            <  1:1115050/70‑2 (MQ=255)           accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            <  1:2408153/70‑2 (MQ=255)           accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            <  1:2414555/70‑2 (MQ=255)           aCCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGAcctc                                                           >  1:574382/1‑69 (MQ=255)              gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGTATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:3742697/71‑1 (MQ=255)              gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:3456577/71‑1 (MQ=255)              gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:49739/71‑1 (MQ=255)               ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt                                                       >  1:504233/1‑71 (MQ=255)                cTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCttt                                                       <  1:1081837/70‑1 (MQ=255)                 ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGc                                                           <  1:3452642/65‑1 (MQ=255)                 ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTAc                                                     >  1:226188/1‑71 (MQ=255)                   ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                     <  1:2860583/69‑3 (MQ=255)                   ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                     <  1:3709759/69‑3 (MQ=255)                   ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                     <  1:1200742/69‑3 (MQ=255)                   ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACgg                                                   >  1:1301709/1‑71 (MQ=255)                      tACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:1418166/63‑1 (MQ=255)                       aCCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCt                                                <  1:2288744/70‑1 (MQ=255)                          cGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                                     >  1:276554/1‑46 (MQ=255)                             ttCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCt                                                                 >  1:4389781/1‑47 (MQ=255)                               cAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGAt                                       >  1:3225039/1‑71 (MQ=255)                               cAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGATATGCTGACCGCCTTTACGGCCTGCGTCGGAt                                       >  1:4576986/1‑71 (MQ=255)                                  aTGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGcc                                    <  1:2214833/71‑1 (MQ=255)                                        gTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtctaa                              <  1:4430476/71‑3 (MQ=255)                                        gTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCa                              <  1:1577361/71‑1 (MQ=255)                                             ggATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTc                         <  1:1007233/71‑1 (MQ=255)                                               aTCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt                       <  1:4524946/71‑1 (MQ=255)                                               aTCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt                       <  1:4348405/71‑1 (MQ=255)                                                tCATTCTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTc                      >  1:4203481/1‑71 (MQ=255)                                                 cATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCg                     >  1:4987754/1‑71 (MQ=255)                                                 cATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCg                     >  1:2239612/1‑71 (MQ=255)                                                  aTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                                        >  1:122792/1‑51 (MQ=255)                                                  aTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCgg                    <  1:2795517/71‑1 (MQ=255)                                                   tttttAAAATTACCGCCGCTATGCTGACCGCCttt                                                       <  1:1184598/35‑1 (MQ=255)                                                   tttttAAAATTACCGCCGCTATGCTGACCGCCttt                                                       <  1:1575907/35‑1 (MQ=255)                                                   tttttAAAATTACCGCCGCTATGCTGACCGCATTTACGGCCTGCGTCGGATGCCTTCTCAc                             <  1:5426690/61‑1 (MQ=255)                                                    ttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGc                   >  1:3910092/1‑70 (MQ=255)                                                     tttAAAATTATCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                                        <  1:4904911/48‑1 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCTCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCg                  >  1:4801684/1‑70 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCGCCTTTAc                                                     >  1:5648954/1‑35 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACTGTCTTCgg                    >  1:4192220/1‑68 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:1497894/1‑62 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:679450/1‑62 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:1160322/1‑62 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:136047/1‑62 (MQ=255)                                                     tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCg                  >  1:70114/1‑70 (MQ=255)                                                      ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:1891397/1‑61 (MQ=255)                                                      ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                 <  1:3303030/70‑1 (MQ=255)                                                      ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                 <  1:1108846/70‑1 (MQ=255)                                                                 gccgcTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAac     <  1:5376772/71‑1 (MQ=255)                                                                   cgcTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAacca   <  1:1263127/71‑1 (MQ=255)                                                                    gcTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccac  <  1:2010800/71‑1 (MQ=255)                                                                       aTGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTc           <  1:2353036/59‑1 (MQ=255)                                                                       aTGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTg        >  1:1894013/1‑62 (MQ=255)                                                                       aTGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTg        >  1:2560812/1‑62 (MQ=255)                                                                       |                                                                    CGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCAC  >  minE/854836‑854973

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