Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,544,074	C→T	*257* (TAG→TAA)	bioH ←	pimeloyl‑ACP methyl ester carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,544,074	0	C	T	100.0%	164.4 / NA	52	*257* (TAG→TAA)	bioH	pimeloyl‑ACP methyl ester carboxylesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/25);  total (27/25)

AAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTT  >  NC_000913/3543939‑3544186                                                                                                                                        |                                                                                                                 aaGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCGGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTac                                                                                                               >  2:551846/1‑139 (MQ=255)       tGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCctc                                                                                                         <  2:1300795/139‑1 (MQ=255)                              gAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTg                                                                                  >  1:1046520/1‑139 (MQ=255)                                              ttGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGa                                                                  <  2:539030/139‑1 (MQ=255)                                                 cAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGc                                                               <  1:398447/139‑1 (MQ=255)                                                 cAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGc                                                               >  2:420186/1‑139 (MQ=255)                                                      atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTc                                                                                                   <  1:864456/98‑1 (MQ=255)                                                      atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTc                                                                                                   >  2:864456/1‑98 (MQ=255)                                                      atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCa                                                                                                  <  2:8315/99‑1 (MQ=255)                                                      atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCa                                                                                                  >  1:8315/1‑99 (MQ=255)                                                          tGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataaat                                                      >  2:290180/1‑139 (MQ=255)                                                          tGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataaat                                                      >  2:933407/1‑139 (MQ=255)                                                                ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCa                                                >  1:1331867/1‑139 (MQ=255)                                                                 tACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCAt                                               >  2:1056116/1‑139 (MQ=255)                                                                  aCCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataaat                                                      <  1:300792/131‑1 (MQ=255)                                                                  aCCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataaat                                                      >  2:300792/1‑131 (MQ=255)                                                                  aCCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATg                                              >  1:726492/1‑139 (MQ=255)                                                                    cGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATggg                                            <  2:990498/139‑1 (MQ=255)                                                                           tCAGGAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaata                                                         >  1:855012/1‑119 (MQ=255)                                                                           tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                               >  2:1154374/1‑81 (MQ=255)                                                                           tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                               <  1:1154374/81‑1 (MQ=255)                                                                           tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGa                                                                                 <  1:371791/95‑1 (MQ=255)                                                                           tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGa                                                                                 >  2:371791/1‑95 (MQ=255)                                                                           tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaata                                                         <  2:855012/119‑1 (MQ=255)                                                                           tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCttt                                     >  2:721060/1‑139 (MQ=255)                                                                             agaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTgg                                   <  2:482398/139‑1 (MQ=255)                                                                             agaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTgg                                   <  1:290180/139‑1 (MQ=255)                                                                               aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                 <  1:420186/139‑1 (MQ=255)                                                                                 gaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGaa                               <  1:1056116/139‑1 (MQ=255)                                                                                     agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGc                                  >  2:409055/1‑132 (MQ=255)                                                                                     agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGc                                  <  1:409055/132‑1 (MQ=255)                                                                                         ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGa                                >  1:354522/1‑130 (MQ=255)                                                                                         ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGa                                <  2:354522/130‑1 (MQ=255)                                                                                            tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGAtt                    >  2:644870/1‑139 (MQ=255)                                                                                              cggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                 >  2:354665/1‑92 (MQ=255)                                                                                              cggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                 <  1:354665/92‑1 (MQ=255)                                                                                              cggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCg                  >  2:182463/1‑139 (MQ=255)                                                                                                  gAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCtgt              <  1:1164606/139‑1 (MQ=255)                                                                                                      cAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAgg          <  1:182463/139‑1 (MQ=255)                                                                                                              ttCTGTCTCGCCATTTTAATCGCCACTCACACCCTCCGCTTCAACGCCACCAGCAGGTGACAAATCTTGACCGGATGCGAAATAAATGGTGAATGGGCCGCTTTTGCGAAGATATATGTTTCCCTGTGAGGCCAAAGtt  >  1:318804/1‑139 (MQ=255)                                                                                                               tCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagca                                                                                      <  1:1305115/54‑1 (MQ=255)                                                                                                               tCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagca                                                                                      >  2:1305115/1‑54 (MQ=255)                                                                                                                 tGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACt                                                                          <  2:1269487/64‑1 (MQ=255)                                                                                                                 tGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACt                                                                          >  1:1269487/1‑64 (MQ=255)                                                                                                                   tctcGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                               >  2:461170/1‑41 (MQ=255)                                                                                                                   tctcGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                               <  1:461170/41‑1 (MQ=255)                                                                                                                           tttCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAACGGCGCATGGGCCGCTTTGGCg                                 >  2:453906/1‑95 (MQ=255)                                                                                                                           tttCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAACGGCGCATGGGCCGCTTTGGCg                                 <  1:453906/95‑1 (MQ=255)                                                                                                                                aaaGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                 <  2:552834/58‑1 (MQ=255)                                                                                                                                aaaGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                 >  1:552834/1‑58 (MQ=255)                                                                                                                                   gCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaat                                                          >  2:599103/1‑62 (MQ=255)                                                                                                                                   gCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaat                                                          <  1:599103/62‑1 (MQ=255)                                                                                                                                        |                                                                                                                 AAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTT  >  NC_000913/3543939‑3544186

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