Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	712,150	A→T	L218Q (CTG→CAG)	ybfF ←	acyl‑CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	712,150	0	A	T	91.5%	186.5 / 8.4	60	L218Q (CTG→CAG)	ybfF	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (27/27);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCTA  >  NC_000913/712054‑712246                                                                                                 |                                                                                                 cGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTa                                                                                               <  1:251330/100‑1 (MQ=255) cGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTa                                                                                               >  1:1432415/1‑100 (MQ=255)    cGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAat                                                                                            >  2:1161217/1‑100 (MQ=255)       aTAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAatcat                                                                                         >  2:1491172/1‑100 (MQ=255)          gcgcGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCAc                                                                                      <  1:539281/100‑1 (MQ=255)            gcgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACgg                                                                                    >  1:1208858/1‑100 (MQ=255)            gcgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACgg                                                                                    >  1:1209394/1‑100 (MQ=255)                aaTACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGACATTGAGCCTGTAAATCATCACGGTACt                                                                                >  1:1418541/1‑100 (MQ=255)                aaTACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACt                                                                                <  2:1505262/100‑1 (MQ=255)                aaTACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACt                                                                                <  1:76681/100‑1 (MQ=255)                    ccGCATCCGGTTTTTCAGCATGGACCTCGTGACCCGCGCGTGCAATCGCATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTg                                                                               <  1:285453/97‑1 (MQ=255)                    ccGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTc                                                                            >  1:858895/1‑100 (MQ=255)                      gCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGc                                                                          <  1:333114/100‑1 (MQ=255)                       cATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCt                                                                         <  1:254333/100‑1 (MQ=255)                       cATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCt                                                                         <  1:649791/100‑1 (MQ=255)                         tCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTaa                                                                       >  2:2310/1‑100 (MQ=255)                          ccGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAAc                                                                      >  2:82024/1‑100 (MQ=255)                             gTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATa                                                                   >  1:1482373/1‑100 (MQ=255)                             gTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATa                                                                   <  1:897139/100‑1 (MQ=255)                              tttttCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATAc                                                                  <  1:1161215/100‑1 (MQ=255)                                 ttCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGga                                                               <  2:1092293/100‑1 (MQ=255)                                  tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                              <  2:1624474/100‑1 (MQ=255)                                     gCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAAt                                                           >  2:1324035/1‑100 (MQ=255)                                       aTGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTg                                                         >  1:302617/1‑100 (MQ=255)                                       aTGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTg                                                         <  1:735127/100‑1 (MQ=255)                                       aTGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTg                                                         <  2:1443684/100‑1 (MQ=255)                                          gACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgccg                                                      >  1:957217/1‑100 (MQ=255)                                               aGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAgg                                                 >  2:40904/1‑100 (MQ=255)                                               aGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCCCGGTACTGCTCGCTAACATACGAAGAATTGCCGCCAgg                                                 <  2:1259048/100‑1 (MQ=255)                                                  gACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGAt                                              >  2:1066485/1‑100 (MQ=255)                                                      cgcgCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAc                                          >  2:229368/1‑100 (MQ=255)                                                         gcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagg                                       <  1:677481/100‑1 (MQ=255)                                                         gcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagg                                       <  2:605439/100‑1 (MQ=255)                                                          ccTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcaggg                                      >  1:1136112/1‑100 (MQ=255)                                                             gCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcag                                   >  2:749388/1‑100 (MQ=255)                                                              cAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcagg                                  >  2:171262/1‑100 (MQ=255)                                                               aaTCCCATGCGCCCGTGCCTGTGGAAATTGCGCCTGTAACTCATCCCGGTACTGCTCGCTCACATCACGAGAATTCCCGCTAGCCATAAAc                                          >  1:1794911/1‑91 (MQ=255)                                                                aTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGt                                <  1:156272/100‑1 (MQ=255)                                                                 tCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTg                               >  1:1882425/1‑100 (MQ=255)                                                                 tCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTg                               <  1:923463/100‑1 (MQ=255)                                                                  cacaTTCGCCCGTTCCTGTGGAAATTGAGCCCGTAAATCATCCCGGTTCTCCACGCTAACATACGGAGACTTGCCGCCAGGCACAAAcacggcatgg                                 >  1:560656/1‑97 (MQ=255)                                                                     aTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATcc                           <  2:878913/100‑1 (MQ=255)                                                                      tGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATccc                          <  1:1693677/100‑1 (MQ=255)                                                                       gcgcCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTg                               <  1:1511371/94‑1 (MQ=255)                                                                              tGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCggg                  >  1:828967/1‑100 (MQ=255)                                                                              tGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCggg                  >  2:1502770/1‑100 (MQ=255)                                                                                  tgtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAttt              <  2:1482381/100‑1 (MQ=255)                                                                                   gtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAtttt             >  1:360690/1‑100 (MQ=255)                                                                                   gtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAtttt             >  2:1141495/1‑100 (MQ=255)                                                                                    tgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtc            <  1:595963/100‑1 (MQ=255)                                                                                     ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGTATAAACAGGGCAGGGTGATCCCATGCCGGGGTTTtct           >  1:1340785/1‑100 (MQ=255)                                                                                     ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct           <  2:1491623/100‑1 (MQ=255)                                                                                     cgAAATAGAGCCTGTAAATCCTCATGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGGTAAACAGGGTAGGGTGATCCCATGCCGGGAtttt             >  1:1330864/2‑98 (MQ=255)                                                                                      gAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtctc          >  2:1680503/1‑100 (MQ=255)                                                                                      gAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtctc          <  1:1549395/100‑1 (MQ=255)                                                                                       aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc         <  1:1018914/100‑1 (MQ=255)                                                                                        aaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTccc        >  2:345855/1‑100 (MQ=255)                                                                                         aTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCa       <  2:1746582/100‑1 (MQ=255)                                                                                            gAGCCTGTAAATCAGCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcc    <  1:1324028/100‑1 (MQ=255)                                                                                              gCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCTa  >  1:106622/1‑100 (MQ=255)                                                                                                 |                                                                                                 CGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCTA  >  NC_000913/712054‑712246

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