Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,816	0	G	A	43.4%	36.5 / 51.2	53	intergenic (+75/+144)	yneM/dgcZ	inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (29/1);  new base A (0/23);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 3.85e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACG  >  NC_000913/1622668‑1622955                                                                                                                                                              |                                                                                                                                            ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑Ga                                                                                                                                            >  2:44708/1‑150 (MQ=255)          tACTGGGAATAATTTTATTTTCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCg                                                                                                                                  >  1:541351/1‑151 (MQ=255)              gggAATAATTTTATTTTCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCa                                                                                                                               <  2:981016/150‑1 (MQ=255)                   taatTTTATTTTCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAAc                                                                                                                         >  2:805358/1‑151 (MQ=255)                    aatTTTATTTTCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAAcc                                                                                                                        >  2:31129/1‑151 (MQ=255)                      ttttattttCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACcgg                                                                                                                      >  2:254834/1‑151 (MQ=255)                       ttttttttCTGTTTTTC‑GGGCCGCTTTTTCCACCCAAAAAGGGG‑ATGATTAATGAA‑GGGAAATAACCCCCCCCCTACCCGCCCCT‑TTGTT‑ACGGTTGTGA‑AAAGGGGGC‑T‑TGTTTTTTAGGAGGGA‑AAAAAAAAACCGCCAGTAAACcgg                                                                                                                      <  2:701421/150‑1 (MQ=255)                       tttattttCGGTTTTTC‑GGCCCGTTTTTTCCAGCCACAATGGGG‑ATGAATAGGGAA‑CGAAGATTCCCCCCCACCTACCCGCCCGT‑TTATA‑ACAGTTTTGA‑ACAGGGAGC‑T‑TGTTTTTTAGACCGGC‑AAAAAAAAACCGCCAGTAAACcgg                                                                                                                      <  1:601560/150‑1 (MQ=255)                          tttttttggttttcgtgCCCGCGTTTTCCACCCAAAAAGGGG‑ATGACTAATAAA‑GGGAAATACCCCCTCCCTTACCCGCCCCA‑TTATT‑ACA‑TTGGAG‑AAAAGGGGC‑CGTTTTTTTTATGACGGA‑AAAAAAAAACCGCCAGTAAACcggcg                                                                                                                    <  1:594939/149‑1 (MQ=255)                            tttttttttttc‑gggCCGCTTTTTCCACCCAAAAAGGGG‑ATGATAAAAGAG‑GGAAAAAACCCCCTCCCCAACCCGCCCCT‑TTTTT‑ACAGTTGTGT‑AAAGGGGGC‑T‑TGTTTTTTAGGAGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                 <  2:347712/143‑1 (MQ=255)                             ttttgtttttc‑gggCCGCTTATTCCACCCA‑AAATGGGAATGA‑TAAAGAA‑CGGGAAAAACCCCCTCCCCAACCGGCCCC‑TTTTT‑AAAATT‑TGTGAAAAGGGGCTC‑TGTTTTTTATGCGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGa                                                                                                                <  1:277811/144‑1 (MQ=255)                               ctggttttc‑tggCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                             >  2:991071/1‑151 (MQ=255)                                 ggttttc‑tggCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt                                                                                                           >  2:607236/1‑151 (MQ=255)                                   ttttc‑tggCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                         >  2:833407/1‑151 (MQ=255)                                   ttttc‑tggCCGCGTATTCCAGCAAAAAAGGGG‑ATAAATAAGAAA‑GGGAGAAAACCCCCCCCCTACCCGCCCCT‑TTTTT‑ACATTTGTGA‑AAAGGGGGC‑T‑TGTTTTTTAGGGGGAC‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                          <  2:895782/150‑1 (MQ=255)                                       c‑tggCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                     >  1:281877/1‑151 (MQ=255)                                         tggCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                    >  1:409331/1‑151 (MQ=255)                                               cgTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                   >  1:1097970/1‑146 (MQ=255)                                               cgTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                  >  2:733328/1‑147 (MQ=255)                                                        gCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtttt                                                                                     >  2:418810/1‑151 (MQ=255)                                                           acaAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCt                                                                                  >  1:425986/1‑151 (MQ=255)                                                            caAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                 >  2:1038223/1‑151 (MQ=255)                                                             aaggggg‑ATGACTAATAAA‑GGGAGATAACCCCTCCCCTACCCGGCCCT‑TTTTA‑AATTTTGTGT‑AAAGGGGGC‑C‑TGTTTTTTAGGAGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                <  2:352954/147‑1 (MQ=255)                                                             aaaTGGG‑TTGA‑TAATTAAGCGGAGAAAACCCCACCCCAACCGGGCCCT‑TTTTT‑ACTGGGTTAC‑ACAAGGGCC‑T‑TTTTTTTTGTGAGGGA‑AAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                    <  1:788138/147‑1 (MQ=255)                                                              aaTGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                               >  2:101926/1‑151 (MQ=255)                                                                  gg‑ATGACAAATGAC‑GGGAAATACCCCCCCCCCAACCCGCCCCT‑TTGTT‑AAGTTTGTGA‑ACAGGGGCC‑T‑TGTTTTTTAGGAGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGc                                                                           <  2:998960/151‑1 (MQ=255)                                                                            aTGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc                                                                   >  2:1083721/1‑150 (MQ=255)                                                                            aTGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAAAGGc                                                                   >  1:520962/1‑150 (MQ=255)                                                                                aagggataaaaccccTCACC‑AACCGGCCCCTTTGTG‑AAATTTGTGT‑AAAAGGGCC‑T‑TGTTTTTTATGAGGGC‑AAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                               <  1:713321/139‑1 (MQ=255)                                                                                   ggaaaaaaccccTCCCCTACCCGCCCCT‑TTTTT‑ACATTTGTGA‑ACAAGGGGT‑T‑TGTTTTTTGGGGCGGG‑AAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                              <  1:719117/140‑1 (MQ=255)                                                                                       atacccccTCCCCTACCCGCCCCT‑TTGTT‑AATGTTGTGT‑AAAGGGGGC‑T‑TGTTTTTTAGGCGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                         <  2:341437/145‑1 (MQ=255)                                                                                        tAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                       >  2:34834/1‑151 (MQ=255)                                                                                        ataccccccccccaaCCCGCCCC‑TTTTTAACATTTGTGA‑AAAGGGGGC‑T‑TGTTTTTTAGGCCGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                           <  2:69064/135‑1 (MQ=255)                                                                                                  ccgttctgcccct‑TTTTT‑ACATTTGTGT‑AAAGGGGGC‑T‑TGTTTTTTAGGACGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                   <  2:191565/136‑1 (MQ=255)                                                                                                     acccGCCCCT‑TTTTT‑ACAGTTGTGA‑ACAGGGGGG‑T‑TTTTTTTTAGGAGGAA‑AAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                            <  1:126177/147‑1 (MQ=255)                                                                                                       ccccccct‑TTTTT‑ACAGTTGTGT‑AGAGGGGCC‑T‑TTTTTTTTATGCCGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                <  2:972532/139‑1 (MQ=255)                                                                                                         ggCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                   >  1:191531/1‑138 (MQ=255)                                                                                                             cc‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTTCTGCACTGATAACGAATCGTTGacac                                   >  1:108539/1‑150 (MQ=255)                                                                                                                ttttt‑ACGGTTGTGA‑AAAGGGGGC‑T‑TTTTTTTTAGGAGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca                                  <  2:670360/146‑1 (MQ=255)                                                                                                                ttGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAg                                 >  1:191602/1‑150 (MQ=255)                                                                                                                 tGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                >  1:516319/1‑150 (MQ=255)                                                                                                                   tt‑ACAGTTGTGT‑ACAGGGGGC‑C‑TTTTTTTTATGAGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                               <  2:600469/149‑1 (MQ=255)                                                                                                                   tt‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAg                              >  1:752649/1‑150 (MQ=255)                                                                                                                        aaagtggtgaacaGGGGC‑C‑TTGTTTTTTAGGCGGCAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCa                            <  2:152831/138‑1 (MQ=255)                                                                                                                          ttGTGA‑ACAGGGGCC‑T‑TTTTTTTTAGGAGGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                        <  2:281932/150‑1 (MQ=255)                                                                                                                           tggtaaaaagggggC‑T‑TGATTTTTATGGCGAA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                        <  1:255483/140‑1 (MQ=255)                                                                                                                                   aaGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                >  1:1087820/1‑150 (MQ=255)                                                                                                                                    agaggC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGCAAACCGGCGGTGAATGCCTGGATGGATAGAATTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat              >  1:401259/4‑151 (MQ=255)                                                                                                                                     aagggcctttttttttATGACGGC‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                  <  2:723487/137‑1 (MQ=255)                                                                                                                                       ggC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga            >  2:1050927/1‑150 (MQ=255)                                                                                                                                           ctttttttttAGGCCGGA‑AAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat          <  2:881546/145‑1 (MQ=255)                                                                                                                                              gATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt       >  1:167487/1‑150 (MQ=255)                                                                                                                                                  tttATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACg  >  2:539548/1‑151 (MQ=255)                                                                                                                                                              |                                                                                                                                            TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTC‑TGGCCGCGTATTTCAGCCACAAATGGG‑ATGACTAATGAA‑CGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAGTTGTGT‑ACAAGGGGC‑C‑TGATTTTTATGACGGC‑GAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACG  >  NC_000913/1622668‑1622955

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