Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	854,912	A→T	13.2%	G26G (GGT→GGA)	yciG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	854,912	0	A	T	13.2%	171.2 / 17.6	68	G26G (GGT→GGA)	yciG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base T (3/6);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

GACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTT  >  Exported/854845‑854982                                                                    |                                                                       taccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                                     >  1:5201179/2‑70 (MQ=255)  accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTccac                                                                    >  1:2807254/1‑69 (MQ=255)  accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                                     <  1:1115050/70‑2 (MQ=255)  accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                                     <  1:2414555/70‑2 (MQ=255)  accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                                     <  1:2408153/70‑2 (MQ=255)  aCCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGAcctc                                                                    >  1:574382/1‑69 (MQ=255)     gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGTATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                                 <  1:3742697/71‑1 (MQ=255)     gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                                 <  1:49739/71‑1 (MQ=255)     gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                                 <  1:3456577/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:3709759/69‑3 (MQ=255)          ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                              <  1:1200742/69‑3 (MQ=255)          ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                              <  1:2860583/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)                      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:2239612/1‑71 (MQ=255)                                        cATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCg                              >  1:4987754/1‑71 (MQ=255)                                         aTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                                                 >  1:122792/1‑51 (MQ=255)                                         aTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCgg                             <  1:2795517/71‑1 (MQ=255)                                          tttttAAAATTACCGCCGCTATGCTGACCGCCttt                                                                <  1:1575907/35‑1 (MQ=255)                                          tttttAAAATTACCGCCGCTATGCTGACCGCCttt                                                                <  1:1184598/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:1160322/1‑62 (MQ=255)                                            tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                                   >  1:679450/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:1108846/70‑1 (MQ=255)                                             ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                          <  1:3303030/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:2560812/1‑62 (MQ=255)                                                              aTGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTg                 >  1:1894013/1‑62 (MQ=255)                                                               tGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                          <  1:2061437/52‑1 (MQ=255)                                                               tGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                          <  1:715450/52‑1 (MQ=255)                                                               tGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                          <  1:4536748/52‑1 (MQ=255)                                                               tGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                          <  1:4342678/52‑1 (MQ=255)                                                                 cTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa               >  1:76538/1‑61 (MQ=255)                                                                 cTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa               >  1:992671/1‑61 (MQ=255)                                                                 cTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa               >  1:4916590/1‑61 (MQ=255)                                                                 cTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa               >  1:455601/1‑61 (MQ=255)                                                                 cTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa               >  1:1712787/1‑61 (MQ=255)                                                                  tGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATg    >  1:3840097/1‑71 (MQ=255)                                                                  tGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATg    >  1:696905/1‑71 (MQ=255)                                                                  tGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATg    >  1:1600322/1‑71 (MQ=255)                                                                   gACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAtt                      <  1:3047110/52‑1 (MQ=255)                                                                   gACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAtt                      <  1:2775555/52‑1 (MQ=255)                                                                    aCCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGtt  <  1:2378291/71‑1 (MQ=255)                                                                    aCCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGtt  >  1:17026/1‑71 (MQ=255)                                                                    |                                                                       GACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTT  >  Exported/854845‑854982

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