Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,0040CT68.3% 32.6 / 32.1 41intergenic (‑111/+180)borD/ybcVprophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (13/15);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTC  >  NC_000913/578760‑579220
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
ccAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCAGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGAAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTg                                                                                                                                                                                                                    >  1:401738/1‑251 (MQ=255)
  ngAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCg                                                                                                                                                                                                                  >  1:266372/2‑251 (MQ=255)
                gatgGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGGTTTGAACTGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAGTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTGAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTg                                                                                                                                                                                                    <  1:204364/251‑1 (MQ=255)
                                  ttGGTGCTACTTCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAATCCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTa                                                                                                                                                                                   <  1:393388/250‑1 (MQ=255)
                                             gctgTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGcc                                                                                                                                                                        <  1:141702/250‑1 (MQ=255)
                                                    gtttgtttTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAAGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAgg                                                                                                                                                                                                                      <  1:400033/197‑1 (MQ=255)
                                                        gtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAgn                                                                                                                                                                         <  1:269147/238‑2 (MQ=255)
                                                                    aaacgTCTGTTGAGCTCATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTGTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGAGTAACTATGTCCACAGCCCGGACTGGGAACTTCTCTGCggg                                                                                                                                                 <  1:140163/246‑1 (MQ=255)
                                                                            gTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGa                                                                                                                                                              >  1:374105/1‑229 (MQ=255)
                                                                            gTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGAGTAGAAAATTAAACAAACCCTAAACAAGGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTcc                                                                                                                                          <  1:431971/249‑1 (MQ=255)
                                                                            gTTGAGCACATCATGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGa                                                                                                                                                              >  1:374110/1‑229 (MQ=255)
                                                                             ttGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTGTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCgg                                                                                                                                        <  1:149414/250‑1 (MQ=255)
                                                                                       ccTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACCAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTACTGTATGCCAGGTGCGATGaa                                                                                                                                                                                                             >  1:61829/1‑171 (MQ=255)
                                                                                          gTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTaaaaa                                                                                                                           >  1:191627/1‑250 (MQ=255)
                                                                                            aataaGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGCTGCTCTTTGAAGTTCGCAGAATCGTCTGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACgg                                                                                                                                                                 >  1:242128/1‑210 (MQ=255)
                                                                                                      gCCAGCGCAGTATCGAGTAGCATTTTTTTCATGGTGTTAGTCCCGATGCTTTGTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTTTATTCCCTGATTAACGATTTCCACAGCCCTGACGGGTAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAggg                                                                                                             <  1:395352/251‑1 (MQ=255)
                                                                                                      gCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAggg                                                                                                             >  1:29770/1‑251 (MQ=255)
                                                                                                        cAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCNGAANCGTANGTNTAGAAAATTAAACAANCCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAANNTNNNTNATGTGCGa                                                                                                                                                                                                                 <  2:318769/150‑1 (MQ=255)
                                                                                                        cAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGa                                                                                                                                                                                                                 >  1:318769/1‑150 (MQ=255)
                                                                                                                   ngAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAg                                                                                                        >  1:358939/2‑243 (MQ=255)
                                                                                                                       tAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTAATAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATAAACTATGTCCACAGCCCTTACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAgg                                                                                                              <  1:220409/233‑1 (MQ=255)
                                                                                                                                   cATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGc                                                                                 >  1:9627/1‑251 (MQ=255)
                                                                                                                                    aTGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATAGTGTTAATATTTAATAATGTATGTCAGGTGCGATTAATCGTCATTGGATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCTTACACGTATTGCATTATGCCa                                                                               <  1:442774/251‑1 (MQ=255)
                                                                                                                                      ggTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCgg                                                                                                                                        <  1:409940/193‑1 (MQ=255)
                                                                                                                                                    ngCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAAc                                                                                                                                                            >  1:266336/2‑159 (MQ=255)
                                                                                                                                                     gCTTTTTGAAGTACGCAGAATCGTAGGTTTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAAATTCATATTGTGAATATTTATTAAGGTATGTCAGGGGCGATGAATCGTCATTGTATTCCCGGATTAACTATTTCCACAGCCCTTTGGGGGAACTTCTCTTCGGGAGTGTCCTGTAAAAATTAAAACGATGCAAGCAGGGTTAAGCGCGTGCGCGTATTGCATTATGCCAACGCCCCGGTGCTg                                                                 <  1:395353/248‑1 (MQ=255)
                                                                                                                                                      cTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCggg                                                                                                                                       <  1:355148/178‑1 (MQ=255)
                                                                                                                                                         nnnnnnnnnngCNNNNNCGTNTGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAAc                                                                                                                                                            <  2:266372/144‑1 (MQ=255)
                                                                                                                                                                   gCNNNNNCGTNTGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAAc                                                                                                                                                            <  2:266336/144‑1 (MQ=255)
                                                                                                                                                                   gCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATTTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCgg                                                 >  1:405546/1‑250 (MQ=255)
                                                                                                                                                                          ngTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTAAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCt                                                                  >  1:80389/2‑226 (MQ=255)
                                                                                                                                                                                tgtAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTAGTAATGTATGTCAGGTGCTATGAATCGTCAGTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAgtgt                                                                                                                                            <  1:411309/147‑1 (MQ=255)
                                                                                                                                                                                    gAAAATTAAACAAACCCTAAACACTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACg                                                            >  1:398070/1‑222 (MQ=255)
                                                                                                                                                                                     aaaaTTAAACAACCCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCGGACGGGGAACTTCTCTTCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCAGTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTgg                              <  1:62363/251‑1 (MQ=255)
                                                                                                                                                                                                      aaaCAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGGTATGATTTAGCGTGGAAAGATGTGTGTAgtg              <  1:342826/250‑1 (MQ=255)
                                                                                                                                                                                                              aGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGCATTAACTATGTCCACAGCCCTGACggg                                                                                                                                                                >  1:219010/1‑97 (MQ=255)
                                                                                                                                                                                                               gTTGAAATTTCATATTGTTAATCTTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAAt      >  1:36478/1‑249 (MQ=255)
                                                                                                                                                                                                                  gAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctn  <  1:87026/251‑2 (MQ=255)
                                                                                                                                                                                                                              tGTTAATATTTATTAATGTATGTCAGNTGCGATNNATNGTCATTGTATTCCCGGATTNACNNNGNNNANCGCCCNGANNGGGAACTTCNCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAgcg                                                                                                      >  2:102040/1‑138 (MQ=255)
                                                                                                                                                                                                                              tGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAgcg                                                                                                      <  1:102040/138‑1 (MQ=255)
                                                                                                                                                                                                                                                 attccAGGTGGGATGAATCGTCATTGTAGTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTtc                                                                                                                                                         <  1:216683/66‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTGCCNGGNGNGANNANTCGTNATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGNGNa                                                                                                                                                              >  2:216683/1‑64 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
CCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTC  >  NC_000913/578760‑579220

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