Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	90,908	0	A	T	30.3%	54.1 / 22.7	34	Q272L (CAG→CTG)	rsmH	16S rRNA m(4)C1402 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (17/6);  new base T (6/4);  total (24/10)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTG  >  NC_000913/90669‑91153                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      ggTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAAc                                                                                                                                                                                                                                            >  1:246321/1‑251 (MQ=255)    gAAAGATAAGTGTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCTGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCTTTAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTTCTCAAAAAACTg                                                                                                                                                                                                                                          <  2:32652/250‑1 (MQ=255)     aaaGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTggg                                                                                                                                                                                                                                        >  1:244310/1‑251 (MQ=255)       aGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGGTAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTg                                                                                                                                                                                                                                      >  1:122857/1‑251 (MQ=255)       aGATAAGTTTAAACATCCCGCGACACGTACCTTCCAGGCGGTGCGCATTTGGGAAAACAGTTAACTGGAGGAAATAGAGCAGGCGCTAAAAAGCTCGCTCAACGAGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTCTTGTGAAACGTTTTATGCGTGAAAACGGCCGCGGTCCGCAAGTTCCGTCAGGGTTACCGATGACTGAAGAGCAGCTCCAAAAACTGGGTg                                                                                                                                                                                                                                      >  2:261863/1‑251 (MQ=255)                   aCATCCCGCGCCCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCACTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGCTGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTgg                                                                                                                                                                                                                          >  1:33787/1‑250 (MQ=255)                      tCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATTACTGAAGAGCTGCTCAAAAAAc                                                                                                                                                                                                                                            <  2:217200/230‑1 (MQ=255)                      tCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGCCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCGGCTCAAAAAAc                                                                                                                                                                                                                                            >  1:217200/1‑230 (MQ=255)                            gACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACt                                                                                                                                                                                                                   >  1:255365/1‑249 (MQ=255)                                 gTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCa                                                                                                                                                                                                              >  1:240960/1‑249 (MQ=255)                                         cAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGcc                                                                                                                                                                                                    >  1:299143/1‑251 (MQ=255)                                                   gcATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCgaagc                                                                                                                                                                                           >  1:147816/1‑249 (MQ=255)                                                                                  agagCAGGCGCTAAAACGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTCTTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTAACGATGACTGACGAGCAGCTCAAAAACCTGGGTGGCCGTAAGCTGCGAGCAATAGGCACGTTAATGCCGGGCGCAGAAGAGTTGGCTGAGAACCCTCGTGACCGTagtt                                                                                                                                                           >  1:125995/1‑251 (MQ=255)                                                                                    agCAGGCGCTAAAAAGCTCGCTAAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTagttc                                                                                                                                                          >  2:71481/1‑250 (MQ=255)                                                                                       aGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTagttcagtt                                                                                                                                                      >  1:77114/1‑251 (MQ=255)                                                                                         gcgcTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCt                                                                                                                                                    <  2:186119/251‑1 (MQ=255)                                                                                            cTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGCAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCAGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTGGGTGGCCGGCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCACGTGCCCGTAGTGCAGTTCTGCg                                                                                                                                                 >  2:262405/1‑251 (MQ=255)                                                                                                 aaGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCTGGCGAAGAAGAGGTGGCTGAGAGCCCTCGTGCCCGTAGTTCAGTTCTGCGTATTg                                                                                                                                            <  2:35338/251‑1 (MQ=255)                                                                                                  aGCTCGCGCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGc                                                                                                                                           <  2:223346/251‑1 (MQ=255)                                                                                                                         gggTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGAACGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATc                                                                                                                      >  1:191630/1‑249 (MQ=255)                                                                                                                             gggCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAA‑CTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATcagcag                                                                                                                 >  2:160835/1‑249 (MQ=255)                                                                                                                                               aGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGc                                                                                              >  1:13690/1‑251 (MQ=255)                                                                                                                                                          gCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAgg                                                                                   >  1:73778/1‑251 (MQ=255)                                                                                                                                                             ggAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATc                                                                                >  1:254829/1‑251 (MQ=255)                                                                                                                                                                              aaCGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAgc                                                               <  1:102574/251‑1 (MQ=255)                                                                                                                                                                                       ggCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAATCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCAt                                                      <  2:240331/250‑1 (MQ=255)                                                                                                                                                                                                aCAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCATGgtg                                             >  2:262133/1‑251 (MQ=255)                                                                                                                                                                                                 cAGCCGCGGTCCGCCAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGGACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGgtgt                                            >  2:246471/1‑251 (MQ=255)                                                                                                                                                                                                           ccGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGACATGCATTGCCTGGTGTTATCGGTGAc                                  >  1:147277/1‑251 (MQ=255)                                                                                                                                                                                                            cGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGAt                                                                             >  2:150934/1‑207 (MQ=255)                                                                                                                                                                                                            cGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCTGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGGATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGAt                                                                             <  1:150934/207‑1 (MQ=255)                                                                                                                                                                                                             gCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGa                                <  1:96041/251‑1 (MQ=255)                                                                                                                                                                                                               aaGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATc                              <  1:293677/251‑1 (MQ=255)                                                                                                                                                                                                                                           agagCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAACACTAGGCAAGTTCATGCAGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTg  >  2:157703/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      GGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTG  >  NC_000913/90669‑91153

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