Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0790CT41.8% 69.1 / 39.1 55L134L (CTG→CTAthiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (4/28);  new base T (23/0);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 9.40e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTT  >  NC_000913/74934‑75228
                                                                                                                                                  |                                                                                                                                                           
aGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGtt                                                                                                                                                         <  1:145010/150‑1 (MQ=255)
    ttGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTtgtttgtt                                                                                                                                                     >  2:561958/1‑150 (MQ=255)
     tgtgGGGCGTCACCGCCATAGTCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGtt                                                                                                                                                     <  2:159878/149‑1 (MQ=255)
          ggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCa                                                                                                                                              <  2:295551/151‑1 (MQ=255)
           gCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTtttttttttttat                                                                                                                                             >  2:910416/1‑146 (MQ=255)
               cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTTATaaa                                                                                                                                          >  1:401907/1‑150 (MQ=255)
                 tCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTTATaaa                                                                                                                                          >  1:964781/1‑148 (MQ=255)
                   gCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTttttttttttttttaa                                                                                                                                          >  1:317745/1‑138 (MQ=255)
                             ttttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTtttttttttttat                                                                                                                                             >  1:160778/1‑128 (MQ=255)
                              tttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTTTTATAAAAGAAgggggaggagg                                                                                                                           >  1:1077197/1‑141 (MQ=255)
                                    tCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTATTTTGTTTTTTTTCTTAAAGAAAGCGaaagaaccatca                                                                                                                      >  1:732513/1‑139 (MQ=255)
                                     ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                                 <  2:486176/151‑1 (MQ=255)
                                     ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                                 <  2:756836/151‑1 (MQ=255)
                                        tAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTTTTTTTTTTTTTTATTAAAGAAGGCGCAGGAGCCCTAATTaaa                                                                                                               >  1:258866/1‑150 (MQ=255)
                                           cAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACg                                                                                                            <  2:984598/149‑1 (MQ=255)
                                            aaCAGCCCCAGCCCCGGTGTCCTGGTGCGCGGACCTTGATAACTCCCCCGCCAGTTTTGATCGCTCTCATCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTa                                                                                                         <  1:991239/151‑1 (MQ=255)
                                             aCAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAc                                                                                                        <  1:5743/151‑1 (MQ=255)
                                             aCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTATTTTGTTTTTTTCTTTAAAGCAAGCGCAGTAGTCACAATCaaaac                                                                                                             >  1:387107/1‑145 (MQ=255)
                                                gACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACg                                                                                                      <  2:28103/149‑1 (MQ=255)
                                                 aCCCAGCCCCGGTGTACTGGTGCGCGGATCCTAATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACg                                                                                                      <  1:721638/148‑1 (MQ=255)
                                                  cccAGCCGCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGa                                                                                                     <  2:413440/148‑1 (MQ=255)
                                                       ccccGGTGTACTGGTGCGCGGCTCCTGATAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCa                                                                                            <  2:153843/151‑1 (MQ=255)
                                                          cGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTCTTTGTTTTTGTTATAAAACAAAGCGGAAGAGGCAAAATTAAACGGGTACAAAGGGTCatt                                                                                          >  2:473885/1‑150 (MQ=255)
                                                           ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCatta                                                                                         <  2:840005/150‑1 (MQ=255)
                                                             tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCattatt                                                                                       <  2:1117658/150‑1 (MQ=255)
                                                                 cTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTc                                                                                      <  1:607903/147‑1 (MQ=255)
                                                                  tGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTACTTTGTTGTTTTCCTAAAAAAAGTGGAAAAAAGCACAAAAAAAAGGGACGAAAGTGTTTTTATTcc                                                                                     >  2:347363/1‑147 (MQ=255)
                                                                    gTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTTTAAAAGAAGGGGGAGTAGCCCTTAACCAAAGGGACGAAAAGGGCATTTTTTCAtccgc                                                                               >  2:600286/1‑151 (MQ=255)
                                                                       cgcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAATTTGTTTTTGTCTTAAATGAAAGAGAAGGAAGCCAAATAAAACAGATCGAAAGTGGCATTATTTCAagcgcccg                                                                            >  2:103228/1‑149 (MQ=255)
                                                                        gcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGGTTTTTTTTTTTCTAAAAAAAAGGGAAGGAACCCAAATTAAAACGGACGAAAAGGTCATTTTTTCAACCGCCcgg                                                                          >  2:1020552/1‑151 (MQ=255)
                                                                          ggtcctggATAAACCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGaa                                                                        <  1:231093/144‑1 (MQ=255)
                                                                            aTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTCATAAAACAAAGGGAAGGAGGCCTAATCAAAAGGGACGAAAAGGTCATTATTTCCACCCCCgggg                                                                         >  1:582065/1‑147 (MQ=255)
                                                                              catGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTGTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGt                                                                     <  1:159565/148‑1 (MQ=255)
                                                                               cTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTTTAAAAGAAAGCGAAGTAGGCATAATCAAAAGGGACGAAAGAGTCATTATTTCAGCCCCCGGGGACGt                                                                     >  1:951164/1‑149 (MQ=255)
                                                                                   tAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAAc                                                                 <  2:1110143/149‑1 (MQ=255)
                                                                                        cccccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATc                                                           <  1:1111594/148‑1 (MQ=255)
                                                                                           ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCt                                                       <  1:672609/151‑1 (MQ=255)
                                                                                            cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTg                                                      <  2:872419/151‑1 (MQ=255)
                                                                                            cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTTTGATAAAAGCAAGCGAAGGAGGCATAATTAAACAGGACGGAAGTGTCATTATTTCCGCCGCCGGGGACCGTTAAAggattcgc                                                         >  1:856764/1‑143 (MQ=255)
                                                                                            cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGCTTTTTGTTTTTTTTTTTTTCTAAACGAAGGCGAAGGAGCCATAATTAAAAGGGACGAAAATGGCCTTATTTTCACCGCCGGGGACGTTAAAGGCCTcccc                                                        >  1:175867/1‑147 (MQ=255)
                                                                                             gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTCTTTTTTCCTTTCTTTAACAAAGGCGCAGAAGACCAAAACAAACAGTGCGAAAGGGGCCGTTTTTCCGCCCCCCGGGACGGTTAACggatcc                                                           >  2:134477/1‑142 (MQ=255)
                                                                                                 gTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCATAAAAGAAAGGGAAGGAGCCATAATCAAAAGGGACGAAAGGGTCATTATTTCAGCCCCCGGGGACGTTAACGGCCTCCGCTGcc                                                    >  2:306704/1‑148 (MQ=255)
                                                                                                         cGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCttt                                          <  2:976357/150‑1 (MQ=255)
                                                                                                             ctcAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTgg                                       <  1:910292/149‑1 (MQ=255)
                                                                                                             ctcAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTCTTTTTATAAAAAAAGGAGAAGGAGGCAAAAATAAAAGGGAAGGAAACGGTATTTTTTCAGCCGCCGGGGAAGGTAAAGGGATCCCCCGCCACACCGCcttt                                         >  1:665580/1‑147 (MQ=255)
                                                                                                               cAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCaaa                                   <  2:530605/151‑1 (MQ=255)
                                                                                                               cAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCaa                                    <  1:1098632/150‑1 (MQ=255)
                                                                                                                  ccAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAg                                <  2:1060786/151‑1 (MQ=255)
                                                                                                                  ccAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTATTTTTTTTTTTTCTTAAAACAAGGGGAAGGAGCCCTAATCAAAAAGTACGGAAAGGGTATTTTTTCAGCCCCCCGGGACGTTAAAGGCATCCCCTGCCCCACCCCCTTTTGCcaa                                   >  1:369531/1‑148 (MQ=255)
                                                                                                                           tCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGGCCGGttt                        <  1:906052/150‑1 (MQ=255)
                                                                                                                                  ttGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggc                 <  2:240308/150‑1 (MQ=255)
                                                                                                                                      ggCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGgcggcg             >  1:992827/1‑150 (MQ=255)
                                                                                                                                          gggTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt          >  1:158537/1‑150 (MQ=255)
                                                                                                                                            gTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAc       >  2:415461/1‑151 (MQ=255)
                                                                                                                                               tttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAgg    >  2:601867/1‑151 (MQ=255)
                                                                                                                                                  cAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGgtt  >  1:649555/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                           
AGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTT  >  NC_000913/74934‑75228

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