Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP028101	336,797	1 bp→TTT	noncoding (27/1558 nt)	C7Y58_01750 ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP028101	336,797	0	A	T	100.0%	48.3 / NA	26	noncoding (27/1558 nt)	C7Y58_01750	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/8);  total (18/8)
*	CP028101	336,797	1	.	T	92.3%	41.0 / 0.1	26	noncoding (27/1558 nt)	C7Y58_01750	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (18/6);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 8.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01
*	CP028101	336,797	2	.	T	92.3%	42.4 / 0.1	26	noncoding (27/1558 nt)	C7Y58_01750	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (18/6);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 8.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

CAAATTCAAGTAGACTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTA‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACTTATTAT  >  CP028101/336714‑337049                                                                                    |||||                                                                                                                                                                                                                                                           cAAATTCAAGTAGACTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAatttttttattt                                                                                         <  1:240454/251‑1 (MQ=2)       cAAGTAGACTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTc                                                                                   <  1:537086/251‑1 (MQ=9)             gACTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCa                                                                             >  2:845133/1‑251 (MQ=9)             gACTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCa                                                                             >  1:337216/1‑251 (MQ=9)              aCTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTGTGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCat                                                                            >  2:568123/1‑251 (MQ=9)               cTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATTAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCatt                                                                           >  2:259627/1‑251 (MQ=9)                 tGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCatttt                                                                         >  2:122931/1‑250 (MQ=9)                      gtgtTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTAttttt                                                                    >  1:303473/1‑251 (MQ=9)                        gtTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGTTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTGATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTa                                                                   >  1:815789/1‑250 (MQ=255)                        gtTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATGTTTTTTTATTTTTTTTCAATATTTTTATTTTTCTTCTTACCATTTATTTTTa                                                                   >  1:619243/1‑250 (MQ=11)                          tAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACtt                                                                >  2:171744/1‑251 (MQ=9)                                ttCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCATTTTTTTTATTTTTTTTCTTACCatttatt                                                                       >  2:32246/1‑238 (MQ=12)                                ttCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCatttatt                                                                       <  1:32246/238‑1 (MQ=12)                                ttCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTAttt                                                          <  2:635266/251‑1 (MQ=9)                                    gTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTcaa                                                      >  2:445653/1‑251 (MQ=9)                                        gCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACaaa                                                  >  2:244789/1‑251 (MQ=9)                                         cGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACaaaa                                                 >  2:451459/1‑251 (MQ=9)                                           ttCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACaaaaaa                                               >  1:513200/1‑251 (MQ=9)                                               tCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAGCAATAGTTTTTTATTTTTTTTAAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAAt                                             >  1:708688/1‑249 (MQ=11)                                               tCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAAt                                             >  1:908524/1‑249 (MQ=12)                                                   gAGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGAc                                        >  2:556454/1‑250 (MQ=255)                                                    aGCCAGGATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCGCATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGAcc                                       <  1:298310/250‑1 (MQ=255)                                                          gATCAAACTCTTCGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAAc                                  <  1:212971/249‑1 (MQ=255)                                                                      cGTTCAATCTTTTTTTTATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTGACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTGTTATTTTTTTTCTTACCATTTGTTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCtt                    >  1:649663/1‑251 (MQ=14)                                                                        caatctttttttt‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTtatat                <  1:126792/244‑1 (MQ=255)                                                                        caatctttttttt‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTtatat                <  2:418048/244‑1 (MQ=255)                                                                        caatctttttttt‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTtatat                <  1:445653/244‑1 (MQ=255)                                                                        caatctttttttt‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATCGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTtatat                <  1:644852/244‑1 (MQ=255)                                                                               tttttt‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAAc         <  2:260179/251‑1 (MQ=255)                                                                               tttttt‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAAc         <  1:508902/251‑1 (MQ=255)                                                                                ttttT‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACt        <  2:439421/251‑1 (MQ=255)                                                                                ttttT‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTT‑ATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACtt       <  1:82495/251‑1 (MQ=255)                                                                                        tAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACTtattat  <  1:122931/251‑1 (MQ=255)                                                                                        tAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATACCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACTtattat  <  2:619243/251‑1 (MQ=255)                                                                                        tAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATGCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACTtattat  <  2:554165/251‑1 (MQ=255)                                                                                        tAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTCTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACTtattat  <  1:799781/251‑1 (MQ=255)                                                                                        tAGCTCATCTTTGCTAGCTTATTTAAAAACAAATTTTATGGGGGCGTTTTTGTAGGGGAGCTATTCTGTGGCGAATGTCCTTGTCACATTGACATAAAGGAGATTAACAGGAGGAATAGTTTTGGGAAAAAATATTTGGGGAGGTTTTTTCAAGTTGGTTATTTTTTTGAGGACCATTTATTTTTACTTTTATGTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACttaggat  <  1:270539/251‑5 (MQ=255)                                                                                    |||||                                                                                                                                                                                                                                                           CAAATTCAAGTAGACTTGCATGTGTTAAGCATTCTGTCAGCGTTCATCCTGAGCCAGGATCAAACTCTTCGTTCAATCTTTTTA‑‑ATAGCTCATCTTTGCTATCTTATTTAACACCAAATTTTATGGTTGCTTTTTTGTCTTTTCTCTATTCTGTTGCTAATGTCCTTGTCACATTGACATCAATTATATTAACATTATTCCTAGTTTTTGTCAACAATATTTTTTTATTTTTTTTCAATTTTTTTATTTTTTTTCTTACCATTTATTTTTACTTTTATTTTCAACAAAAAAATCTGACCTTAACAGCCAGATTTTCTTATATAATTAACTTATTAT  >  CP028101/336714‑337049

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