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

TGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGCA  >  NC_000913/3543945‑3544201                                                                                                                                  |                                                                                                                                tGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCctc                                                                                                                        <  2:500280/139‑1 (MQ=255) tGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCctc                                                                                                                        >  1:1273996/1‑139 (MQ=255)      ataGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCtt                                                                                                                   >  1:1652150/1‑139 (MQ=255)        aGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCa                                                                                                                 >  1:1499539/1‑139 (MQ=255)          ccAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAAc                                                                                                               <  1:1306869/139‑1 (MQ=255)          ccAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAAc                                                                                                               >  1:149532/1‑139 (MQ=255)             aaaTAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGcc                                                                                                            >  1:442306/1‑139 (MQ=255)                taatGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGccacc                                                                                                         >  1:388692/1‑139 (MQ=255)                      ttGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAgg                                                                                                   <  2:333217/139‑1 (MQ=255)                        gAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTg                                                                                                 >  1:1523951/1‑139 (MQ=255)                                aGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACt                                                                                         >  2:739614/1‑139 (MQ=255)                                aGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACt                                                                                         >  1:406074/1‑139 (MQ=255)                                    ccTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGc                                                                                     >  1:811097/1‑139 (MQ=255)                                     ctctTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTc                                                                                                                  <  1:1431015/109‑1 (MQ=255)                                     ctctTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTc                                                                                                                  >  2:1431015/1‑109 (MQ=255)                                                atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                             <  1:1322222/87‑1 (MQ=255)                                                atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTaca                                                                                                                             >  2:1322222/1‑87 (MQ=255)                                                atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaat                                                                         >  1:582544/1‑139 (MQ=255)                                                atcatGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaat                                                                         >  1:466800/1‑139 (MQ=255)                                                  catGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGc                                                                              >  1:902139/1‑132 (MQ=255)                                                  catGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGc                                                                              <  2:902139/132‑1 (MQ=255)                                                    tGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGaaataaat                                                                     >  2:547740/1‑139 (MQ=255)                                                     gcgcTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATg                                                                    <  2:1652150/139‑1 (MQ=255)                                                                ggCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGcc                                                         >  2:426684/1‑139 (MQ=255)                                                                ggCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGcc                                                         >  2:374778/1‑139 (MQ=255)                                                                ggCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGcc                                                         >  1:1489822/1‑139 (MQ=255)                                                                   tGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCgg                                                                                  <  1:597265/111‑1 (MQ=255)                                                                   tGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCgg                                                                                  >  2:597265/1‑111 (MQ=255)                                                                   tGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      <  1:547740/139‑1 (MQ=255)                                                                   tGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      <  2:466800/139‑1 (MQ=255)                                                                    gTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCtt                                                     <  2:218440/139‑1 (MQ=255)                                                                     tCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCTCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCttt                                                    >  2:1170374/1‑139 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:1499539/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:1024627/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:1523951/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  1:77590/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  1:701051/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:582544/139‑1 (MQ=255)                                                                         aagaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  1:1262042/139‑1 (MQ=255)                                                                          agaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGa                                               <  2:406074/139‑1 (MQ=255)                                                                          agaagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGa                                               <  2:535460/139‑1 (MQ=255)                                                                            aagaGTTGTCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAg                                             <  1:983141/139‑1 (MQ=255)                                                                            aagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAg                                             <  1:483464/139‑1 (MQ=255)                                                                            aagaGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAg                                             <  2:1489822/139‑1 (MQ=255)                                                                               agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGc                                                                                     >  2:974473/1‑96 (MQ=255)                                                                               agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGc                                                                                     <  1:974473/96‑1 (MQ=255)                                                                               agTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGata                                          >  1:400389/1‑139 (MQ=255)                                                                                 ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACt                                                                                         <  2:943684/90‑1 (MQ=255)                                                                                 ttGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACt                                                                                         >  1:943684/1‑90 (MQ=255)                                                                                           cgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCtg                              >  2:722701/1‑139 (MQ=255)                                                                                           cgAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCtg                              >  1:476380/1‑139 (MQ=255)                                                                                             aGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                >  2:456934/1‑87 (MQ=255)                                                                                             aGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGAt                                                                                <  1:456934/87‑1 (MQ=255)                                                                                                        ttCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGc                                                 <  2:496083/107‑1 (MQ=255)                                                                                                        ttCTGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGc                                                 >  1:496083/1‑107 (MQ=255)                                                                                                           tGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  2:1452793/105‑1 (MQ=255)                                                                                                           tGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                >  1:1452793/1‑105 (MQ=255)                                                                                                           tGTCTCGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTAt              >  1:471076/1‑139 (MQ=255)                                                                                                             tctcGCCATTTCAAAAGCCACTTAGACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATcc            >  2:1122850/1‑139 (MQ=255)                                                                                                             tctcGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                <  1:919012/103‑1 (MQ=255)                                                                                                             tctcGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCg                                                >  2:919012/1‑103 (MQ=255)                                                                                                                cGCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGc         <  2:149532/139‑1 (MQ=255)                                                                                                                 gCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGc                                                                                     <  1:969332/62‑1 (MQ=255)                                                                                                                 gCCATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGc                                                                                     >  2:969332/1‑62 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      <  2:1131955/92‑1 (MQ=255)                                                                                                                  ccATTTCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCt                                                      >  1:1131955/1‑92 (MQ=255)                                                                                                                       tCAAAAGCCACTTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGca  <  1:722701/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                TGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGATGCGAAATAAATGGCGCATGGGCCGCTTTGGCGAAGATATATGATTCGCTGTGAGGCCAAAGTTTATCCAGCATCGGCA  >  NC_000913/3543945‑3544201

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