Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,183,795 A→G T851A (ACC→GCC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,7950AG100.0% 210.9 / NA 57T851A (ACC→GCC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/30);  total (27/30)

GCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAA  >  NC_000913/4183651‑4183942
                                                                                                                                                |                                                                                                                                                   
gCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTg                                                                                                                                               <  2:1016769/151‑1 (MQ=255)
    tCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACAt                                                                                                                                           >  1:31385/1‑151 (MQ=255)
       tGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATccc                                                                                                                                        >  2:292580/1‑151 (MQ=255)
        gCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCg                                                                                                                                       >  1:781785/1‑151 (MQ=255)
         ccGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCg                                                                                                                                       >  1:691199/1‑150 (MQ=255)
           gTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCTCATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCGGGGTCCGGAAGAGATCGCCGCTGACATCCCGAAc                                                                                                                                    <  2:186679/151‑1 (MQ=255)
            tGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACg                                                                                                                                   >  2:279000/1‑151 (MQ=255)
              gAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTg                                                                                                                                 >  2:498236/1‑151 (MQ=255)
                     tACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGaa                                                                                                                           <  2:1061609/150‑1 (MQ=255)
                      aCAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCTCCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAg                                                                                                                          <  2:724150/150‑1 (MQ=255)
                       cAACTTCGAAGACTCCATCCGCAGATCCGAGCGTGTTGTTCCTGAACACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGGCACCAAGCGGGGTCGGGAAGAGATCGCCGCTGACATCCCGTTCGTGGGTGATGCt                                                                                                                        <  2:731908/151‑1 (MQ=255)
                        aaCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGCCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTg                                                                                                                       >  2:307669/1‑151 (MQ=255)
                            tCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTgcgc                                                                                                                    <  1:441476/150‑1 (MQ=255)
                             cGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGctc                                                                                                                  <  2:188451/151‑1 (MQ=255)
                             cGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGct                                                                                                                   <  2:125387/150‑1 (MQ=255)
                              gAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGctct                                                                                                                 <  1:116468/151‑1 (MQ=255)
                                     ccATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAc                                                                                                           <  1:213294/150‑1 (MQ=255)
                                     ccATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACt                                                                                                          <  2:345911/151‑1 (MQ=255)
                                       aTCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTgg                                                                                                        >  1:713536/1‑151 (MQ=255)
                                           tCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAt                                                                                                      >  2:1053844/1‑149 (MQ=255)
                                            cGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATg                                                                                                     <  1:745879/149‑1 (MQ=255)
                                                 ccGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCgg                                                                                              >  1:183385/1‑151 (MQ=255)
                                                 ccGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCg                                                                                               <  2:917425/150‑1 (MQ=255)
                                                    aGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTAt                                                                                           <  2:173390/151‑1 (MQ=255)
                                                            gttCAGGAAGACCGGTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGATGACATCCCGAACGTGGGTGAAGCTGCGCTCTTCAAACTGGATGAATCCGGTATCGTTTAc                                                                                    <  2:323703/150‑1 (MQ=255)
                                                               cAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAtt                                                                                 <  1:13420/150‑1 (MQ=255)
                                                                    aGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCg                                                                           <  1:6931/151‑1 (MQ=255)
                                                                     gACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCgg                                                                          <  1:843855/151‑1 (MQ=255)
                                                                      aCCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCg                                                                           >  2:564624/1‑149 (MQ=255)
                                                                       ccGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                         >  1:1024747/1‑150 (MQ=255)
                                                                       ccGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                         >  2:376318/1‑150 (MQ=255)
                                                                       ccGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGa                                                                         >  1:879532/1‑150 (MQ=255)
                                                                             caccacCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAcc                                                                  >  2:931295/1‑151 (MQ=255)
                                                                             caccacCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAcc                                                                  >  1:119281/1‑151 (MQ=255)
                                                                                  ccatccatATTCAGGAACTGGCGTGTGTGCCCCTTCACTCCAAGCTGGGTCCGAAAGAGATCGCCGCTGACATCCCGAACGTGCGTGGAGCTGCGCTCTCCATACGGGATGAATCCGGTATCGTTTACATTGGTGCGGTAGTGACCGGTg                                                              <  2:346663/150‑1 (MQ=255)
                                                                                   caTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGt                                                               <  2:667603/148‑1 (MQ=255)
                                                                                   caTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTCGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTgg                                                             >  2:959955/1‑150 (MQ=255)
                                                                                    aTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCg                                                           >  2:823517/1‑151 (MQ=255)
                                                                                    aTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCg                                                           <  2:651182/151‑1 (MQ=255)
                                                                                               ggAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggtt                                                <  2:644519/151‑1 (MQ=255)
                                                                                                  aCTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttggt                                             <  1:143220/151‑1 (MQ=255)
                                                                                                       cGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt                                        >  1:832288/1‑151 (MQ=255)
                                                                                                            gtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACg                                    <  1:1044815/150‑1 (MQ=255)
                                                                                                              gtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCg                                 >  2:1012607/1‑151 (MQ=255)
                                                                                                                cccGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGa                                >  1:609795/1‑150 (MQ=255)
                                                                                                                 ccGTGACACCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaa                               >  1:1068489/1‑150 (MQ=255)
                                                                                                                      acacCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGa                         <  2:707583/151‑1 (MQ=255)
                                                                                                                        acCAAGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGaaa                       <  2:526857/151‑1 (MQ=255)
                                                                                                                            aGCTGGGTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCa                   <  1:761243/151‑1 (MQ=255)
                                                                                                                               tGGGTCCGGAAGAGATCGCCGCTGACATCCCGGACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGc                 >  2:613679/1‑150 (MQ=255)
                                                                                                                                gggTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTg               <  2:613270/151‑1 (MQ=255)
                                                                                                                                gggTCCGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTg               <  2:869466/151‑1 (MQ=255)
                                                                                                                                    ccGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGTAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAccc           >  2:265738/1‑151 (MQ=255)
                                                                                                                                    ccGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAccc           >  1:160602/1‑151 (MQ=255)
                                                                                                                                     cGGAAGAGATCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc          >  2:898086/1‑151 (MQ=255)
                                                                                                                                           agaTCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGCGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTAAGCTGACCCCaga       >  2:999403/1‑148 (MQ=255)
                                                                                                                                             aTCGCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGaaa  <  1:220390/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAA  >  NC_000913/4183651‑4183942

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