Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,097,134	G→A	Q147* (CAG→TAG)	glpF ←	glycerol facilitator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,097,134	0	G	A	97.8%	163.2 / ‑3.7	45	Q147* (CAG→TAG)	glpF	glycerol facilitator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/23);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

TGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCG  >  REL606/4097062‑4097204                                                                         |                                                                       tGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTacac                                                                     >  1:1636971/1‑76 (MQ=255)  gCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTAcaca                                                                    <  1:716373/76‑1 (MQ=255)  gCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTAcaca                                                                    >  1:484597/1‑76 (MQ=255)     aTCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACaaaa                                                                 <  2:110957/76‑1 (MQ=255)      tCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAAt                                                                >  2:143958/1‑76 (MQ=255)       cGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAAtt                                                               >  2:1315225/1‑76 (MQ=255)       cGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAAtt                                                               >  2:1322346/1‑76 (MQ=255)        gTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTg                                                              <  1:88847/76‑1 (MQ=255)        gTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTg                                                              >  2:1259329/1‑76 (MQ=255)        gTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTg                                                              <  1:158103/76‑1 (MQ=255)         tCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGa                                                             >  1:1666341/1‑76 (MQ=255)          ccGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGat                                                            >  1:18343/1‑76 (MQ=255)            gTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGatat                                                          <  1:1232040/76‑1 (MQ=255)             ttAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATg                                                         >  1:579791/1‑76 (MQ=255)             ttAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATg                                                         <  2:1963682/76‑1 (MQ=255)              tAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGa                                                        >  1:940064/1‑76 (MQ=255)                 cGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGa                                                     <  1:1276033/76‑1 (MQ=255)                    cAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTa                                                  <  1:1897837/76‑1 (MQ=255)                    cAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTa                                                  >  2:1272118/1‑76 (MQ=255)                    cAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTa                                                  <  1:558601/76‑1 (MQ=255)                      ggATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAgg                                                <  1:864485/76‑1 (MQ=255)                       gATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAggg                                               <  1:615409/76‑1 (MQ=255)                               ccATCAGAATAGCGGTAATCACCATCTTAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTaga                                       >  2:899089/1‑76 (MQ=255)                                cATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTagag                                      <  1:535429/76‑1 (MQ=255)                                  tCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGaa                                    <  1:335740/76‑1 (MQ=255)                                  tCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGaa                                    >  2:219927/1‑76 (MQ=255)                                        tAGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgcc                              <  1:612792/76‑1 (MQ=255)                                        gaGCGGTAATCGCCATCTCATGTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgcc                              <  2:1231801/75‑1 (MQ=255)                                         aGCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgcca                             >  1:1520606/1‑76 (MQ=255)                                          gCGGTAATCACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgccag                            <  1:509104/76‑1 (MQ=255)                                                  cACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCa                    >  1:333663/1‑76 (MQ=255)                                                  cACCATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGAAAGTAGAGAAAGTGCCAGCCAGATCa                    >  2:1018035/1‑76 (MQ=255)                                                    ccATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAac                  >  1:1423861/1‑76 (MQ=255)                                                     cATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAaca                 <  1:1024523/76‑1 (MQ=255)                                                     cATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAaca                 <  2:1269267/76‑1 (MQ=255)                                                     cATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAaca                 >  2:1000348/1‑76 (MQ=255)                                                       tctcAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACt               >  1:1149005/1‑76 (MQ=255)                                                        ctcAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCATATCAACACtt              <  1:430481/76‑1 (MQ=255)                                                        ctcAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACtt              >  2:1119450/1‑76 (MQ=255)                                                            aCTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCaa          <  1:1357920/76‑1 (MQ=255)                                                               gCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGc       >  1:337807/1‑76 (MQ=255)                                                                cGAAAGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCt      <  1:806308/76‑1 (MQ=255)                                                                  aaaGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGc    >  2:1838125/1‑76 (MQ=255)                                                                   aaGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGcc   >  2:352928/1‑76 (MQ=255)                                                                    aGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCcg  <  1:187474/76‑1 (MQ=255)                                                                    aGCCTACACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCcg  <  1:5179/76‑1 (MQ=255)                                                                         |                                                                       TGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCG  >  REL606/4097062‑4097204

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑