Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	44.4%	67.3 / 17.7	55	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (7/23);  new base G (24/0);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 2.69e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925751‑926049                                                                                                                                                       |                                                                                                                                                         gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGTGGGGGa                                                                                                                                                          <  2:796092/150‑1 (MQ=255)   gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTc                                                                                                                                                       <  1:172480/151‑1 (MQ=255)      cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGTTC‑g                                                                                                                                                     >  1:168542/1‑149 (MQ=255)      cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑Gaa                                                                                                                                                   <  1:502064/151‑1 (MQ=255)         tGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GGGGGGGGGGtt                                                                                                                                                        >  1:229342/1‑144 (MQ=255)            cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCg                                                                                                                                            <  1:625360/151‑1 (MQ=255)                     aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTTG‑GAA‑ccccccgacgtt                                                                                                                                      >  1:19495/1‑142 (MQ=37)                               ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTT‑GGA‑CCCCCCAAACGTTGCCGta                                                                                                                               >  1:863903/1‑145 (MQ=255)                               ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GTGGGGGGGGGTT‑GCA‑CACCCCCTACGTAGCCGTCTac                                                                                                                            >  1:643932/1‑148 (MQ=25)                                 cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTC‑GAA‑CCCCCCCTACGTTTTCGta                                                                                                                               >  2:295646/1‑143 (MQ=25)                                  acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGTTT‑GAA‑ACCCCCCTACGTTGTCGGATacacac                                                                                                                        >  1:640919/1‑149 (MQ=25)                                  acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGTTC‑GAA‑CCCCCCCTACGTTTCCGTATACACACtt                                                                                                                      >  1:320236/1‑151 (MQ=255)                                   caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GAGGGGGGGGTTT‑GAA‑ACCCCCATAAGTTGCCGTATaaaccactt                                                                                                                     >  2:140327/1‑146 (MQ=25)                                     ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTcc                                                                                                                   <  2:126940/151‑1 (MQ=255)                                        aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAgg                                                                                                                <  1:42619/151‑1 (MQ=255)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGTTC‑CAA‑ACCCCCATACGTTTCCGTATACACAATTTCCCgg                                                                                                                >  1:348910/1‑150 (MQ=21)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTT‑GAA‑CCCCCCATACGTTGCCGGATACACACttttccgg                                                                                                                >  1:972242/1‑143 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCtt                                                                                                      <  1:228923/150‑1 (MQ=255)                                                               gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                         <  1:89454/151‑1 (MQ=255)                                                                 tCACTGCCGTACAGACAGATAAAATTCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCGTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                        <  1:81209/150‑1 (MQ=255)                                                                   aCTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GATGGGGGGGATA‑GAA‑ACCCCCCTAAATTGTCCTATAACAACCTTCCCGGCGGGCTCCTTCCGCCAcccaggccc                                                                                       >  1:729316/1‑141 (MQ=12)                                                                    cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCCCTACGTTTCCGTATACACACTTTCCAGGCGGGCTCCTTCAGCCACTCGGACAcc                                                                                      >  1:325119/1‑149 (MQ=255)                                                                     tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GTGGGGGGGGGTCTGAA‑CCCCC‑CTACCTTGGCCCATCCACACCTTCCCGGGGGGCTCCCTCAGcccccccggcacccc                                                                                    >  1:983407/1‑137 (MQ=12)                                                                       ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGTTT‑GAA‑ACCCCGATACGTTGCCGTATACCCACTTTTCAGGCGTGGTCCTTCCGCCACTCCGGCACCCCAcc                                                                                 >  1:267401/1‑151 (MQ=17)                                                                       ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTC‑GAA‑CCCCCCATACGTTGCCGTATACACACTTTTCAGGCGGGCTCCTTCAGCCACTCCGACACCTCAcc                                                                                 >  1:485809/1‑151 (MQ=17)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATCCACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCaaa                                                                              <  2:869896/150‑1 (MQ=255)                                                                            cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTC‑GCA‑CCCCCCATACGTTGCCGTATACACACTTTCCCGGAGTGGTCCTTCAGCCCCCCCGAGACCTCACCcaa                                                                              >  1:124900/1‑149 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GAGGGGGGGGTTT‑CAA‑CCCCCCATACGTTTCCGTATTAAAACTTTTCAGGCGTGCTCCCTCACCCCCTCCGGCACCCCCCCCAAAtttttt                                                                       >  2:220621/1‑151 (MQ=9)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GAGGGGGGGGGTT‑GAA‑CCCCCCCTACGCTTCCGTATACACACTTTTCAGGGGGGGTCCTTCATCCCCCCGGACACCCCACCCAAATGtttt                                                                       >  2:424136/1‑151 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GAGGGGGGGGGTC‑GAA‑ACCCCCAAACGTTGCCGTATACACACTTTCCAGGGGGGGCCCCTCTCCCCCTCCGACACCCCCCCCAAATGtttt                                                                       >  2:479577/1‑151 (MQ=14)                                                                                  gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GAGGGGGGGGTTC‑GAA‑CCCCCCACACGTTGCCGTGTACAAACTTTCCCGGCGGGCTCCTTCCGCCCCCCGGGCCCCCCCCCAAAATGtttt                                                                       >  1:104317/1‑150 (MQ=14)                                                                                   aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GAGGGGGGGGGTT‑GAA‑ACCCCCATAAGTTTGCGGATAAAAACATTTCAGGGGGGGGCCTTTCGCCAAATGGAAACCTCACCAAAAtttttt                                                                       >  2:521448/1‑149 (MQ=9)                                                                                    tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGAGTC‑TGA‑ACCCCGATACGTTGCCGTATACACACTTTTCAGGCGTGGCCCTTCCTACCCCCGGACACCCCACCCAAAAGttttttt                                                                    >  1:539937/1‑148 (MQ=14)                                                                                     aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTa                                                                   <  1:169835/151‑1 (MQ=255)                                                                                                aaaaaGCTCGCCCTTTCGTACGCGCTCTTCTTTAAAAATGGCGGT‑GACGGGGCCACTC‑GAA‑CCCCCTATACGTTCCGGTATACACACTTTCCAGGCGTTCTCCTTCAGCCTCTCGGACACCTCGCCAAATTGTTTTGCTACCAATCCT‑CAt                                                       <  2:257413/151‑1 (MQ=255)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATg                                                      <  1:921747/151‑1 (MQ=255)                                                                                                  aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GAGGGGGGGGGTC‑TTA‑CCCCCCCTACGTTTTCCTCTACACACTTTCCAAGGGTGGTTCCTCAGCCACTCGGGCACCCCACCCCACTGTTTTGTTACCCAAACT‑CCTgg                                                     >  2:370551/1‑151 (MQ=9)                                                                                                  aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTT‑GGA‑CCCCCCAGACGTTTTCGGATAAACACTTTTCAGGCGGGCTCCCTCAGCCACTCGGACACCCCCCCCAATttttctttcttgcttt                                                             >  2:471266/1‑132 (MQ=2)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGc                                                <  1:836449/149‑1 (MQ=255)                                                                                                            cTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACg                                            <  2:404918/150‑1 (MQ=255)                                                                                                              ttCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACggg                                          <  1:649520/150‑1 (MQ=255)                                                                                                               tCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGTTT‑GAACCCCCCGCTAC‑TTTCCGTATAAACCATTTCCCAGCGGGGTCCTTCAGCCCCTCGGGCACCCCCCCAAAATGTTTTTGTACCAAACCTACATGGG‑GGCAACgagg                                         >  2:766393/1‑150 (MQ=9)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGG‑TTC‑GAA‑CCCCCCATACGTTGCCGTATACACACTTTTCAGGCGTGCTCCTTCAGCCACCCGGACACCTCACCAAAATTTTTTTTTACCCAACCC‑CCTGGGGGGGAACgggggg                                       >  1:13705/1‑147 (MQ=11)                                                                                                                     gAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGTTT‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGGCACCCCACCAAATTGTTTTGCTACCAAAACT‑CATGGGGGGGAACGGGGCGccact                                  >  1:371991/1‑151 (MQ=17)                                                                                                                             cttTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAg                              <  1:328745/147‑1 (MQ=255)                                                                                                                                tAAATATGGCGGT‑GATGGGGGGATTC‑GAA‑ACCCCCATACGTTGCCGTATCCACCCTTTCCCAGCGTGCGCCTTCAGCCACTCCGACAACTCCCCCAATTTTTTTTCTACCACACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGGGtt                        >  1:202877/1‑150 (MQ=255)                                                                                                                                tAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                        <  2:604314/150‑1 (MQ=255)                                                                                                                                tAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCCATACGTTTCCGTATACACACTTTCCAGGCGGGGTCCTTCAGCCACTCGGACACCTCACCAAATTTTTTTGCTACCAAACCT‑CATTGGGGGCAACGGGGCGCTACTATAGGGAGtt                        >  1:163781/1‑150 (MQ=18)                                                                                                                                     atGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                   <  2:320186/150‑1 (MQ=255)                                                                                                                                         cGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc              <  1:34102/151‑1 (MQ=255)                                                                                                                                         cGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc              <  1:340621/151‑1 (MQ=255)                                                                                                                                         cGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc              <  1:92805/151‑1 (MQ=255)                                                                                                                                              gAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt           <  2:497501/150‑1 (MQ=255)                                                                                                                                               aGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc          >  1:406615/1‑150 (MQ=255)                                                                                                                                                gggggggATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa         >  1:612430/1‑150 (MQ=255)                                                                                                                                                      gATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa  <  2:693889/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                         GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGATTC‑GAA‑CCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCT‑CATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925751‑926049

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