Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	854,891	A→G	20.6%	F33F (TTT→TTC)	yciG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	854,891	0	A	G	20.6%	130.9 / 28.8	68	F33F (TTT→TTC)	yciG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base G (6/8);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

TGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGA  >  Exported/854822‑854958                                                                      |                                                                    tGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                     >  1:1199558/1‑70 (MQ=255)   aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                                    <  1:4841389/69‑1 (MQ=255)      tattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                               >  1:4708965/1‑71 (MQ=255)        ttaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                              <  1:4865645/70‑1 (MQ=255)           ccaccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgc                                                           <  1:1450435/70‑1 (MQ=255)             accGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                                                        <  1:4981776/71‑1 (MQ=255)               cGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                      <  1:2394448/71‑1 (MQ=255)                 cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATg                                                    >  1:1814136/1‑71 (MQ=255)                    tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCt                                                  >  1:599170/1‑70 (MQ=255)                       ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTAccccc                                                          >  1:2537712/3‑59 (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:2408153/70‑2 (MQ=255)                         accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                             <  1:1115050/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)                                tttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                                                 <  1:3374548/43‑1 (MQ=255)                                 ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTa                                                               >  1:223896/1‑44 (MQ=255)                                 ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                      <  1:1200742/69‑3 (MQ=255)                                 ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                      <  1:3709759/69‑3 (MQ=255)                                 ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                      <  1:2860583/69‑3 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                                                  >  1:2976150/1‑41 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                     >  1:3626140/1‑38 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                     >  1:2882276/1‑38 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                               >  1:3777310/1‑44 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                               >  1:4983566/1‑44 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                               >  1:4340307/1‑44 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                               >  1:4627231/1‑44 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                              >  1:1949034/1‑45 (MQ=255)                                 ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                              >  1:4608047/1‑45 (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:4348405/71‑1 (MQ=255)                                                             aTCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt        <  1:4524946/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:136047/1‑62 (MQ=255)                                                                   tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt           >  1:1497894/1‑62 (MQ=255)                                                                   tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt           >  1:1160322/1‑62 (MQ=255)                                                                   tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt           >  1:679450/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)                                                                      |                                                                    TGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGA  >  Exported/854822‑854958

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