Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,063,726 T→G D17A (GAT→GCT)  cbpA ← DnaK co‑chaperone; curved DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,063,7260TG100.0% 206.8 / NA 57D17A (GAT→GCT) cbpADnaK co‑chaperone; curved DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/30);  total (27/30)

ATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAATTTCATATGCCAACACA  >  NC_000913/1063578‑1063872
                                                                                                                                                    |                                                                                                                                                  
aTCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGcg                                                                                                                                                  >  1:538795/1‑151 (MQ=255)
  cATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGcgtc                                                                                                                                                >  2:425550/1‑151 (MQ=255)
      cTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGcgtccgtc                                                                                                                                            >  2:638281/1‑151 (MQ=255)
        cAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGGTGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGcgtccgtc                                                                                                                                            <  1:284442/149‑1 (MQ=255)
         aGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGt                                                                                                                                         >  1:735234/1‑151 (MQ=255)
              gACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCAc                                                                                                                                    >  1:860887/1‑151 (MQ=255)
               aCGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGCCTTGATTGTCTTGAGAGCGTCCGGCGGTTTCACg                                                                                                                                   >  1:150319/1‑151 (MQ=255)
                cGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGc                                                                                                                                  >  1:800507/1‑151 (MQ=255)
                    gttCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCat                                                                                                                              <  1:373653/151‑1 (MQ=255)
                    gttCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCat                                                                                                                              <  1:429009/151‑1 (MQ=255)
                      tcatcaCTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCatg                                                                                                                             <  1:690311/150‑1 (MQ=255)
                       catcaCTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCatgat                                                                                                                           <  1:499498/151‑1 (MQ=255)
                            cTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCg                                                                                                                       >  2:122709/1‑150 (MQ=255)
                              tAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGta                                                                                                                     >  1:628476/1‑150 (MQ=255)
                              tAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGt                                                                                                                      <  2:859884/149‑1 (MQ=255)
                                   cTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTaataat                                                                                                                <  1:230932/150‑1 (MQ=255)
                                      cccAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCtt                                                                                                            <  2:544265/151‑1 (MQ=255)
                                       ccAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCttt                                                                                                           <  2:871377/151‑1 (MQ=255)
                                            cTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTc                                                                                                      >  2:17660/1‑151 (MQ=255)
                                              tCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTcc                                                                                                     <  1:177167/150‑1 (MQ=255)
                                                aGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATa                                                                                                  <  1:920027/151‑1 (MQ=255)
                                                aGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATa                                                                                                  <  2:566367/151‑1 (MQ=255)
                                                        ctTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTAtc                                                                                          >  2:44996/1‑151 (MQ=255)
                                                           tGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTcg                                                                                       >  2:190130/1‑151 (MQ=255)
                                                                  ggCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAAt                                                                                >  2:914864/1‑151 (MQ=255)
                                                                      tCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCaa                                                                             <  2:167673/150‑1 (MQ=255)
                                                                              cGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCTCGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACCCAAATTg                                                                    <  1:531127/151‑1 (MQ=255)
                                                                                    tttGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGa                                                               >  1:443426/1‑150 (MQ=255)
                                                                                     ttGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGaa                                                              <  1:250873/150‑1 (MQ=255)
                                                                                     ttGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGa                                                               <  1:488967/149‑1 (MQ=255)
                                                                                          gACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTg                                                        <  2:840248/151‑1 (MQ=255)
                                                                                            cATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTa                                                      <  1:123546/151‑1 (MQ=255)
                                                                                               cAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAg                                                    <  1:351705/150‑1 (MQ=255)
                                                                                                aGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAgg                                                   >  2:843213/1‑150 (MQ=255)
                                                                                                 ggATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGt                                                  >  1:483864/1‑150 (MQ=255)
                                                                                                 ggATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTa                                                 <  1:542145/151‑1 (MQ=255)
                                                                                                     ggTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACt                                              <  1:616008/150‑1 (MQ=255)
                                                                                                                gCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAgtgt                                   >  2:243159/1‑150 (MQ=255)
                                                                                                                gCAAGTCGACGATAAGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAgtgt                                   >  2:509384/1‑150 (MQ=255)
                                                                                                                  aaGTCGACGATAGGCGGTCTTGATTGTCTTGATAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCCCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAg                                 <  1:22662/150‑1 (MQ=255)
                                                                                                                   aGTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAgg                                >  1:501814/1‑150 (MQ=255)
                                                                                                                    gTCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGt                              <  1:445561/151‑1 (MQ=255)
                                                                                                                     tCGACGATAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGt                              >  2:375068/1‑150 (MQ=255)
                                                                                                                      cgacgaTAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTa                             <  2:398971/150‑1 (MQ=255)
                                                                                                                      cgacgaTAGGCGGGCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTa                             >  1:318340/1‑150 (MQ=255)
                                                                                                                         cgaTAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATcc                         >  1:598078/1‑151 (MQ=255)
                                                                                                                          gaTAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCt                        >  2:651568/1‑151 (MQ=255)
                                                                                                                            tAGGCGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTc                       >  1:712663/1‑150 (MQ=255)
                                                                                                                                cGGTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAAt                  <  1:53576/151‑1 (MQ=255)
                                                                                                                                 ggTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCTCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAAtt                 <  2:452561/151‑1 (MQ=255)
                                                                                                                                 ggTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAAtt                 >  1:699317/1‑151 (MQ=255)
                                                                                                                                 ggTCTTGATTGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAAtt                 <  2:452240/151‑1 (MQ=255)
                                                                                                                                  gTCTTGATTCTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAAttt                <  2:266240/151‑1 (MQ=255)
                                                                                                                                         ttGTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATACTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAATTTCATATGc         >  2:866271/1‑151 (MQ=255)
                                                                                                                                           gTCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAATTTCATATGCCa       <  1:434800/151‑1 (MQ=255)
                                                                                                                                            tCTTGAGAGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAATTTCATATGCCa       >  2:44020/1‑150 (MQ=255)
                                                                                                                                                gagaGCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAATTTCATATGCCAAcaca  <  1:336399/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGCCCATGATGGCGTAATAATCCTTTAATTCCATAGCGTTATCTCGCGTAAATCAACACAAATTGAAGGAACCCCTGTAAGGTAACTCCTATAAGTGTAGGGTAATCCTCAAAATTTCATATGCCAACACA  >  NC_000913/1063578‑1063872

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