Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,996,372	C→T	*302* (TAG→TAA)	insD1 ←	IS2 transposase TnpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,996,372	0	C	T	100.0%	216.7 / NA	68	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/31);  total (37/31)

AAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTG  >  NC_000913/2996242‑2996496                                                                                                                                   |                                                                                                                             aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTc                                                                                                                      >  1:1192129/1‑139 (MQ=255)      atCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATCGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGAc                                                                                                                 >  1:628653/1‑139 (MQ=255)        cccTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATCGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACAt                                                                                                               <  2:628653/139‑1 (MQ=255)            tAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGt                                                                                                           >  1:184470/1‑139 (MQ=255)                  aTCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCAc                                                                                                     >  1:1419192/1‑139 (MQ=255)                   tCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    >  1:1628784/1‑139 (MQ=255)                    cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         >  2:156318/1‑133 (MQ=255)                    cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         <  1:156318/133‑1 (MQ=255)                    cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                   >  1:799712/1‑139 (MQ=255)                       gAACCTTATGAAATAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTGATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAc                                                                                                >  1:846951/1‑139 (MQ=255)                        aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                               <  1:358798/139‑1 (MQ=255)                         aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAGAAGTGACAAGCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                               >  1:1107538/1‑138 (MQ=255)                         aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTCAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                               <  2:1107538/138‑1 (MQ=255)                         aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                          <  1:935909/127‑1 (MQ=255)                         aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                          >  2:935909/1‑127 (MQ=255)                         aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAccc                                                                                              >  2:1464222/1‑139 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATAttt                                                                                                                       <  1:386787/110‑1 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATAttt                                                                                                                       >  2:386787/1‑110 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         <  2:1223423/124‑1 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         >  1:1223423/1‑124 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         >  2:1335182/1‑124 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         <  1:1335182/124‑1 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTa                                                                                          >  1:27096/1‑139 (MQ=255)                             tataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTa                                                                                          >  1:1090816/1‑139 (MQ=255)                              ataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATTAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    <  1:247164/128‑1 (MQ=255)                              ataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTAc                                                                                         >  2:362782/1‑139 (MQ=255)                              ataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATAACt                                                                                                    >  2:247164/1‑128 (MQ=255)                                     aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    >  1:789616/1‑121 (MQ=255)                                     aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    <  2:789616/121‑1 (MQ=255)                                     aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                  <  2:799712/139‑1 (MQ=255)                                      aCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCa                                                                                                                    <  1:1441101/104‑1 (MQ=255)                                      aCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCa                                                                                                                    >  2:1441101/1‑104 (MQ=255)                                          gcgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                             <  2:1419192/139‑1 (MQ=255)                                                          tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                             <  1:362782/139‑1 (MQ=255)                                                          tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                             >  2:877845/1‑139 (MQ=255)                                                          tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                             >  1:892568/1‑139 (MQ=255)                                                          tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                             >  1:1502344/1‑139 (MQ=255)                                                          tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                             <  2:1090816/139‑1 (MQ=255)                                                            tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  2:1628784/139‑1 (MQ=255)                                                                          aaCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  1:636660/125‑1 (MQ=255)                                                                          aaCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           >  2:636660/1‑125 (MQ=255)                                                                           aCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                            <  1:1464222/139‑1 (MQ=255)                                                                             cACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                          <  2:27096/139‑1 (MQ=255)                                                                               cTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTAc                                                                                         >  2:147138/1‑90 (MQ=255)                                                                               cTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTAc                                                                                         <  1:147138/90‑1 (MQ=255)                                                                               cTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                            <  1:719675/135‑1 (MQ=255)                                                                               cTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                            >  2:719675/1‑135 (MQ=255)                                                                               cTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACt                                        <  1:877845/139‑1 (MQ=255)                                                                                       aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                        >  2:1303406/1‑99 (MQ=255)                                                                                       aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                        <  1:1303406/99‑1 (MQ=255)                                                                                       aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAt                                >  1:207569/1‑139 (MQ=255)                                                                                          gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                             <  2:1192129/139‑1 (MQ=255)                                                                                          gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                             <  2:829398/139‑1 (MQ=255)                                                                                          gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                             <  2:184470/139‑1 (MQ=255)                                                                                          gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                             >  1:301361/1‑139 (MQ=255)                                                                                            tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAccc                                                <  1:329207/118‑1 (MQ=255)                                                                                            tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAccc                                                >  2:329207/1‑118 (MQ=255)                                                                                            tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                           >  2:998658/1‑139 (MQ=255)                                                                                            tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                           >  2:998216/1‑139 (MQ=255)                                                                                              aaTGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                            >  2:1363770/1‑136 (MQ=255)                                                                                              aaTGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                            <  1:1363770/136‑1 (MQ=255)                                                                                                tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGt                       >  1:1044406/1‑139 (MQ=255)                                                                                                  gAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCAATAACCCAGCGCACTATGCGGATGCCATTCGTta                     <  2:1502344/139‑1 (MQ=255)                                                                                                     gTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATaa                  >  1:410463/1‑139 (MQ=255)                                                                                                            ccAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGa           <  2:1044406/139‑1 (MQ=38)                                                                                                             cAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                                    >  2:1087098/1‑97 (MQ=255)                                                                                                             cAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                                    <  1:1087098/97‑1 (MQ=255)                                                                                                                     tATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTg  >  1:180753/1‑139 (MQ=11)                                                                                                                                   |                                                                                                                             AAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTG  >  NC_000913/2996242‑2996496

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