Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	854,861	C→T	10.2%	A43A (GCG→GCA)	yciG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	854,861	0	C	T	10.2%	170.2 / 8.7	59	A43A (GCG→GCA)	yciG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/23);  new base T (2/4);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGG  >  minE/854791‑854924                                                                       |                                                                atTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAccc                                                                 >  1:346464/1‑71 (MQ=255) atTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAccc                                                                 >  1:151715/1‑71 (MQ=255)  ttAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCg                                                                >  1:513926/1‑71 (MQ=255)    aataatTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCt                                                              >  1:651470/1‑71 (MQ=255)       aatTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                                             >  1:890762/1‑69 (MQ=255)        atTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                                             >  1:716124/1‑68 (MQ=255)        atTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                           <  1:1857160/70‑1 (MQ=255)              ggTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                                                    <  1:1114954/71‑1 (MQ=255)               gTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTc                                                            <  1:1045767/62‑1 (MQ=255)               gTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTc                                                            <  1:2331677/62‑1 (MQ=255)                  aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAccc                                                                 <  1:4013686/54‑1 (MQ=255)                  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)                     cAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                                                    >  1:719115/1‑64 (MQ=255)                      aGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATg                                                       >  1:1518414/1‑60 (MQ=255)                      aGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                                            <  1:2826311/71‑1 (MQ=255)                        gATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                           >  1:2523876/1‑54 (MQ=255)                          tttGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTg                                                >  1:1619656/1‑63 (MQ=255)                            tGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                                         >  1:1166275/1‑52 (MQ=255)                               cTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                           <  1:832163/47‑1 (MQ=255)                               cTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                                         >  1:2209277/1‑49 (MQ=255)                                 gATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGcgcg                                                   <  1:2607053/53‑1 (MQ=255)                                 gATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                  >  1:3311015/1‑70 (MQ=255)                                  aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                <  1:1728286/71‑1 (MQ=255)                                  aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                >  1:2608063/1‑71 (MQ=255)                                  aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                  <  1:1100104/69‑1 (MQ=255)                                  aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                  <  1:3055749/69‑1 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                                         >  1:171145/1‑45 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                           >  1:2845877/1‑59 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                           >  1:1897317/1‑59 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                               >  1:2845390/1‑71 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                               <  1:1195799/71‑1 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                               >  1:2621739/1‑71 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                               >  1:880963/1‑71 (MQ=255)                                   tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                               <  1:775819/71‑1 (MQ=255)                                   cttattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                           <  1:1845132/42‑1 (MQ=255)                                    ttattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                               >  1:2959527/1‑70 (MQ=255)                                    ttattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                               >  1:3186597/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:976829/1‑63 (MQ=255)                                                cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccg                          >  1:3628490/1‑63 (MQ=255)                                                   tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCa                                         >  1:2943820/1‑45 (MQ=255)                                                   tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                  <  1:3636562/52‑1 (MQ=255)                                                      ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                               >  1:3810581/3‑52 (MQ=255)                                                      ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                               >  1:1305664/3‑52 (MQ=255)                                                           gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt       <  1:1599748/71‑1 (MQ=255)                                                            ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt      <  1:4034793/71‑1 (MQ=255)                                                            ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt      <  1:3075387/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:335368/56‑1 (MQ=255)                                                              ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTAt                   <  1:2961839/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)                                                                       |                                                                ATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGG  >  minE/854791‑854924

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