Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,4530AC22.6% 104.5 / 25.3 62intergenic (‑97/+19)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base C (8/6);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCC  >  NC_000913/4365315‑4365586
                                                                                                                                          |                                                                                                                                     
caccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                   >  2:560737/1‑143 (MQ=255)
  ccgcGACAGACGCGGTATTCGAACTTTGTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGccc                                                                                                                                 <  2:997930/140‑1 (MQ=255)
    acGACAGGCTCGGTCCTCGAACTTTTTTCTCCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTAGACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcg                                                                                                                               <  1:322060/143‑1 (MQ=255)
    acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGcaatg                                                                                                                               <  1:745751/143‑5 (MQ=255)
     cGACAGACGCGGTATTCGAACTTTTTTCATAAAGCATAATCACATCGCAATGATTCTTTCGTCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgc                                                                                                                              <  1:1004821/143‑1 (MQ=255)
       aCAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctt                                                                                                                            <  2:12007/143‑3 (MQ=255)
        cAGACGCGGTATTCGAACTTTTTGCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAt                                                                                                                           <  2:803450/143‑1 (MQ=255)
            cgcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGAt                                                                                                                       <  1:769181/143‑1 (MQ=255)
                gTATTCGAACTTTTTTCATCAGGCATAAACACCTCGCAAGGACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                   <  2:320060/143‑1 (MQ=255)
                 tATTCGAACTTTTTTCATCAAGCACCAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGCACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                                                                  <  1:32955/143‑1 (MQ=255)
                    tCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGCAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATTAccg                                                                                                               >  1:24087/1‑141 (MQ=255)
                            gtttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                       <  1:220868/142‑1 (MQ=255)
                                  tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGATCATTACAGCATGAAGCTACCCa                                                                                                 >  1:649612/1‑143 (MQ=255)
                                  tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTGAGTCGGCGGAAAAGAAAGCGGAGCGTGAACGCCCCGCGATTGACCACTACAGCAACAAGCTACCCa                                                                                                 >  2:327166/1‑143 (MQ=255)
                                  tcaAGCATAAACAACTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGCGACCGCCCCGCGACTGATCATTACAGCATGAAGCTACCCa                                                                                                 >  2:364020/1‑143 (MQ=255)
                                    aaGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTCAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                                                               >  1:319029/1‑143 (MQ=255)
                                     aGCATAAACACCCCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCa                                                                                              <  1:983160/143‑1 (MQ=255)
                                      gCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCa                                                                                              >  1:1011187/1‑142 (MQ=255)
                                      gCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCa                                                                                              <  2:1011187/142‑1 (MQ=255)
                                         tAAACGCCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaa                                                                                          <  1:688560/143‑1 (MQ=255)
                                             cacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGATCATTACAGCATGAAGCTACCCAGCAcgaccccg                                                                                      >  2:949198/1‑138 (MQ=255)
                                               ccgcgcgAGACCCGACGGAACCGAAAGCGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaa                                                                                    <  2:739400/140‑1 (MQ=255)
                                               ccTCGCTAGAACAGATGGAACCGCAAGAGAGCGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaa                                                                                    <  2:266535/143‑1 (MQ=255)
                                               ccTCGCAAGGACAGATGGAACCTCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaa                                                                                    <  1:101656/143‑1 (MQ=255)
                                                cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                   <  1:755909/143‑1 (MQ=255)
                                                  cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGCTCCTTACAGCATGAAGCCACCCAGCAcgacaccgagcca                                                                                 >  2:545996/1‑139 (MQ=255)
                                                   gCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                           >  2:312449/1‑97 (MQ=255)
                                                   gCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                           <  1:312449/100‑4 (MQ=255)
                                                          cAGATGGAACCGCAAGGGAGAGGTCACTATACAACGGCCGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCa                                                                         <  1:423726/143‑1 (MQ=255)
                                                            gagggAACCGCAAGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                                                       <  1:344317/140‑1 (MQ=255)
                                                            gATGGCACCGCAACAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                                                       <  1:9254/143‑1 (MQ=255)
                                                             gtgaaaCCGCCAGGTAGAGGTCACCATAGAAGGGAGGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCCTTACCGCCTGATCCTTCCTAGCACGAAGCCGAAGCAAACGGCCAggg                                                                      <  2:312625/139‑1 (MQ=255)
                                                              tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                     <  1:560737/143‑1 (MQ=255)
                                                                     gCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                              >  1:385565/1‑143 (MQ=255)
                                                                                gTCACTATACAACGGACGGGGGACGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCAGCGTCCGGTCAGTACAGCATGGAGCTACCCAGCACGCAGCCGAAGCAAACGGCCCGGGCAACACCCAGCGGAGCCgg                                                   >  2:305121/1‑143 (MQ=255)
                                                                                gTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCgg                                                   >  1:253863/1‑143 (MQ=255)
                                                                                   aCTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt                                                >  2:142024/1‑143 (MQ=255)
                                                                                    cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg                                               >  1:572531/1‑143 (MQ=255)
                                                                                              gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                          <  1:989355/58‑1 (MQ=255)
                                                                                              gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                          >  2:989355/1‑58 (MQ=255)
                                                                                              gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAt                                     <  2:24087/143‑1 (MQ=255)
                                                                                              gacggGGGAAGGATGTTAGGCGCCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAt                                     <  2:851498/143‑1 (MQ=255)
                                                                                                   gggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                >  1:935059/1‑143 (MQ=255)
                                                                                                    ggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGa                               <  1:194686/143‑1 (MQ=255)
                                                                                                         gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACg                          >  1:303831/1‑143 (MQ=255)
                                                                                                          aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGGCGAAGCACGGATGGTTGAAGACGa                         >  2:630518/1‑143 (MQ=255)
                                                                                                          aTGTCAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGa                         <  2:319029/143‑1 (MQ=255)
                                                                                                           tGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGaa                        >  1:834303/1‑143 (MQ=255)
                                                                                                           tGTTAGGCTGCGCGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCTGGATGAAGAACGGATGGTTGAAGACGaa                        <  2:410318/143‑1 (MQ=255)
                                                                                                              tAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAttt                     >  1:655129/1‑143 (MQ=255)
                                                                                                               aGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTa                    >  2:907730/1‑143 (MQ=255)
                                                                                                                ggcggcggGAAAGAAAGCGGGGCGTGAGCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                   <  1:614776/141‑1 (MQ=255)
                                                                                                                gTCGGCGGGAACGACAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                   <  1:796847/143‑1 (MQ=255)
                                                                                                                gTCGGCGGGAAAGAAGGCGGGGCGTGAACGCCCCGCGATTGATCATGACAGCATGAAGCTCCCCAGCCCGACGCCGTAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                   <  2:763219/143‑1 (MQ=255)
                                                                                                                        gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACgg           <  1:142024/143‑1 (MQ=255)
                                                                                                                         aaagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGACCGGATGGTTGAAGACGAATTTACCGATACggg          >  1:335886/1‑143 (MQ=255)
                                                                                                                             atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCACACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATCTACCGATACAGgtagt      >  1:312701/4‑143 (MQ=255)
                                                                                                                              tcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAGGACGAGTTTACCGATACGGgtagta     >  2:685726/3‑143 (MQ=255)
                                                                                                                              tcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta     >  1:247889/3‑143 (MQ=255)
                                                                                                                              tcgCGGGGCGTTCACGCCCCGCGATGGATGATTACCGCATGAAGCTACCCAGCAAGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGCGAACCCGGGAGGAAGAGCGGACGGTGGAAGGCGCAGATAGGGATACGGggagta     >  1:259784/3‑143 (MQ=255)
                                                                                                                                 cGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGCACCCGGGATGAAGAACGGACGGTCGAAGACGAATCTACCGATACGGGCAGTAccc  >  2:587699/1‑143 (MQ=255)
                                                                                                                                 cGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTAACCAGAACGAAGAAGAAGCAAACGGGTGGGTCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAccc  >  1:471252/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                     
CACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCC  >  NC_000913/4365315‑4365586

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