Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,604 G→T 12.5% intergenic (+59/‑237) PP_2661 → / → PP_2662 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,048,6040GT12.5% 146.8 / 10.3 56intergenic (+59/‑237)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base T (3/4);  total (27/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 = 6.61e-01

GGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAG  >  NC_002947/3048461‑3048749
                                                                                                                                               |                                                                                                                                                   
ggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGccc                                                                                                                                                >  2:189075/1‑149 (MQ=255)
ggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCAAGGAGATTCCGCTCAGATAATGAAATTGGGGCTGCTTTGCGGCCCTTCGCGCGCTCGccc                                                                                                                                                >  5:75683/1‑149 (MQ=255)
     tgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCCGCGaa                                                                                                                                           <  1:290327/149‑1 (MQ=255)
     tgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGaa                                                                                                                                           <  2:138334/149‑1 (MQ=255)
      gCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGTGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                          >  8:287041/1‑149 (MQ=255)
       cAGGCACGGCGGGCCTGGCAGGAACCGCACCTTTATCTAATGCCCGCTTGCCAGCGTTCGCAGGCAATGGCTTCTTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCGGCCCTTCGCGGGCTCGCCCGCGAAg                                                                                                                                          <  6:75683/148‑1 (MQ=255)
        aGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAggg                                                                                                                                        <  1:189075/149‑1 (MQ=255)
           cACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGctg                                                                                                                                     <  8:294445/149‑3 (MQ=255)
               gCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAg                                                                                                                                 <  2:47928/149‑1 (MQ=255)
                   gcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                             <  2:37608/149‑1 (MQ=255)
                   gcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                             >  6:343111/1‑149 (MQ=255)
                    gcttgcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcggc                                                                                                                            >  6:152451/5‑149 (MQ=255)
                     ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCCGCGaa                                                                                                                                           >  1:353267/1‑133 (MQ=255)
                     ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCCGCGaa                                                                                                                                           <  2:353267/133‑1 (MQ=255)
                     ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGcc                                                                                                                           >  6:114969/1‑149 (MQ=255)
                        gcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGccc                                                                                                                          <  5:114969/147‑1 (MQ=255)
                             aaCCGCACCGTTATGAAATGCCCGCTTACCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttca                                                                                                                   <  1:305471/149‑1 (MQ=255)
                                   aCCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAAt                                                                                                             >  5:22108/1‑149 (MQ=255)
                                         aTGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGATGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                                                                                                      >  6:16361/1‑149 (MQ=255)
                                           gAAATGCCCCCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCTGGAGATTCCGCTCAGATAAAGTAATTGGGGCTGCTTTGCACCCCTACGCTGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGATCATTAAATCGAATTg                                                                                                     >  8:62260/1‑148 (MQ=255)
                                              aTGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAc                                                                                                 <  5:377523/149‑1 (MQ=255)
                                                 cccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCg                                                                                              >  8:379985/1‑149 (MQ=255)
                                                   cGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTTGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGCCCCCGGTTCATTCAATGGGGATAGATACCCGc                                                                                             >  1:139635/1‑149 (MQ=255)
                                                         ccAGCGTTCGCAGGAAGTGGCGTCCTGAAGCCAGGAGATTCCGCGCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCAc                                                                                      <  8:67081/149‑1 (MQ=255)
                                                               ttCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcctgc                                                                                <  7:287041/149‑1 (MQ=255)
                                                                   cAGGCACTGGCTTCATTAAGCCAGGAGATTCCGCTCAGATAATGTATTTGGGGCTGCTTTGCAGCCCTTCGCGGGGGAGCCCGCGAAGGGCCTCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGa                                                                            <  2:139635/149‑1 (MQ=255)
                                                                    aGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAg                                                                           <  5:343111/149‑1 (MQ=255)
                                                                      gCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGTGCGAGACCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCTTATTGATACACGCCTGCCACTCCTGCCTGAGc                                                                          >  2:122853/1‑148 (MQ=255)
                                                                        aCTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCa                                                                        >  7:167748/1‑148 (MQ=255)
                                                                        aCTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCa                                                                        >  2:305047/1‑148 (MQ=255)
                                                                            gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATaa                                                                     <  1:179824/147‑1 (MQ=255)
                                                                            gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGACCCATAAGt                                                                   >  8:71967/1‑149 (MQ=255)
                                                                                   gAAGCCAGGAGATTCCGCTCAGATAAGGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGCAGGGCCGCAAGGCGGCCCCGGTTCAGTCAATCGGATTGATACCCGCCTGCCACTCCTGCCTGCGCCATAAGTCTCGAAc                                                            >  6:154285/1‑149 (MQ=255)
                                                                                     aGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCa                                                          <  4:182615/149‑1 (MQ=255)
                                                                                         aGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAccc                                                                                               >  3:263879/1‑108 (MQ=255)
                                                                                          ggAGATTCCTCTCAGTTATTGTAATTGGGGCTGCTTTGCAGCCCTGCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtc                                                                  <  2:289956/136‑1 (MQ=255)
                                                                                          ggAGATTCCCCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTAATTCAATCGTAATGATACACGCCTGCCACTCCTGCCTGAGCCATAACtc                                                                  >  1:289956/1‑136 (MQ=39)
                                                                                                 ccGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCa                                              <  8:167748/148‑1 (MQ=255)
                                                                                                  cGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAg                                             <  7:379985/148‑1 (MQ=255)
                                                                                                     tCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTAc                                         >  4:207834/1‑149 (MQ=255)
                                                                                                       ataTAATTTATTTGGGGCTGCTTTGCAGCCCTTCGGGGGCGAGCCCGCTAAGGGCCGCAAGGCGGCCCGGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTAc                                         <  6:22108/145‑1 (MQ=255)
                                                                                                             tGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGAACCGGGAcc                                  >  4:246093/1‑148 (MQ=255)
                                                                                                               tAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCaaa                               >  5:237198/1‑149 (MQ=255)
                                                                                                               tAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGGCTCGAACCAGGACCCAAGCAGGTACCGGGAccc                                 >  4:313924/1‑146 (MQ=255)
                                                                                                                 aTTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATc                             >  6:173076/1‑149 (MQ=255)
                                                                                                                         tgctTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGcc                     >  5:318860/1‑149 (MQ=255)
                                                                                                                                    ccTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCg          <  6:318860/149‑1 (MQ=255)
                                                                                                                                       tCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATc                             <  7:194014/127‑1 (MQ=255)
                                                                                                                                       tCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATc                             >  8:194014/1‑127 (MQ=255)
                                                                                                                                        cgcgGGCGAGCCCGCTAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGcc        <  7:62260/147‑1 (MQ=255)
                                                                                                                                        cgcgGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGccc       <  3:207834/148‑1 (MQ=255)
                                                                                                                                         gcgGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGcccc      >  7:227600/1‑148 (MQ=255)
                                                                                                                                          cgGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaa    <  7:71967/149‑1 (MQ=255)
                                                                                                                                           gggCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaaa   <  5:152451/149‑1 (MQ=255)
                                                                                                                                             gCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAg  <  1:305047/149‑1 (MQ=255)
                                                                                                                                             gCGAGCCCGCGAAGGGACTCAAGGCGGCCCAGGTTCATTCAATCGTATTTATACCCTCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaaa   <  5:16361/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAG  >  NC_002947/3048461‑3048749

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: