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

AGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGCACCACTTTG  >  NC_000913/3543940‑3544209                                                                                                                                       |                                                                                                                                        aGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                                     <  2:279672/139‑1 (MQ=255)   cccTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAcc                                                                                                                                   <  2:372043/139‑1 (MQ=255)             aGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCa                                                                                                                         >  1:548385/1‑139 (MQ=255)                        tGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCct                                                                                                                                 >  2:149063/1‑120 (MQ=255)                        tGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCct                                                                                                                                 <  1:149063/120‑1 (MQ=255)                             gAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagca                                                                                                             <  2:346358/135‑1 (MQ=255)                             gAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACcagca                                                                                                             >  1:346358/1‑135 (MQ=255)                                         ccTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGc                                                                                             >  2:1274214/1‑139 (MQ=255)                                          ctctTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAcc                                                                                                                                   >  2:1429922/1‑100 (MQ=255)                                          ctctTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAcc                                                                                                                                   <  1:1429922/100‑1 (MQ=255)                                               gCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                        <  1:693045/138‑1 (MQ=255)                                               gCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                        >  2:693045/1‑138 (MQ=255)                                                 aGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCg                                                                                     <  1:154776/139‑1 (MQ=255)                                                       catGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACTCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataa                                                                               >  1:947331/1‑139 (MQ=255)                                                              tttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGgcgc                                                                        <  2:947331/139‑1 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                                     <  2:1028455/77‑1 (MQ=255)                                                               ttACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                                     >  1:1028455/1‑77 (MQ=255)                                                                tACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGcc                                                                                            >  2:1336671/1‑117 (MQ=255)                                                                tACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGcc                                                                                            <  1:1336671/117‑1 (MQ=255)                                                                 aCCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATg                                                                     >  2:518355/1‑139 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                                     >  1:691626/1‑66 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                                     <  2:691626/66‑1 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                                        <  2:1037686/79‑1 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCaa                                                                                                                        >  1:1037686/1‑79 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                        <  1:553665/111‑1 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                        >  2:553665/1‑111 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaa                                                                                   >  1:1311074/1‑116 (MQ=255)                                                                          tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaa                                                                                   <  2:1311074/116‑1 (MQ=255)                                                                            agaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACa                                                                                                      <  1:870213/95‑1 (MQ=255)                                                                            agaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACa                                                                                                      >  2:870213/1‑95 (MQ=255)                                                                             gaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGc                                                         <  1:819376/139‑1 (MQ=255)                                                                             gaagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGc                                                         <  1:438467/139‑1 (MQ=255)                                                                              aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                        <  1:1087295/139‑1 (MQ=255)                                                                               agatgggTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAACCTCTGCTTCAAAGCCACCAGCAGGTGACAAACCTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTGGGCGa                                                       <  2:743018/133‑1 (MQ=255)                                                                               agaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGa                                                       <  1:518355/139‑1 (MQ=255)                                                                               agaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGa                                                       <  2:548385/139‑1 (MQ=255)                                                                                      ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                        <  1:1274296/131‑1 (MQ=255)                                                                                      ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                        >  2:1274296/1‑131 (MQ=255)                                                                                        ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                        <  2:994145/129‑1 (MQ=255)                                                                                        ggCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                        >  1:994145/1‑129 (MQ=255)                                                                                           tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                        <  1:767129/126‑1 (MQ=255)                                                                                           tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                        >  2:767129/1‑126 (MQ=255)                                                                                           tGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGAtt                                           >  2:1192044/1‑139 (MQ=255)                                                                                             cggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAccc                                                                                                                                  <  2:669630/50‑1 (MQ=255)                                                                                             cggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACAccc                                                                                                                                  >  1:669630/1‑50 (MQ=255)                                                                                              ggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaa                                                                                   >  2:102213/1‑96 (MQ=255)                                                                                              ggcgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaa                                                                                   <  1:102213/96‑1 (MQ=255)                                                                                                       gTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGcc                               <  1:1192044/139‑1 (MQ=255)                                                                                                            gTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGt                          >  1:377689/1‑139 (MQ=255)                                                                                                                                    cACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGCACCACTTTg  <  1:1274214/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        AGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGCACCACTTTG  >  NC_000913/3543940‑3544209

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