Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 425,351 A→C 7.0% intergenic (‑39/+18) PP_0349 ← / ← PP_0350 membrane protein/ferrichrome‑iron receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947425,3510AC 7.0% 148.7 / 2.9 58intergenic (‑39/+18)PP_0349/PP_0350membrane protein/ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base C (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

ATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCCAG  >  NC_002947/425209‑425500
                                                                                                                                              |                                                                                                                                                     
aTAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcggc                                                                                                                                                <  2:150351/150‑1 (MQ=255)
     gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccgg                                                                                                                                           <  2:80727/150‑3 (MQ=255)
       gCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGct                                                                                                                                         >  2:418236/1‑150 (MQ=255)
               gTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGccc                                                                                                                                              >  2:257278/1‑137 (MQ=255)
               gTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGccc                                                                                                                                              <  1:257278/137‑1 (MQ=255)
                    tGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGCCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTg                                                                                                                            >  2:704165/1‑150 (MQ=255)
                      ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCTTGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTa                                                                                                                          >  1:1231747/1‑150 (MQ=255)
                      ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTa                                                                                                                          >  2:750543/1‑150 (MQ=255)
                      ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTa                                                                                                                          <  2:515516/150‑1 (MQ=255)
                       cAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACATGTaa                                                                                                                         <  2:342417/150‑1 (MQ=255)
                        aGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAg                                                                                                                        >  2:634065/1‑150 (MQ=255)
                          cAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTg                                                                                                                      <  2:984281/150‑1 (MQ=255)
                            aCAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGAAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGgc                                                                                                                    <  2:1248083/150‑1 (MQ=255)
                             cagcagGAACAGGGTACAGAGCAGGCTGCTCCAGGTGTGGATAAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGgcg                                                                                                                   <  2:88769/150‑1 (MQ=255)
                              agcagGAACAGGGTACAGACCAGGCTGCTCCAGGTTTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCCCCGCACCGATTGTTGTCGATAATTCGTCGTAACCCGTACGCGGTCCCAGCTCACTCAGAACTTGTAAGTGgcgc                                                                                                                  <  1:543138/150‑1 (MQ=255)
                              agcagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGgcgc                                                                                                                  >  1:574338/1‑150 (MQ=255)
                               gcagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTgg                                                                                                                     <  2:264052/146‑1 (MQ=255)
                               gcagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTGGTAAGTgg                                                                                                                     >  1:264052/1‑146 (MQ=255)
                                   gAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAg                                                                                                                                  >  2:697983/1‑129 (MQ=255)
                                   gAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAg                                                                                                                                  <  1:697983/129‑1 (MQ=255)
                                         ggTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACg                                                                                                       <  2:1179973/150‑1 (MQ=255)
                                           tACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCt                                                                                                     >  1:1157095/1‑150 (MQ=255)
                                                aCCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCt                                                                                                >  2:461544/1‑150 (MQ=255)
                                                  cAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCa                                                                                              <  2:1111680/150‑1 (MQ=255)
                                                       tgctCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGcc                                                                                         <  1:750543/150‑1 (MQ=255)
                                                               gtgtGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTg                                                                                                                            <  1:1263846/107‑1 (MQ=255)
                                                               gtgtGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTg                                                                                                                            >  2:1263846/1‑107 (MQ=255)
                                                                    gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                            >  2:637421/1‑150 (MQ=255)
                                                                    gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                            >  1:344051/1‑150 (MQ=255)
                                                                    gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                            >  2:711830/1‑150 (MQ=255)
                                                                    gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                            >  1:1218437/1‑150 (MQ=255)
                                                                     aTCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtagt                                                                           <  1:352803/150‑1 (MQ=255)
                                                                      tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAgtagta                                                                          <  2:1157095/150‑1 (MQ=255)
                                                                         aGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGa                                                                       <  1:637421/150‑1 (MQ=255)
                                                                           gACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccggg                                                                                                                                          >  1:377579/1‑78 (MQ=38)
                                                                                              ggCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATc                                                  <  1:794614/150‑1 (MQ=255)
                                                                                                 tctTCATCACGGTACAGATTGTTGTCGGTTGAACCACGGGGCCGAAAGGCGGACCCAGCTAAATTAGATCTTGTAAGTGGCGCTGTCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCtt                                               <  2:1231747/150‑1 (MQ=255)
                                                                                                   ttCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCa                                                                                                                                            >  2:759915/1‑55 (MQ=255)
                                                                                                   ttCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCa                                                                                                                                            <  1:759915/55‑1 (MQ=255)
                                                                                                   ttCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTc                                                                    >  1:587479/1‑127 (MQ=255)
                                                                                                   ttCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTc                                                                    <  2:587479/127‑1 (MQ=255)
                                                                                                      atcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTca                                          >  1:130131/1‑150 (MQ=255)
                                                                                                         aCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcagca                                       <  2:344051/150‑1 (MQ=255)
                                                                                                                  aTTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAgg                                                       <  2:1198755/125‑1 (MQ=255)
                                                                                                                  aTTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAgg                                                       >  1:1198755/1‑125 (MQ=255)
                                                                                                                      ttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACACGAACGGCGGTTTGCCGGGTCCAACAGGTTGGTGGCggg                          >  2:916451/1‑149 (MQ=255)
                                                                                                                        gtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGa                        <  1:634065/150‑1 (MQ=255)
                                                                                                                         tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCGCTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGt                                                                     <  1:1115750/104‑1 (MQ=255)
                                                                                                                         tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCGCTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGt                                                                     >  2:1115750/1‑104 (MQ=255)
                                                                                                                              tGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcga                  <  1:418236/150‑1 (MQ=255)
                                                                                                                                   gctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCgg                           <  1:319504/133‑1 (MQ=255)
                                                                                                                                   gctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCgg                           >  2:319504/4‑136 (MQ=255)
                                                                                                                                   cccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac             >  1:5913/1‑150 (MQ=255)
                                                                                                                                    ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac             <  2:1172086/149‑1 (MQ=255)
                                                                                                                                    ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac             >  1:1172086/1‑149 (MQ=255)
                                                                                                                                     cGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGAcgc           <  1:711830/150‑1 (MQ=255)
                                                                                                                                            gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtg                                  <  2:715390/120‑1 (MQ=255)
                                                                                                                                            gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtg                                  >  1:715390/1‑120 (MQ=255)
                                                                                                                                              aaGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTTCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCCAg  >  2:821693/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
ATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCCAG  >  NC_002947/425209‑425500

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