Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,200,626	A→G	M21V (ATG→GTG)	yjaH →	DUF1481 domain‑containing protein YjaH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,200,626	0	A	G	100.0%	204.9 / NA	56	M21V (ATG→GTG)	yjaH	DUF1481 domain‑containing protein YjaH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/31);  total (25/31)

CGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGATAATGTGCATCTGC  >  NC_000913/4200496‑4200761                                                                                                                                   |                                                                                                                                        cGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTgctggctg                                                                                                                                <  1:376268/140‑1 (MQ=255)    aGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGgcg                                                                                                                             <  1:624460/140‑1 (MQ=255)    aGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGgcg                                                                                                                             <  1:453752/140‑1 (MQ=255)    aGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGgcg                                                                                                                             <  2:657888/140‑1 (MQ=255)     gCTAACGTACCCGCATTCTTTGTTGTCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTTGCAGGGAACAGTTTTAATGAATGGGTGGTTTCGCCCCTTGTGTCTTTCTGGCGTCGTTCATTGGTgctagctttagt                                                                                                                            >  2:582518/1‑135 (MQ=255)        aaCGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGctct                                                                                                                         <  2:705969/140‑1 (MQ=255)               cGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTtctc                                                                                                                  >  1:107156/1‑140 (MQ=255)               cGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTtctc                                                                                                                  <  1:979759/140‑1 (MQ=255)               cGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTtctc                                                                                                                  <  1:209107/140‑1 (MQ=255)                ggCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCa                                                                                                                 <  2:1001925/140‑1 (MQ=255)                ggCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCa                                                                                                                 >  2:227373/1‑140 (MQ=255)                 gCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCAc                                                                                                                >  1:207133/1‑140 (MQ=255)                 gCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCAc                                                                                                                <  1:745437/140‑1 (MQ=255)                 gCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCAc                                                                                                                <  2:214544/140‑1 (MQ=255)                   aTTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTg                                                                                                              >  1:177773/1‑140 (MQ=255)                      tGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCt                                                                                                           <  1:493674/140‑1 (MQ=255)                      tGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCt                                                                                                           <  1:890442/140‑1 (MQ=255)                       gTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTg                                                                                                          <  2:268231/140‑1 (MQ=255)                         ttCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAATTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTa                                                                                                        <  2:496817/140‑1 (MQ=255)                            tGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTc                                                                                                     >  1:807305/1‑140 (MQ=255)                              gCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCAt                                                                                                   <  2:946713/140‑1 (MQ=255)                                aaGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATaa                                                                                                 <  2:153064/140‑1 (MQ=255)                                aaGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATaa                                                                                                 <  2:107156/140‑1 (MQ=255)                                     aCTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCtt                                                                                            <  2:656939/140‑1 (MQ=255)                                         aaaGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACt                                                                                        >  1:701763/1‑140 (MQ=255)                                          aaGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACtt                                                                                       <  1:147298/140‑1 (MQ=255)                                             aCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTtcct                                                                                    <  2:641753/140‑1 (MQ=255)                                              cGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTtcctc                                                                                   <  1:330748/140‑1 (MQ=255)                                              cGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTtcctc                                                                                   <  1:50315/140‑1 (MQ=255)                                                     aagtaagATTGCGTGGCAGTGAACAGTTTTAACGAGGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTAc                                                                            <  2:906905/140‑1 (MQ=255)                                                        taagATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGc                                                                         >  2:12744/1‑140 (MQ=255)                                                             ttGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTg                                                                    >  1:429908/1‑140 (MQ=255)                                                             ttGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTg                                                                    >  1:457587/1‑140 (MQ=255)                                                               gCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGa                                                                  >  1:274781/1‑140 (MQ=255)                                                               gCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGa                                                                  >  2:867522/1‑140 (MQ=255)                                                                 gTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAtt                                                                <  2:98753/140‑1 (MQ=255)                                                                      aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg                                                           >  2:686266/1‑140 (MQ=255)                                                                      aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg                                                           >  2:353923/1‑140 (MQ=255)                                                                           aCAGTTTTAACGAAGGGGTGGTTTCCCCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc                                                      <  2:800780/140‑1 (MQ=255)                                                                             aGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCa                                                    <  1:534533/140‑1 (MQ=255)                                                                               ttttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAgg                                                  >  2:28004/1‑140 (MQ=255)                                                                                 ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGgc                                                <  1:803545/140‑1 (MQ=255)                                                                                  tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCTTTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGgcg                                               <  2:807305/140‑1 (MQ=255)                                                                                  tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGgcg                                               >  2:604799/1‑140 (MQ=255)                                                                                    aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGcgcg                                             >  2:370053/1‑140 (MQ=255)                                                                                    aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGcgcg                                             >  2:609995/1‑140 (MQ=255)                                                                                    aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGcgcg                                             >  2:205572/1‑140 (MQ=255)                                                                                               gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTg                                  >  1:784714/1‑140 (MQ=255)                                                                                                    aCCCCTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGaa                             >  2:535651/1‑140 (MQ=255)                                                                                                      cccTTTTGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGaaaa                           >  1:527274/1‑140 (MQ=255)                                                                                                         ttttGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGAc                        >  1:488023/1‑140 (MQ=255)                                                                                                         ttttGTCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGAc                        >  2:565033/1‑140 (MQ=255)                                                                                                              tCTTTCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCCCTCCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAgcgg                   <  2:693248/140‑1 (MQ=255)                                                                                                                tttCTGGCGTCAATATTTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATACCTCTTCACTTCCCCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAgcggcg                 <  1:495129/140‑1 (MQ=255)                                                                                                                 ttCTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGa                >  2:608199/1‑140 (MQ=255)                                                                                                                   cTGGCGTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTTCCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGATa              <  2:177773/140‑1 (MQ=255)                                                                                                                    gggCGTCGACCATTGGTGCTGGCTGGCGCTCTGCTTCTCACGGCCTGTAGTCATACCTCTACCCTTCCCCCCTTTCCCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGATaa             <  2:701763/139‑1 (MQ=255)                                                                                                                        gTCGATCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGATAAtgtg         <  1:651260/140‑1 (MQ=255)                                                                                                                            aTCATTGGTGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGATAATGTGCATc     >  2:1031819/1‑140 (MQ=255)                                                                                                                               aTTGGTGTTGGTTGGCGCTTCGCTTCTCCTTGCTTTTGGTCTTAACTCTTCTCTTCCTCCATTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGATAATGTGCATCTGc  <  2:207133/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        CGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCTGGCGCTCTGCTTCTCACTGCCTGTAGTCATAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAGGGCGCGGTACGCATCTGGCGAAAAGACAGCGGCGATAATGTGCATCTGC  >  NC_000913/4200496‑4200761

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