Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	220,004	C→A	9.1%	Q8504K (CAG→AAG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	220,004	0	C	A	9.1%	164.4 / 3.0	58	Q8504K (CAG→AAG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base A (2/4);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGCAAAGGTAACGACATCCTG  >  NC_002947/219857‑220150                                                                                                                                                      |                                                                                                                                                   aTTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCa                                                                                                                                                   <  5:43045/149‑1 (MQ=255) aTTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCa                                                                                                                                                   <  4:133751/149‑1 (MQ=255)    gATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGt                                                                                                                                                 <  6:350405/148‑1 (MQ=255)      tAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACa                                                                                                                                              >  6:392962/1‑149 (MQ=255)       aGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACAt                                                                                                                                             >  6:294736/1‑149 (MQ=255)           gTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGTGGTCGTCAGCGAGGGCTACCAGGCCATGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAATTCGGTGCTGTCAATACAAGCAATGTTCACCAGTACATAAc                                                                                                                                          >  1:102523/1‑149 (MQ=255)           gTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCAcc                                                                                                                                         <  7:19228/149‑1 (MQ=255)           gTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCAcc                                                                                                                                         <  8:212731/149‑1 (MQ=255)              aCGGCAACGACATCATCTTCGGTGACCTGATCACGCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGa                                                                                                                                       <  7:196975/148‑1 (MQ=255)              aCGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAg                                                                                                                                      <  2:200533/149‑1 (MQ=255)                ggCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCa                                                                                                                                    <  4:246114/149‑1 (MQ=255)                 gCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCAc                                                                                                                                   <  4:90275/149‑1 (MQ=255)                                 cGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGa                                                                                                                     <  1:4051/147‑1 (MQ=255)                                           aTCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTATTTTTTTTTTTGATAAGTAGATCACCGAGATCTACACCGAGTTCGATATCTCCGGt                                                                                                           <  2:399042/147‑1 (MQ=255)                                           aTCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCAATTTTTTTTTTTTTGATAAGTAGATCACCGAGATCTACACCGAGTTCGATATCTCCGGt                                                                                                           <  8:214235/147‑1 (MQ=255)                                               cccTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGGAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAg                                                                                                     <  3:160826/149‑1 (MQ=255)                                               cccTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACAAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAg                                                                                                     >  7:344041/1‑149 (MQ=255)                                                ccTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAg                                                                                                     >  4:180860/1‑148 (MQ=255)                                                 cTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGa                                                                                                   >  5:46431/1‑149 (MQ=255)                                                  tCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGAc                                                                                                  >  3:266501/1‑149 (MQ=255)                                                  tCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCAACGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGAc                                                                                                  >  5:231583/1‑149 (MQ=255)                                                    aaCGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAAGGTACACCAGTACATCACCGAGAACTACAACGAGTTCGATATCTCCGGTGCCAAGGACgg                                                                                                >  5:310504/1‑149 (MQ=255)                                                       ggT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCaa                                                                                             <  7:286833/149‑1 (MQ=255)                                                        gt‑gtCGTCAGCGATGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAAt                                                                                            <  1:431353/149‑1 (MQ=255)                                                            tcgtcAGCGAGTGCTACCAGTCCCTGCAGACGTATGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTTTCAATTATATTATTTTTGATAAGTAGATCACCGAGATCTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGa                                                                                          <  4:270431/148‑1 (MQ=255)                                                              gtcAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACa                                                                                        <  3:180860/148‑1 (MQ=255)                                                              gtcAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACAt                                                                                       >  3:357390/1‑149 (MQ=255)                                                               tcAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAGTGGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACAt                                                                                       >  1:308376/1‑149 (MQ=255)                                                                cAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACAcc                                                                                                                             >  6:299705/1‑109 (MQ=255)                                                                cAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACAcc                                                                                                                             <  5:299705/109‑1 (MQ=255)                                                                 aGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTAGAACTCGGTGCTGTCACTACCATCAAAGTTCACAAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCt                                                                                    >  2:226894/1‑149 (MQ=255)                                                                  gCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGACGGTGCTGTCACTACCACCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATATCCGGTGCCa                                                                                                       >  1:199382/1‑129 (MQ=255)                                                                    gAGGGCTACCATGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTc                                                                                 <  5:127201/149‑1 (MQ=255)                                                                                gCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGAAATACTGTCGGGtggcaatgg                                                                     >  3:405310/1‑149 (MQ=255)                                                                                gCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGAGATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGtggcaatgg                                                                     >  7:19216/1‑149 (MQ=255)                                                                                     gCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTggcaatggcaat                                                                 >  1:319645/1‑148 (MQ=255)                                                                                      cAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGTTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTggcaatggcaat                                                                 <  7:322962/147‑1 (MQ=255)                                                                                      cAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTggcaatggcaat                                                                 <  5:294736/147‑1 (MQ=255)                                                                                      cAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTggcaatggcaat                                                                 <  5:392962/147‑1 (MQ=255)                                                                                           gtacgtTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGAGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATc                                                           >  5:161125/1‑148 (MQ=255)                                                                                              cgtTTCACAGAAG‑ATCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCctc                                                        <  5:409816/148‑1 (MQ=255)                                                                                              cgtTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCctct                                                       >  2:414646/1‑149 (MQ=255)                                                                                                     caGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTAACATAAATGTTCAACAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACTGCAATGACATCCTGTCGGGTGGCAATGGCAATTACATCCTATTCGGcc                                                 >  4:396456/1‑148 (MQ=255)                                                                                                      agaag‑aGCGGTGTTGAAGTCGGTGCTGTCACTAATATTATTTTTGATAAGTAGATCACCGAGAACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCaa                                               <  2:102523/149‑1 (MQ=255)                                                                                                               ggTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAg                                       <  7:341950/149‑1 (MQ=255)                                                                                                                     gAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGacac                                 >  8:409196/1‑149 (MQ=255)                                                                                                                                tGTCACTACAATCAATTTTCACCAGTAGATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAAtggtg                      <  6:310504/149‑1 (MQ=255)                                                                                                                                  tCACTATCAGTAATTTTCAACAGTAGATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTggcaatggcaat                                                                 <  4:313026/104‑1 (MQ=255)                                                                                                                                  tCACTACCAGCAATTTTCACCAGTACATCACCGAGAACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAAtggtg                      <  2:392113/147‑1 (MQ=255)                                                                                                                                  tCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTggcaatggcaat                                                                 >  3:313026/1‑104 (MQ=255)                                                                                                                                   cACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGc                    >  2:434191/1‑148 (MQ=255)                                                                                                                                     gtgacAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGCaaa                 <  3:334192/145‑1 (MQ=255)                                                                                                                                      tACCAGCAATGTTCACCAGTACATAACCGAGCAATACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGCAAAg                >  7:104499/1‑149 (MQ=255)                                                                                                                                          aGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGCAAAGGTaa            >  1:361421/1‑149 (MQ=255)                                                                                                                                              aTGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGATGAATGGTGGCAAAGGTAACGa         >  3:292565/1‑148 (MQ=255)                                                                                                                                                gTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGCAAAGGTAACGACAt      >  8:57885/1‑149 (MQ=255)                                                                                                                                                    aCCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACa                                                                                        >  4:222923/1‑63 (MQ=255)                                                                                                                                                    aCCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGCAAAGGTAACGACATCctg  >  6:442688/1‑149 (MQ=255)                                                                                                                                                      cAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCg                                                    >  4:47117/1‑97 (MQ=255)                                                                                                                                                      cAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCg                                                    <  3:47117/97‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ATTGATAGTGGTAACGGCAACGACATCATCTTCGGTGACCTGATCACCCTCAACGGT‑GTCGTCAGCGAGGGCTACCAGGCCCTGCAGACGTACGTTGCACAGAAG‑AGCGGTGTTGAAGTCGGTGCTGTCACTACCAGCAATGTTCACCAGTACATCACCGAGCACTACACCGAGTTCGATATCTCCGGTGCCAAGGACGGCAATGACATCCTGTCGGGTGGCAATGGCAATGACATCCTCTTCGGCCAAGGTGGTAGCGACACGCTGAATGGTGGCAAAGGTAACGACATCCTG  >  NC_002947/219857‑220150

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