Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,144 T→G 8.0% intergenic (+61/+45) PP_1587 → / ← dapC NhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,780,1440TG 8.0% 131.3 / 2.5 50intergenic (+61/+45)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (33/11);  new base G (2/2);  total (37/13)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

TCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCC  >  NC_002947/1779999‑1780293
                                                                                                                                                 |                                                                                                                                                      
tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgc                                                                                                                                                    <  1:2522895/150‑1 (MQ=255)
 cAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                                                                   <  1:3359729/150‑1 (MQ=255)
   gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                                                                 <  1:4681936/150‑1 (MQ=255)
   gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGTAATGCGCTGCGCATGAAATCAACGCt                                                                                                                                                 >  2:4312087/1‑150 (MQ=255)
    cACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTg                                                                                                                                                <  1:4033824/150‑1 (MQ=255)
       cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAc                                                                                                                                             <  1:930868/150‑1 (MQ=255)
       cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACACGTAAGTCGCGGCGCTTGAAATAAGCGCTGTAc                                                                                                                                             >  1:3686333/1‑150 (MQ=255)
        tggtggGGGATGTGGTGGTGCGGCTAGTACTGGGTGAGCTGGTTTTGAGTGAGGCGATCCTGGGTTCAGTCATGTAGACCCTTTCAAGGCCCTTGGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGATATCAGCGCTGTAcc                                                                                                                                            <  2:203192/150‑1 (MQ=255)
           tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                                        <  1:553067/150‑1 (MQ=255)
           tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                                        <  1:1557468/150‑1 (MQ=255)
             ggggTTGTGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGTGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                                      <  1:3883468/150‑1 (MQ=255)
             ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCCGGTACTGCGCGGCGGGTGGAATCAGCGCGGTGACTGTgg                                                                                                                                      >  1:2401824/1‑150 (MQ=255)
             ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACCGGGACAGCGCTGCGCTTGAAATCACCGCGGTAACGGTgg                                                                                                                                      >  1:2401825/1‑150 (MQ=255)
             ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                                                      <  1:2035980/150‑1 (MQ=21)
                 aTGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGc                                                                                                                                  <  1:2669386/150‑1 (MQ=255)
                   gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCgg                                                                                                                                <  1:4352409/150‑1 (MQ=255)
                   gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCgg                                                                                                                                <  2:441015/150‑1 (MQ=21)
                     ggtggtGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAACAGTGGGAGCGGGc                                                                                                                              >  1:3313696/1‑150 (MQ=255)
                     ggtggtGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCAGTACCAGTGGGAGCGGGc                                                                                                                              >  1:3313594/1‑150 (MQ=255)
                                   tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTACGCGGCCGCGCCCGCTCACAAAGGTACAGCGCAGCGCTTGAAATCAGAGCAGTACCAGTGGGAGCGGGCGTGGCCGCGAAggg                                                                                                                >  1:2976216/1‑150 (MQ=255)
                                        ggTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCGCCCACCAGTACAGCGCGGGGCGTGAAAAAAACCCGGTAACGGTGGGGGCGGGGGGGGCCGCGGAGGGGGTgt                                                                                                           >  1:1034453/1‑149 (MQ=255)
                                              cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACAGGGACTGGGCAGCGCGTGTAATAAGCCCGGTAACGGTTGGGGGGGGGCGGGCCCCCGAGGGGGTGTTCAGCg                                                                                                     >  1:2357828/1‑150 (MQ=255)
                                                 gATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGGACAGCGCTGCGCGTGAAAAAAGCGCAGTTCCAGTGGGGGCGGGCGTGCCCGCGAACGGGTTGGTCAGCGGtt                                                                                                  >  1:1154511/1‑150 (MQ=255)
                                                         tgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGa                                                                                          >  2:3511887/1‑150 (MQ=255)
                                                             gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTTAAATCAGCGCTTTAACTTTGGGAGCGGGCGGGCCCGCGAACGGGTTGAACACCGGTTTTGCAGGAAGgc                                                                                      >  2:3937483/1‑150 (MQ=38)
                                                                                aGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACAGGTAATGCCCAGCCCGTGCAATAAGCGCTGTACCAGTGGGAGCGGCCGTGCCCCCGAAGGGGTTGAACACCGGATTAGAAGGAAGGCCCGAAAGCGCAacgcctca                                                                   >  2:329338/1‑147 (MQ=255)
                                                                                          aTCAAGGCCCTTCGCGGGCGCGCCCGCTCCCAAAGGTACTGCGCTGCGCGTGAAATCCCCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGGAGGGGGTGATCAGCCGTGTTGCGGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCAtt                                                         >  2:4978599/1‑150 (MQ=37)
                                                                                                         gggCGCGCCCGCTCCCACACGTACTGCCCCGCGCGTGAAATACCCCCCGTACCAGGGGGAGCGGGCGTGCCCCCGCAGGGGGTGATCAGCCGGTTTTCAGGGAGGCGCGAAAGCGCACCGCCCCCTCGATGCATTCCGCCCCCGGTGCaa                                          >  2:1403564/1‑150 (MQ=255)
                                                                                                             gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg                                      >  2:3522682/1‑150 (MQ=255)
                                                                                                             acgcCCGCTCCCACCGGTACAGCGCTGATTTCAAGCGCAGCCCAGTACCTGTGGGAGCCGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCCGGAAGGCGCGAATGCGCACCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg                                      >  2:4473450/2‑150 (MQ=18)
                                                                                                             acgcCCGCTCCCACACGTACAGAGCTGATTTCAATCGCAGCGCAGTAACTGTGGGAGCGGGCGTGCCCGCGAATGGGTTGAACAGCGGTTTTTCAGGAATGCCCCAATGCGCTCCGCCGCTTCGATGCAATCGGCCACCGGTGCAACCaa                                      >  2:4473106/2‑149 (MQ=19)
                                                                                                                 ccGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                                  >  2:301536/1‑150 (MQ=255)
                                                                                                                       ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                            >  1:2106983/1‑150 (MQ=255)
                                                                                                                        cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                           >  2:874701/1‑150 (MQ=255)
                                                                                                                        cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                           >  2:5035270/1‑150 (MQ=255)
                                                                                                                         acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                          >  1:1121758/1‑150 (MQ=255)
                                                                                                                         acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                          >  1:3826808/1‑150 (MQ=255)
                                                                                                                         acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                          >  1:1317749/1‑150 (MQ=255)
                                                                                                                         acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                          >  1:2996155/1‑150 (MQ=255)
                                                                                                                           aGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgc                        >  1:1126461/1‑150 (MQ=255)
                                                                                                                              tACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGcc                     >  2:3754009/1‑150 (MQ=255)
                                                                                                                                 tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg                  >  2:4001991/1‑150 (MQ=255)
                                                                                                                                 tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg                  >  2:2409102/1‑150 (MQ=255)
                                                                                                                                 tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGAGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg                  >  2:4520175/1‑150 (MQ=255)
                                                                                                                                   cgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCa                >  1:3960333/1‑150 (MQ=255)
                                                                                                                                      tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc             >  1:2713673/1‑150 (MQ=255)
                                                                                                                                           ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt        >  1:3462977/1‑150 (MQ=255)
                                                                                                                                                aTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAccc   >  1:1941426/1‑150 (MQ=255)
                                                                                                                                                aTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAccc   >  1:1941249/1‑150 (MQ=255)
                                                                                                                                                aTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAccc   >  2:4586828/1‑150 (MQ=255)
                                                                                                                                                 tCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcccc  >  1:2963416/1‑150 (MQ=255)
                                                                                                                                                 tCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcccc  >  1:2963355/1‑150 (MQ=255)
                                                                                                                                                 tCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcccc  >  1:2963177/1‑150 (MQ=255)
                                                                                                                                                 tCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcccc  >  2:3122192/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCC  >  NC_002947/1779999‑1780293

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