Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,152	0	A	G	22.8%	148.5 / 17.1	79	intergenic (+135/‑63)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (43/18);  new base G (0/18);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 3.30e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3781055‑3781250                                                                                                  |                                                                                                    gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑gg                                                                                                  >  1:593845/1‑100 (MQ=255) gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑gg                                                                                                  <  2:624632/100‑1 (MQ=255)  agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑ggg                                                                                                 <  1:1396575/100‑1 (MQ=255)  agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑ggg                                                                                                 <  2:1955447/100‑1 (MQ=255)   gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑gggg                                                                                                <  2:1953106/100‑1 (MQ=255)    agGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGa                                                                                               >  2:1927838/1‑100 (MQ=255)     ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGaa                                                                                              >  1:1273175/1‑100 (MQ=255)     ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGaa                                                                                              >  1:2037824/1‑100 (MQ=255)     ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGa                                                                                               >  2:253713/1‑99 (MQ=255)          ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGc                                                                                         >  2:2878475/1‑100 (MQ=255)          ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGc                                                                                         <  1:2040155/100‑1 (MQ=255)                ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACt                                                                                   >  1:2659654/1‑100 (MQ=255)                      gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCAc                                                                             >  2:1179850/1‑100 (MQ=255)                        atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCa                                                                           >  2:2773266/1‑100 (MQ=255)                         tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCaa                                                                          >  2:2868719/1‑100 (MQ=255)                            aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAAttt                                                                       >  2:188386/1‑100 (MQ=255)                               cccccccgctGAAACGCCACCGAACTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGAGGGCGGGGGGGG‑GGGGAAAAGGCCGCACTGCTCACCAAttt                                                                       <  2:390653/89‑1 (MQ=255)                               ccccCACTCTAAAACGACACCGCACTCTTTTTTCCCCCCTCGCCCCCCGGGGGAGAGGCCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAAtt                                                                        <  2:1275075/96‑1 (MQ=37)                               cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  1:1341434/1‑100 (MQ=255)                               cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  1:375348/1‑100 (MQ=255)                               cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  1:286973/1‑100 (MQ=255)                               cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  1:2341804/1‑100 (MQ=255)                               cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  1:1720682/1‑100 (MQ=255)                                cccccctcgGAAACGACACCGCTCTCTTTTTTTCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGTGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                   <  2:1038496/95‑1 (MQ=255)                                 cccACTCTGAAACGCCACCGCACTTTTTTTTTCCCCCCCGCCCCCCGGGGGGGAGGCCGGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                  <  2:375346/98‑1 (MQ=255)                                  tacctctGAAACGACCCCGCCTTTTTTTTTCCCCCCCCGCCCCTCGGGGGGGGGGGCGGGGGGGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                  <  1:1076042/96‑1 (MQ=255)                                  tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                 >  1:1903064/1‑100 (MQ=255)                                  cccctcgaaaaCGCCACCGCCCTTTTTTTTTCCCCCCCCCCCCTCGGGGGGGGGGGCGGGGGGGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                  <  1:2773282/91‑1 (MQ=255)                                   ccctctAAAACGCCCCCCCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGGGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                                                <  1:2294458/98‑1 (MQ=255)                                   cACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                                                >  2:707240/1‑100 (MQ=255)                                    aCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt                                                               >  1:2467207/1‑100 (MQ=255)                                         gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                          >  1:1628875/1‑100 (MQ=255)                                         gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                          >  1:2045014/1‑100 (MQ=255)                                         gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                          >  1:2981420/1‑100 (MQ=255)                                         aaaccgcccccGCCCTCTTTTTTCCCCCCCCCCCCCCCGGGGGAAGGGCCCGGGAGGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                  <  1:1814392/83‑1 (MQ=255)                                         aaaagaaccccGCTCTCTTTTTTCCCCCCCCCCCCCTCCGGGGAGGGGCCCGGGGTAG‑GGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                          <  1:1028941/91‑1 (MQ=255)                                           aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                        >  2:377398/1‑100 (MQ=255)                                             cgccacCGCTCTTTTTTTTTCTCCTCCGCCCTCCGGGGGAGGGGGCGGGGGTGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                       <  1:606716/96‑1 (MQ=255)                                             cGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                      >  1:2320135/1‑100 (MQ=255)                                              gACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccca                                                     >  2:2421633/1‑99 (MQ=255)                                                cacCGCACTCTTTTTTTCCCCCCCGCCCCCCGGGGGGGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                                    <  2:664878/99‑1 (MQ=38)                                                  ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                                 >  2:2788920/1‑100 (MQ=255)                                                  ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                                 >  1:1463504/1‑100 (MQ=255)                                                    gccctctTTTTTTCTCCCTCGCCCCCCCGGGGAGGGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                                                   <  2:2037812/93‑1 (MQ=255)                                                     cACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                              >  1:384383/1‑100 (MQ=255)                                                      aCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                             >  1:504147/1‑100 (MQ=255)                                                         ctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                          >  1:1310175/1‑100 (MQ=255)                                                         ctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                          >  2:960403/1‑100 (MQ=255)                                                              tttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCGGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                                     >  1:2085764/1‑100 (MQ=255)                                                              tttCCCCCCCCCCCCTCGGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                                     <  1:2786692/100‑1 (MQ=38)                                                                  ccgatctccccTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                                 <  2:1296855/93‑1 (MQ=255)                                                                    ccccGCCCCTCCGGGGAGGGGCCGGGGGTGG‑GGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                               <  2:1859796/97‑1 (MQ=255)                                                                    ccTCGCCCCTCGGGGGGGGGGGCCGGGGTGG‑GTGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                                 <  2:593842/98‑1 (MQ=255)                                                                      tCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                             >  1:2469556/1‑100 (MQ=255)                                                                       ccccccccGGGGGGGGGGGCGGGGGGGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                            <  1:565114/98‑1 (MQ=255)                                                                         ccccTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTagatcgg                          >  1:1296862/1‑94 (MQ=255)                                                                          cgcTCGGGGGGGAGGCCCGGGTGCG‑GGGGAAAAGGCCGCACTGCTCACCAAGTTTCTGCTATCCTGCCCCCTTAGTGCGGGGGCAGGAtagca                                <  2:1457657/91‑5 (MQ=12)                                                                          cccTCCGGGGAGGGGCCGGGGGTGG‑GGTGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACACCg                         <  1:377400/100‑1 (MQ=255)                                                                           ccTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                        >  1:1751649/1‑100 (MQ=255)                                                                            cTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                       >  2:1317695/1‑100 (MQ=255)                                                                              ccGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                     >  2:982123/1‑100 (MQ=255)                                                                               cGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt                    >  1:1522330/1‑100 (MQ=255)                                                                               aggggCGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt                    <  2:1848335/99‑1 (MQ=255)                                                                                  ggggggggCGGGGGGGG‑GGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                 <  1:707243/95‑1 (MQ=255)                                                                                  ggAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                 >  2:1356270/1‑100 (MQ=255)                                                                                  ggAGA‑GGCCGGGGGGAGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                 <  1:657643/100‑1 (MQ=255)                                                                                   gggggggCGGGGGTGG‑GGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg                <  1:253715/96‑1 (MQ=255)                                                                                       gggCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc            <  2:1341427/100‑1 (MQ=255)                                                                                       gggCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc            >  1:2371604/1‑100 (MQ=255)                                                                                        ggCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTG‑CCCCGCACTAAGGGGGCAGTATGCTAAACATCCTACTTTACGAACCAg          >  1:2938538/1‑100 (MQ=255)                                                                                         gCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg          >  2:2904903/1‑100 (MQ=255)                                                                                          cgggggTGG‑GGGGAAATGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa         <  2:1751848/98‑1 (MQ=255)                                                                                          ccGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa         >  2:794614/1‑100 (MQ=255)                                                                                           cGGGGTGG‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTATGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg          <  2:1600134/98‑1 (MQ=255)                                                                                            ggggggA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg          <  2:606116/97‑1 (MQ=255)                                                                                             gggTGG‑GGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt      <  2:2981404/100‑1 (MQ=255)                                                                                             gggTGA‑GGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt      <  2:1628866/100‑1 (MQ=255)                                                                                                 gA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCg  <  1:1179856/100‑1 (MQ=255)                                                                                                 gA‑GGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCg  <  1:1927849/100‑1 (MQ=255)                                                                                                  |                                                                                                    GAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGA‑GGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3781055‑3781250

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