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%	179.2 / NA	56	*257* (TAG→TAA)	bioH	pimeloyl‑ACP methyl ester carboxylesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/27);  total (29/27)

TGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGCA  >  NC_000913/3543945‑3544201                                                                                                                                  |                                                                                                                                tGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCctc                                                                                                                        >  2:255603/1‑139 (MQ=255)             aaaTAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGcc                                                                                                            >  2:1633071/1‑139 (MQ=255)                   tGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagc                                                                                                      >  2:895914/1‑139 (MQ=255)                   tGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGACATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagc                                                                                                      >  1:1083812/1‑139 (MQ=255)                    gCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagca                                                                                                     <  1:544689/139‑1 (MQ=255)                        gAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                                <  2:677213/124‑1 (MQ=255)                        gAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                                >  1:677213/1‑124 (MQ=255)                                                atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaat                                                                         >  1:864693/1‑139 (MQ=255)                                                     gcgcTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                             >  2:610530/1‑82 (MQ=255)                                                     gcgcTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                             <  1:610530/82‑1 (MQ=255)                                                        cTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGgcg                                                                 >  2:423872/1‑139 (MQ=255)                                                          ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCa                                                               >  1:1036213/1‑139 (MQ=255)                                                          ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCa                                                               >  2:1655002/1‑139 (MQ=255)                                                                     tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                                <  1:24417/79‑1 (MQ=255)                                                                     tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                                >  2:24417/1‑79 (MQ=255)                                                                     tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCttt                                                    >  2:142793/1‑139 (MQ=255)                                                                     tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      >  1:472808/1‑137 (MQ=255)                                                                     tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      <  2:472808/137‑1 (MQ=255)                                                                        gaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGc                                                 <  1:423872/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:86278/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  1:1229227/139‑1 (MQ=255)                                                                         aagaagaGTTGACTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:1083812/139‑1 (MQ=255)                                                                            aagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAg                                             <  1:1655002/139‑1 (MQ=255)                                                                            aagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACACCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      >  1:1332503/1‑130 (MQ=255)                                                                            aagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACACCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      <  2:1332503/130‑1 (MQ=255)                                                                                   ggatGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataa                                                                       <  2:1009050/103‑1 (MQ=255)                                                                                   ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAAc                                                                                          <  2:1657175/87‑1 (MQ=255)                                                                                   ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAAc                                                                                          >  1:1657175/1‑87 (MQ=255)                                                                                   ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCg                                                                                   <  2:362322/94‑1 (MQ=255)                                                                                   ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCg                                                                                   >  1:362322/1‑94 (MQ=255)                                                                                   ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataa                                                                       >  1:1009050/1‑106 (MQ=255)                                                                                      tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAcc                                                                                                                           <  2:83413/51‑1 (MQ=255)                                                                                      tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAcc                                                                                                                           >  1:83413/1‑51 (MQ=255)                                                                                      tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCct                                                                                                                         <  2:849543/53‑1 (MQ=255)                                                                                      tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCct                                                                                                                         >  1:849543/1‑53 (MQ=255)                                                                                      tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGAtt                                   >  2:197761/1‑139 (MQ=255)                                                                                         ggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACGCCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATggg                                                           <  1:1367553/112‑1 (MQ=255)                                                                                         ggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATggg                                                           >  2:1367553/1‑112 (MQ=255)                                                                                                        ttCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                <  2:1531399/76‑1 (MQ=255)                                                                                                        ttCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                >  1:1531399/1‑76 (MQ=255)                                                                                                          cTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                <  1:1495485/74‑1 (MQ=255)                                                                                                          cTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                >  2:1495485/1‑74 (MQ=255)                                                                                                                cGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGgcgc                                                                <  2:674791/84‑1 (MQ=255)                                                                                                                cGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGgcgc                                                                >  1:674791/1‑84 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGAc                                                                                               <  2:1587666/51‑1 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGAc                                                                                               >  1:1587666/1‑51 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGgc                                                                  >  1:843717/1‑80 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGgc                                                                  <  2:843717/80‑1 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                >  1:1384264/1‑98 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:1384264/98‑1 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCAt       <  2:864693/139‑1 (MQ=255)                                                                                                                     tttCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAg                                             >  1:717839/1‑98 (MQ=255)                                                                                                                     tttCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAg                                             <  2:717839/98‑1 (MQ=255)                                                                                                                       tCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGcc                       <  2:209145/118‑1 (MQ=255)                                                                                                                       tCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGcc                       >  1:209145/1‑118 (MQ=255)                                                                                                                       tCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGca  <  2:1036213/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                TGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGCA  >  NC_000913/3543945‑3544201

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