Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 961,054 A→G intergenic (‑68/‑340) katA ← / → ssuB vegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome961,0540AG100.0% 247.7 / NA 64intergenic (‑68/‑340)katA/ssuBvegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/29);  total (35/29)

GGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAG  >  BS168_reference_genome/960918‑961188
                                                                                                                                        |                                                                                                                                      
ggCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGat                                                                                                                                      <  7:195143/139‑1 (MQ=255)
     tCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGAttatta                                                                                                                                 <  3:121690/139‑1 (MQ=255)
       aTTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATtattat                                                                                                                                >  8:162804/1‑138 (MQ=255)
       aTTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATtattat                                                                                                                                >  7:152253/1‑138 (MQ=255)
        ttGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATaa                                                                                                                              <  6:22838/139‑1 (MQ=255)
               ttGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGat                                                                                                                                      >  1:91870/1‑124 (MQ=255)
               ttGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGat                                                                                                                                      <  2:91870/124‑1 (MQ=255)
                  aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTAt                                                                                                                     <  7:381272/138‑1 (MQ=255)
                  aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTAt                                                                                                                     <  1:315206/138‑1 (MQ=255)
                       ctcCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTAt                                                                                                               >  5:218084/1‑139 (MQ=255)
                          cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTAtt                                                                                                              <  1:333639/137‑1 (MQ=255)
                          cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTaa                                                                                                            <  8:450481/139‑1 (MQ=255)
                          cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTaa                                                                                                            <  8:60944/139‑1 (MQ=255)
                           aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAAt                                                                                                           <  8:489423/139‑1 (MQ=255)
                           aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAAt                                                                                                           <  1:19152/139‑1 (MQ=255)
                           aaCCCGAGCGCCCCAGCTAGCTGTCAGATTATTTGCACGCATCTTACCACCTCAAGGAATTTATAAGAACATACTATCAGCTCTAAGTAAAAAGTCAATATTATTTTCTGATTATTATAAAATAGTTATCTCTATTAAt                                                                                                           <  6:369053/139‑1 (MQ=255)
                             ccGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCg                                                                                                         <  2:30770/139‑1 (MQ=255)
                                gAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCg                                                                                                         <  1:117401/136‑1 (MQ=255)
                                gAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCg                                                                                                         >  2:117401/1‑136 (MQ=255)
                                 aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATGACCTGTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAATGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGat                                                                                                       <  6:235706/137‑1 (MQ=255)
                                 aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGat                                                                                                                                      <  6:410786/106‑1 (MQ=255)
                                 aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGat                                                                                                                                      >  5:410786/1‑106 (MQ=255)
                                  gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGatata                                                                                                    <  5:33641/139‑1 (MQ=255)
                                  gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGatata                                                                                                    <  6:218084/139‑1 (MQ=255)
                                      cccAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTAtct                                                                                                                   <  1:415231/120‑1 (MQ=255)
                                      cccAGCTAGTTGTCAGTTTATTTGAACTAATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTAtct                                                                                                                   >  2:415231/1‑120 (MQ=255)
                                              gTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAt                                                                                        <  8:152253/139‑1 (MQ=255)
                                                    aGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAttgttg                                                                                   >  8:24427/1‑138 (MQ=255)
                                                    aGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAttgttg                                                                                   >  1:345112/1‑138 (MQ=255)
                                                    aGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATAATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAAttgttg                                                                                   >  6:269101/1‑138 (MQ=255)
                                                     gTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGaa                                                                                 <  8:1023/139‑1 (MQ=255)
                                                            tGAACTCAAGTTATCAACACCTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCt                                                                           <  6:372867/138‑1 (MQ=255)
                                                              aaCTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCttg                                                                         <  4:393949/138‑1 (MQ=255)
                                                              aaCTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCttg                                                                         <  8:335854/138‑1 (MQ=255)
                                                                cTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCttgttg                                                                      >  3:205581/1‑139 (MQ=255)
                                                                   aTGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTaa                                                                   >  4:389626/1‑139 (MQ=255)
                                                                    tGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAACCCTTGTTGTaa                                                                   >  8:63399/1‑138 (MQ=255)
                                                                        aTCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCt                                                              >  1:395437/1‑139 (MQ=255)
                                                                         tCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCt                                                              <  7:217755/138‑1 (MQ=255)
                                                                           aCCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCtttt                                                           >  3:428628/1‑139 (MQ=255)
                                                                           aCCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCtttt                                                           >  6:384295/1‑139 (MQ=255)
                                                                             ctctTGGAATTTATAAGAACATAATATCCTCTCTAAGAAAAAAGACAATATTATTTTATGATTATTATAAAATACGTAACTCTATTAATCGATATATAGTATCGAAATTGTTGAACACCCGTGTTGTACCCGCTGTTTa                                                         >  4:129934/1‑139 (MQ=255)
                                                                             ctctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTa                                                         <  8:327383/139‑1 (MQ=255)
                                                                             ctctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTa                                                         >  2:38954/1‑139 (MQ=255)
                                                                               ctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCTAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTAt                                                        >  2:145892/1‑138 (MQ=255)
                                                                               ctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTAt                                                        >  6:131397/1‑138 (MQ=255)
                                                                                 tGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGa                                                      >  4:419626/1‑138 (MQ=255)
                                                                                       ttATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTtaa                                               >  4:385695/1‑139 (MQ=255)
                                                                                          taAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTtaataa                                            >  3:358682/1‑139 (MQ=255)
                                                                                              aaCATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAAtttt                                        >  6:311166/1‑139 (MQ=255)
                                                                                                        catcTCTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTAGCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGatat                               >  5:39474/1‑138 (MQ=255)
                                                                                                          tctcTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGATAAACCTTGTTGTAACAGCTATTTATGACACTTAATAATTTTTCGAGAaataat                            >  1:227532/1‑139 (MQ=255)
                                                                                                          tctcTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat                            >  6:311536/1‑139 (MQ=255)
                                                                                                          tctcTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat                            >  6:275249/1‑139 (MQ=255)
                                                                                                          tctcTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat                            >  5:184216/1‑139 (MQ=255)
                                                                                                          tctcTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat                            >  8:444152/1‑139 (MQ=255)
                                                                                                          tctcTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat                            >  1:329137/1‑139 (MQ=255)
                                                                                                          tctcTAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACATTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat                            >  5:531277/1‑139 (MQ=255)
                                                                                                           ctctAAGTAAAAAGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAtt                           >  7:27976/1‑139 (MQ=255)
                                                                                                                     aaaGTCAATATTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTtctc                 >  4:16235/1‑139 (MQ=255)
                                                                                                                            atatTATTTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGc          <  5:275249/139‑1 (MQ=255)
                                                                                                                               ttattTTATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGa        <  5:269101/138‑1 (MQ=255)
                                                                                                                                  ttttATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACa    <  4:428628/139‑1 (MQ=255)
                                                                                                                                   tttATGATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAAc     >  7:31329/1‑137 (MQ=255)
                                                                                                                                     tatgATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAg  >  8:169449/1‑138 (MQ=255)
                                                                                                                                     tatgATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAg  <  1:145892/138‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAG  >  BS168_reference_genome/960918‑961188

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