Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,198,755 C→T 7.0% pseudogene (2282/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,198,7550CT 7.0% 181.3 / 3.8 57pseudogene (2282/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base T (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCC  >  NC_000913/2198607‑2198898
                                                                                                                                                    |                                                                                                                                               
cTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGc                                                                                                                                                 <  1:299042/149‑1 (MQ=255)
     cGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGa                                                                                                                                            >  3:345998/1‑149 (MQ=255)
      ggCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGaa                                                                                                                                           <  1:441741/149‑1 (MQ=255)
         ctgctgCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGTGGTGAACGc                                                                                                                                        >  2:690631/1‑149 (MQ=255)
          tgctgcCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGct                                                                                                                                       <  3:582290/149‑1 (MQ=255)
            ctgcctgcGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGctct                                                                                                                                     <  1:659955/149‑1 (MQ=255)
              gcctgcGGCGCTCGGCAAAGCCAGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGc                                                                                                                                   <  2:426520/149‑1 (MQ=255)
                ctgcGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCtt                                                                                                                                 >  2:686628/1‑149 (MQ=255)
                ctgcGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCtt                                                                                                                                 <  2:382340/149‑1 (MQ=255)
                       gcTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCg                                                                                                                          >  2:376026/1‑149 (MQ=255)
                           ggCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATcc                                                                                                                      <  1:686628/149‑1 (MQ=255)
                             cAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGTGGTGAACGCTCTGCTTGCGCTCGATCCCt                                                                                                                    <  1:690631/149‑1 (MQ=255)
                             cAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCt                                                                                                                    <  3:233329/149‑1 (MQ=255)
                             cAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCt                                                                                                                    <  1:286452/149‑1 (MQ=255)
                                                gATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACaaa                                                                                                 >  3:32064/1‑149 (MQ=255)
                                                gATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACaaa                                                                                                 >  4:639056/1‑149 (MQ=255)
                                                gATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACaaa                                                                                                 >  2:521380/1‑149 (MQ=255)
                                                      gCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATccc                                                                                           <  4:102853/149‑1 (MQ=255)
                                                        ccGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCa                                                                                         >  3:180553/1‑149 (MQ=255)
                                                            gcTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCActct                                                                                     >  2:84829/1‑149 (MQ=255)
                                                              tGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTc                                                                                   >  1:370083/1‑149 (MQ=255)
                                                                  ttGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGc                                                                               >  1:355193/1‑149 (MQ=255)
                                                                     cGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCtt                                                                            <  1:84829/149‑1 (MQ=255)
                                                                            ttACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAg                                                                     >  4:101488/1‑149 (MQ=255)
                                                                               ccGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCa                                                                  <  3:177396/149‑1 (MQ=255)
                                                                                      cGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGctct                                                           <  4:180553/149‑1 (MQ=255)
                                                                                       ggTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTg                                                          <  3:639056/149‑1 (MQ=255)
                                                                                            tcaGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATACTCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGAcgc                                                     >  3:183522/1‑149 (MQ=255)
                                                                                            tcaGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGAcgc                                                     >  3:233028/1‑149 (MQ=255)
                                                                                                cccttACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGATCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGccc                                                 >  4:60021/4‑149 (MQ=255)
                                                                                                   ttACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTAt                                                                        >  2:54038/1‑123 (MQ=255)
                                                                                                   ttACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTAt                                                                        <  1:54038/123‑1 (MQ=255)
                                                                                                      cTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTAtt                                           <  2:370083/149‑1 (MQ=255)
                                                                                                        gCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTaa                                         >  1:459067/1‑149 (MQ=255)
                                                                                                         cAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTaaa                                        >  4:124849/1‑149 (MQ=255)
                                                                                                          aaGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTaaaa                                       >  4:339225/1‑149 (MQ=255)
                                                                                                          aaGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTaaaa                                       >  4:591094/1‑149 (MQ=255)
                                                                                                           aGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGcc                                                  >  1:672413/1‑137 (MQ=255)
                                                                                                           aGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGcc                                                  <  2:672413/137‑1 (MQ=255)
                                                                                                               aTCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCaaa                                  >  3:262769/1‑149 (MQ=255)
                                                                                                                tCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAAt                                 <  4:345998/149‑1 (MQ=255)
                                                                                                                  gCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGc                               >  2:265428/1‑149 (MQ=255)
                                                                                                                    ccGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGAGCCCGGTATTAAAAGCAAATGccc                             >  2:602994/1‑149 (MQ=255)
                                                                                                                      gACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCaa                           >  4:402606/1‑149 (MQ=255)
                                                                                                                        cGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATc                         <  2:402637/149‑1 (MQ=255)
                                                                                                                          ttATAACCAGCCGGAAGTAACCGATGTGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTgg                                                         <  3:437654/115‑1 (MQ=255)
                                                                                                                          ttATAACCAGCCGGAAGTAACCGATGTGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTgg                                                         >  4:437654/1‑115 (MQ=255)
                                                                                                                          ttATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCAc                       <  3:591094/149‑1 (MQ=255)
                                                                                                                           tataACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACt                      <  4:32064/149‑1 (MQ=255)
                                                                                                                             taACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGc                    <  3:261857/149‑1 (MQ=255)
                                                                                                                             taACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGc                    <  1:376026/149‑1 (MQ=255)
                                                                                                                              aaCCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGcc                   <  4:599032/149‑1 (MQ=255)
                                                                                                                               aCCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGa                                                                                                                         >  3:658257/1‑46 (MQ=255)
                                                                                                                               aCCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGa                                                                                                                         <  4:658257/46‑1 (MQ=255)
                                                                                                                                   gCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGCCAGATa              <  1:602994/149‑1 (MQ=255)
                                                                                                                                           tAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCcctc      <  3:339225/149‑1 (MQ=255)
                                                                                                                                             aCCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCcctcct    >  2:465970/1‑149 (MQ=255)
                                                                                                                                               cGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCtcctcc  <  1:265428/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAACCAGCCGGAAGTAACCGATGCGGTGAACGCTCTGCTTGCGCTCGATCCCTTAGATAATCACCCGACAAAAATCCCCACTCTTCCGGCCTTTTATCAGCCATCGCTCTGGACGCGCCCGGTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCC  >  NC_000913/2198607‑2198898

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