Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	350,174	C→T	A185T (GCT→ACT)	frmB ←	predicted esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	350,174	0	C	T	95.5%	160.4 / ‑1.4	44	A185T (GCT→ACT)	frmB	predicted esterase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (23/19);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGTTGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGCAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACATATTC  >  REL606/350101‑350246                                                                          |                                                                         cgTTGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTaa                                                                        <  1:665186/76‑1 (MQ=255) cgTTGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTaa                                                                        >  1:39265/1‑76 (MQ=255)  gTTGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTaaa                                                                       <  2:156225/76‑1 (MQ=255)  gTTGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTaaa                                                                       <  1:460638/76‑1 (MQ=255)  gTTGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGCaaa                                                                       <  1:1329473/76‑1 (MQ=255)   ttGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAg                                                                      <  1:1363260/76‑1 (MQ=255)      aCCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCt                                                                   <  2:191682/76‑1 (MQ=255)       ccTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCtt                                                                  >  1:11354/1‑76 (MQ=255)       ccTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCtt                                                                  <  2:87176/76‑1 (MQ=255)       ccTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCtt                                                                  <  2:612472/76‑1 (MQ=255)       ccTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCtt                                                                  <  2:1108227/76‑1 (MQ=255)           gtgAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTg                                                              <  2:309712/76‑1 (MQ=255)           gtgAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTg                                                              <  1:1584348/76‑1 (MQ=255)                  aaGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCAc                                                       >  2:729313/1‑76 (MQ=255)                    gACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACgg                                                     <  2:986729/76‑1 (MQ=255)                    gACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACgg                                                     >  1:826009/1‑76 (MQ=255)                       tCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCAc                                                  >  2:742116/1‑76 (MQ=255)                          ccGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTg                                               <  2:1364849/76‑1 (MQ=255)                            gggTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCg                                             >  1:1737359/1‑76 (MQ=255)                                cGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATgg                                         >  1:1458453/1‑76 (MQ=255)                                cGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATgg                                         >  1:85267/1‑76 (MQ=255)                                 gTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATggg                                        >  1:340689/1‑76 (MQ=255)                                 gTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATggg                                        >  1:898395/1‑76 (MQ=255)                                   aaTCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGga                                      >  1:1509351/1‑76 (MQ=255)                                    aTCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGgag                                     <  1:411003/76‑1 (MQ=255)                                     tCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGgaga                                    >  2:1355888/1‑76 (MQ=255)                                     tCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGgaga                                    >  2:963687/1‑76 (MQ=255)                                     tCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGgaga                                    <  1:1425595/76‑1 (MQ=255)                                     tCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGgaga                                    >  2:1751382/1‑76 (MQ=255)                                               cATCTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGa                          >  1:1692327/1‑76 (MQ=255)                                                  cTTTATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGaaaa                       >  1:712230/1‑76 (MQ=255)                                                   tttATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAAc                      <  1:1679689/76‑1 (MQ=255)                                                   tttATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAAc                      <  1:918209/76‑1 (MQ=255)                                                     tATTTTCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGc                    >  2:1044922/1‑76 (MQ=255)                                                         ttCAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCgag                <  1:705830/76‑1 (MQ=255)                                                           cAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGAc              >  1:502684/1‑76 (MQ=255)                                                           cAGCAAGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGAc              >  1:1339759/1‑76 (MQ=255)                                                               aaGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTg          <  2:1242351/76‑1 (MQ=255)                                                               aaGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTg          >  2:435951/1‑76 (MQ=255)                                                               aaGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTg          <  1:651232/76‑1 (MQ=255)                                                               aaGATATGCAGCAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTg          <  2:40999/76‑1 (MQ=255)                                                                aGATATGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGa         <  2:1330655/76‑1 (MQ=255)                                                                  atatGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACa       >  1:149352/1‑76 (MQ=255)                                                                   tatGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACat      >  2:970254/1‑76 (MQ=255)                                                                    atGCAGTAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACata     >  1:273447/1‑76 (MQ=255)                                                                    aggcagtaAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACata     <  1:1033670/74‑1 (MQ=255)                                                                      gcagtaAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACATAtt   <  2:1491424/76‑1 (MQ=255)                                                                       cagtaAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACATATTc  >  1:345793/1‑76 (MQ=255)                                                                          |                                                                         CGTTGACCTTGTGAAATAAGACTCACCGGGTCGTAATCCAACCAGGCATCTTTATTTTCAGCAAGATATGCAGCAAAGGCTTGCTGTCCCCACGGCACTTGCGATGGGGAGACAATGGGCGAAAACGCCGAGACGCTGACATATTC  >  REL606/350101‑350246

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