Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	458,679	G→A	24.2%	G73S (GGC→AGC)	nei →	endonuclease VIII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	458,679	0	G	A	24.2%	116.4 / 35.4	62	G73S (GGC→AGC)	nei	endonuclease VIII
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base A (9/6);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

ACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGTATTGCGGGTAAAACT  >  minE/458610‑458743                                                                      |                                                                 aCGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACg                                                                  >  1:4078035/1‑70 (MQ=255)  cgcgTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGc                                                                <  1:485015/71‑1 (MQ=255)  cgcgTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGc                                                                <  1:786279/71‑1 (MQ=255)          aaGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGgcg                                                        >  1:584805/1‑71 (MQ=255)             gCGTTGTTAAGTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGt                                                     <  1:2614537/71‑1 (MQ=255)               gttgttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTgg                                                   <  1:2609345/71‑1 (MQ=255)                ttgttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCt                                                            <  1:900660/61‑1 (MQ=255)                 tgttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACAGCGt                                                              >  1:3728603/1‑58 (MQ=255)                  gttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGcgcg                                                      <  1:1984306/65‑1 (MQ=255)                  gttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGcgcg                                                      <  1:278789/65‑1 (MQ=255)                  gttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 <  1:764531/70‑1 (MQ=255)                  gttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 <  1:726015/70‑1 (MQ=255)                  gttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 <  1:228711/70‑1 (MQ=255)                  gttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 <  1:3581037/70‑1 (MQ=255)                   ttAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGTGGTTGa                                               >  1:2857348/1‑71 (MQ=255)                    tAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGa                                               >  1:3300484/1‑70 (MQ=255)                     aaCTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTc                                                             >  1:817490/1‑55 (MQ=255)                     aaCTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGAt                                              <  1:4017845/70‑1 (MQ=255)                      aCTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGt                                                              >  1:3243925/1‑53 (MQ=255)                       cTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 <  1:689621/65‑1 (MQ=255)                       cTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATAcc                                           >  1:2534678/1‑71 (MQ=255)                       cTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATAcc                                           <  1:1519788/71‑1 (MQ=255)                         cATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGgcgc                                                       >  1:190566/1‑57 (MQ=255)                         cATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGTGGTTGATAc                                            >  1:1159373/1‑68 (MQ=255)                          aTTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGc                                        >  1:1912568/1‑71 (MQ=255)                           tttttCCAACGATTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 <  1:1383588/61‑1 (MQ=255)                             tttCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGa                                               >  1:2032326/1‑61 (MQ=255)                             tttCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATa                                             >  1:3870192/1‑63 (MQ=255)                             tttCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGTGGTTGATAc                                            >  1:4092873/1‑64 (MQ=255)                              ttCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCg                                       >  1:285633/1‑68 (MQ=255)                                   aCGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATAc                                            >  1:91789/1‑58 (MQ=255)                                   aCGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCg                               >  1:2700318/1‑71 (MQ=255)                                   aCGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCg                               >  1:2100371/1‑71 (MQ=255)                                   aCGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCg                               >  1:3390973/1‑71 (MQ=255)                                      aCTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 >  1:2153801/1‑50 (MQ=255)                                      aCTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCg                               >  1:3998636/1‑68 (MQ=255)                                      aCTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAg                            <  1:684118/71‑1 (MQ=255)                                       cTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTgg                                                   <  1:3610366/47‑1 (MQ=255)                                          aaCGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGAc                          <  1:3560877/69‑1 (MQ=255)                                          aaCGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCa                        >  1:312997/1‑71 (MQ=255)                                           aCGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGa                                               >  1:2671389/1‑47 (MQ=255)                                           aCGCTATACAGCCATAATCAGCTTTACGGCTTCTTGCGCGTGGTTGAt                                              >  1:3509644/1‑48 (MQ=255)                                            cGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGtt                                                 >  1:3227923/1‑44 (MQ=255)                                            cGCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGTGGTTGa                                               <  1:546486/46‑1 (MQ=255)                                            cGCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGTGGTTGa                                               <  1:3209988/46‑1 (MQ=255)                                            cGCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGTGGTTGATACCGGCGAag                                    >  1:2581938/1‑57 (MQ=255)                                             gCTCTACAGCCATAATCAGCTTTACAGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCAcgc                     <  1:1664200/71‑1 (MQ=255)                                              ctctACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCAcgc                     >  1:1017858/1‑70 (MQ=255)                                              ctctACAGCCATAATAAGCTTTACAGCGTCTGGCGCGTGGTTGATa                                             >  1:623732/1‑46 (MQ=38)                                                   cAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAg                  <  1:582347/68‑1 (MQ=255)                                                   cAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGt                 <  1:3284561/69‑1 (MQ=255)                                                      ccATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGAt                                              >  1:2433998/1‑37 (MQ=255)                                                         tAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGa                           >  1:2197948/1‑53 (MQ=255)                                                         tAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGTATTGCgg          >  1:309368/1‑70 (MQ=255)                                                           aTCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCg                                       <  1:3733303/39‑1 (MQ=255)                                                           aTCAGCTTTACAGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCa                             >  1:2945619/1‑49 (MQ=255)                                                           aTCAGCTTTACAGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCa                             >  1:1607667/1‑49 (MQ=255)                                                           aTCAGCTTTACAGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGTATTGCGGGt        >  1:4034533/1‑70 (MQ=255)                                                            tCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGTATTGCGGGTaa      <  1:2982955/71‑1 (MQ=255)                                                            tCAGCTTTACAGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCg                               <  1:2967050/46‑1 (MQ=255)                                                            tCAGCTTTACAGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGa                           <  1:1419186/50‑1 (MQ=255)                                                              aGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGTATTGCGGGTa       >  1:3142784/1‑68 (MQ=255)                                                                cTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGTATTGCGGGTAAAACt  <  1:2343878/71‑1 (MQ=255)                                                                      |                                                                 ACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGAGTATTGCGGGTAAAACT  >  minE/458610‑458743

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