Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	854,876	T→C	11.5%	Q38Q (CAA→CAG)	yciG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	854,876	0	T	C	11.5%	128.1 / 9.8	52	Q38Q (CAA→CAG)	yciG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/25);  new base C (2/4);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

AATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCAC  >  minE/854808‑854946                                                                     |                                                                       aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGtt                                                                       >  1:857602/1‑70 (MQ=255) aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGt                                                                        >  1:2351074/1‑69 (MQ=255)   tCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCg                                                                    <  1:1569631/71‑1 (MQ=255)   tCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCg                                                                    >  1:407248/1‑71 (MQ=255)     aGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                                                                  <  1:2826311/71‑1 (MQ=255)         tttGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                                                                      >  1:1619656/1‑63 (MQ=255)                gATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                        >  1:3311015/1‑70 (MQ=255)                 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                                      >  1:2608063/1‑71 (MQ=255)                 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                                      <  1:1728286/71‑1 (MQ=255)                 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                        <  1:1100104/69‑1 (MQ=255)                 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                        <  1:3055749/69‑1 (MQ=255)                  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                                                 >  1:2845877/1‑59 (MQ=255)                  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                                                 >  1:1897317/1‑59 (MQ=255)                  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     <  1:1195799/71‑1 (MQ=255)                  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     >  1:2845390/1‑71 (MQ=255)                  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     >  1:2621739/1‑71 (MQ=255)                  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     <  1:775819/71‑1 (MQ=255)                  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     >  1:880963/1‑71 (MQ=255)                   ttattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     >  1:3186597/1‑70 (MQ=255)                   ttattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     >  1:2959527/1‑70 (MQ=255)                         ccaccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgcc                                              >  1:3602794/1‑71 (MQ=255)                           accGCTTTTTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                                            <  1:2905950/71‑1 (MQ=255)                            ccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                                            >  1:1423511/1‑70 (MQ=255)                              gCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt                                                             >  1:1558315/1‑51 (MQ=255)                              gCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                  <  1:1423181/62‑1 (MQ=255)                               cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccg                                                >  1:3628490/1‑63 (MQ=255)                               cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccg                                                >  1:976829/1‑63 (MQ=255)                                  tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCa                                                               >  1:2943820/1‑45 (MQ=255)                                  tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                        <  1:3636562/52‑1 (MQ=255)                                     ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                                     >  1:1305664/3‑52 (MQ=255)                                     ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                                     >  1:3810581/3‑52 (MQ=255)                                          gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                             <  1:1599748/71‑1 (MQ=255)                                           ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt                            <  1:3075387/71‑1 (MQ=255)                                           ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt                            <  1:4034793/71‑1 (MQ=255)                                            cTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCAtt                                                             <  1:930354/37‑1 (MQ=25)                                            cTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt                                                    <  1:3144155/46‑1 (MQ=255)                                             ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                     <  1:2300730/44‑1 (MQ=255)                                             ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTAt                                         <  1:2961839/56‑1 (MQ=255)                                             ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTAt                                         <  1:335368/56‑1 (MQ=255)                                               ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCtttttt                         <  1:1365589/70‑4 (MQ=255)                                               ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTAc                          <  1:2023585/69‑1 (MQ=255)                                               ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACgg                        >  1:3378161/1‑71 (MQ=255)                                                        cTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCg               <  1:3970837/71‑1 (MQ=255)                                                        cTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCg               <  1:64729/71‑1 (MQ=255)                                                          tCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa             >  1:854799/1‑71 (MQ=255)                                                            aGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATg           <  1:1524835/71‑1 (MQ=255)                                                             gATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTATGGCCTGCGTCGGATGc          <  1:3409029/71‑1 (MQ=255)                                                              aTGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGcc         >  1:834722/1‑71 (MQ=255)                                                                gcgcGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCg                  <  1:2461211/60‑1 (MQ=255)                                                                gcgcGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCtt       >  1:3489900/1‑71 (MQ=255)                                                                 cgcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCtt       <  1:3442763/70‑1 (MQ=255)                                                                  gcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCtt       <  1:3753653/69‑1 (MQ=255)                                                                  gcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtct     <  1:2298393/71‑1 (MQ=255)                                                                   cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtctc    <  1:3162706/71‑1 (MQ=255)                                                                     ttGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCAc  <  1:1642648/71‑1 (MQ=255)                                                                     ttGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCAc  <  1:1705479/71‑1 (MQ=255)                                                                     |                                                                       AATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCAC  >  minE/854808‑854946

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