Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	92,879	C→T	Q148* (CAG→TAG)	yabB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	92,879	0	C	T	96.6%	215.4 / ‑1.0	59	Q148* (CAG→TAG)	yabB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (31/26);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AACCTGGCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGCAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAAT  >  REL606/92806‑92948                                                                          |                                                                      aacCTGGCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAg                                                                     >  2:945453/1‑76 (MQ=255) aacCTGGCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAg                                                                     >  2:903145/1‑76 (MQ=255) aacCTGGCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAg                                                                     >  2:741711/1‑76 (MQ=255)     tGGCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACt                                                                 >  2:483957/1‑76 (MQ=255)      ggCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACtt                                                                >  1:307846/1‑76 (MQ=255)      ggCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACtt                                                                <  2:1609591/76‑1 (MQ=255)      ggCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGCAGGACtt                                                                >  1:1298245/1‑76 (MQ=255)      ggCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGCAGGACtt                                                                >  1:1298256/1‑76 (MQ=255)       gCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTg                                                               <  2:333903/76‑1 (MQ=255)       gCATCAACAGGGCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTg                                                               <  1:1826025/76‑1 (MQ=255)         aTCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGtc                                                             >  2:1615470/1‑76 (MQ=255)           cAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGtctc                                                           <  1:1693601/76‑1 (MQ=255)             aCAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCta                                                         <  1:1498233/76‑1 (MQ=255)              cAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCtat                                                        >  1:1780167/1‑76 (MQ=255)               aGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCtata                                                       >  2:445965/1‑76 (MQ=255)                ggTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATaa                                                      <  1:1007500/76‑1 (MQ=255)                  tCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATaaaa                                                    >  1:655385/1‑76 (MQ=255)                   cAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAat                                                   >  1:1022134/1‑76 (MQ=255)                   cAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAat                                                   >  1:703790/1‑76 (MQ=255)                    aaGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAatg                                                  >  2:1852769/1‑76 (MQ=255)                      ggAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAatgat                                                <  2:540231/76‑1 (MQ=255)                       gAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAatgatg                                               >  2:1501717/1‑76 (MQ=255)                        aaGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATgg                                              >  1:1814624/1‑76 (MQ=255)                        aaGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATgg                                              >  2:1874788/1‑76 (MQ=255)                          gATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGaa                                            >  1:1161850/1‑76 (MQ=255)                            tcncGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGaaaa                                          >  2:1405519/4‑76 (MQ=255)                            tatCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGaaaa                                          >  1:103054/1‑76 (MQ=255)                            tatCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGaaaa                                          <  2:730158/76‑1 (MQ=255)                              tCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACt                                        <  1:17997/76‑1 (MQ=255)                               cGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACta                                       >  1:1205856/1‑76 (MQ=255)                               cGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGCCTCTATAAAATGATGGAAAACta                                       >  1:33579/1‑76 (MQ=255)                                gACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACtat                                      <  1:808886/76‑1 (MQ=255)                                    cAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAAc                                  <  2:942554/76‑1 (MQ=255)                                    cAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAAc                                  <  2:506253/76‑1 (MQ=255)                                      gagCAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACAt                                <  1:1286080/76‑1 (MQ=255)                                         cAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACAtact                             >  2:20389/1‑76 (MQ=255)                                         cAGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACAtact                             >  1:865629/1‑76 (MQ=255)                                          aGTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACAtacta                            >  2:897459/1‑76 (MQ=255)                                           gTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACAtactac                           <  1:550192/76‑1 (MQ=255)                                           gTTGGCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACAtactac                           <  1:335898/76‑1 (MQ=255)                                              ggCTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGt                        <  2:1580844/76‑1 (MQ=255)                                                cTACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGtgc                      >  1:1786720/1‑76 (MQ=255)                                                 tACCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGtgct                     <  1:133518/76‑1 (MQ=255)                                                  aCCGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGtgctg                    >  2:1874101/1‑76 (MQ=255)                                                    cGGAGACTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTgctgct                  >  1:550914/1‑76 (MQ=255)                                                      gagaCTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTgg                >  1:1419575/1‑76 (MQ=255)                                                        gaCTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGAt              <  1:262500/76‑1 (MQ=255)                                                          cTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGa            >  1:1550701/1‑76 (MQ=255)                                                          cTTATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGa            >  1:1650684/1‑76 (MQ=255)                                                           ttATCGGAGCGACTGTAGGTCTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGan           <  1:1615616/76‑2 (MQ=255)                                                           ttATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGaa           <  1:415161/76‑1 (MQ=255)                                                            tATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAg          <  2:267532/76‑1 (MQ=255)                                                            tATCGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAg          <  1:895930/76‑1 (MQ=255)                                                               cGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCg       >  2:1886186/1‑76 (MQ=255)                                                               cGGAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCg       <  1:1846132/76‑1 (MQ=255)                                                                ggAGCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGt      <  1:1874789/76‑1 (MQ=255)                                                                   gCGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTaa   <  1:915249/76‑1 (MQ=255)                                                                    cGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAAt  <  2:339183/76‑1 (MQ=255)                                                                    cGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAAt  >  1:730855/1‑76 (MQ=255)                                                                    cGACTGTAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAAt  >  1:46297/1‑76 (MQ=255)                                                                          |                                                                      AACCTGGCATCAACAGGTCAAGGAAGATATCGACGCAGAGCAGTTGGCTACCGGAGACTTATCGGAGCGACTGCAGGACTTGTCTCTATAAAATGATGGAAAACTATAAACATACTACGGTGCTGCTGGATGAAGCCGTTAAT  >  REL606/92806‑92948

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