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

AGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCA  >  NC_000913/3543940‑3544195                                                                                                                                       |                                                                                                                          aGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       <  1:1721389/139‑1 (MQ=255)      tGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCctc                                                                                                                  >  2:954871/1‑139 (MQ=255)             aGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCa                                                                                                           >  2:1436602/1‑139 (MQ=255)                  aaaTAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTa                                                                                                                         >  2:64172/1‑120 (MQ=255)                  aaaTAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGcc                                                                                                      <  1:954871/139‑1 (MQ=255)                  aaaTAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGcc                                                                                                      >  2:913872/1‑139 (MQ=255)                        tGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagc                                                                                                >  2:368948/1‑139 (MQ=255)                         gCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagca                                                                                               <  1:368948/139‑1 (MQ=255)                             gAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTg                                                                                           >  1:54165/1‑139 (MQ=255)                                tCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACa                                                                                        <  2:203195/139‑1 (MQ=255)                                          ctctTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGcc                                                                              >  1:417354/1‑139 (MQ=255)                                                     atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaat                                                                   >  2:1797508/1‑139 (MQ=255)                                                        atGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataaa                                                                >  1:1472263/1‑139 (MQ=255)                                                         tGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataaat                                                               <  2:417354/139‑1 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                          <  2:41646/90‑1 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                          >  1:41646/1‑90 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcag                                                                                                 >  1:50762/1‑99 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcag                                                                                                 <  2:50762/99‑1 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                          <  2:137268/122‑1 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                          >  1:137268/1‑122 (MQ=255)                                                                 aCCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCg                                                                       >  2:70412/1‑123 (MQ=255)                                                                 aCCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCGCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCg                                                                       <  1:70412/123‑1 (MQ=255)                                                                    gggCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAAAGGCGCATGGGc                                                    <  2:54165/139‑1 (MQ=255)                                                                       cTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGc                                                 >  1:988649/1‑139 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       >  1:874014/1‑66 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       <  2:874014/66‑1 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCttt                                              >  1:1607977/1‑139 (MQ=255)                                                                            agaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       >  2:1167168/1‑64 (MQ=255)                                                                            agaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       <  1:1167168/64‑1 (MQ=255)                                                                              aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                          <  2:1201793/139‑1 (MQ=255)                                                                              aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                          <  2:1464293/139‑1 (MQ=255)                                                                              aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                          <  2:988649/139‑1 (MQ=255)                                                                              aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                          <  1:1492806/139‑1 (MQ=255)                                                                              aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                          <  1:1254128/139‑1 (MQ=255)                                                                               agaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       >  1:1514515/1‑61 (MQ=255)                                                                               agaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       <  2:1514515/61‑1 (MQ=255)                                                                               agaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGa                                         <  1:913872/139‑1 (MQ=255)                                                                                    agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                                        <  2:1830767/71‑1 (MQ=255)                                                                                    agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                                        >  1:1830767/1‑71 (MQ=255)                                                                                    agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                >  1:333287/1‑127 (MQ=255)                                                                                    agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                <  2:333287/127‑1 (MQ=255)                                                                                    agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGata                                    >  1:1118194/1‑139 (MQ=255)                                                                                      ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcag                                                                                                 >  1:723171/1‑76 (MQ=255)                                                                                      ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcag                                                                                                 <  2:723171/76‑1 (MQ=255)                                                                                      ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                          <  1:454303/131‑1 (MQ=255)                                                                                      ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                          >  2:454303/1‑131 (MQ=255)                                                                                        ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       <  2:737779/52‑1 (MQ=255)                                                                                        ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                       >  1:737779/1‑52 (MQ=255)                                                                                        ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAATGCCACCAGCAGGTGACAAAACTCGGcc                                                                              >  1:492163/1‑93 (MQ=255)                                                                                        ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAATGCCACCAGCAGGTGACAAAACTCGGcc                                                                              <  2:492163/93‑1 (MQ=255)                                                                                                 gAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCtgt                       <  1:1797508/139‑1 (MQ=255)                                                                                                   gTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTc                                                                                                            >  1:1288835/1‑52 (MQ=255)                                                                                                   gTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTc                                                                                                            <  2:1288835/52‑1 (MQ=255)                                                                                                       gTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACCCCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGcc                 <  2:1472263/139‑1 (MQ=255)                                                                                                       gTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGcc                 <  1:1436602/139‑1 (MQ=255)                                                                                                       gTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGcc                 <  2:1515339/139‑1 (MQ=255)                                                                                                          aaGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCaaa              <  1:1126476/139‑1 (MQ=255)                                                                                                                tGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                                        >  1:1386716/1‑43 (MQ=255)                                                                                                                tGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACg                                                                                                        <  2:1386716/43‑1 (MQ=255)                                                                                                                      gCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCa  >  2:81645/1‑139 (MQ=255)                                                                                                                                  gCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGcc                                                                              >  1:1821399/1‑51 (MQ=255)                                                                                                                                  gCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGcc                                                                              <  2:1821399/51‑1 (MQ=255)                                                                                                                                       |                                                                                                                          AGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCA  >  NC_000913/3543940‑3544195

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