Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 425,348 A→C 9.8% 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,3480AC 9.8% 155.1 / 6.6 58intergenic (‑39/+18)PP_0349/PP_0350membrane protein/ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base C (3/3);  total (29/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 = 9.72e-01

AGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCC  >  NC_002947/425204‑425491
                                                                                                                                                |                                                                                                                                               
aGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCTATTGTTGTCGGTTGGAACCCGGGGCCGCAAggc                                                                                                                                             <  1:237580/149‑1 (MQ=255)
     aGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGcc                                                                                                                                         <  6:304905/148‑1 (MQ=255)
       gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGccc                                                                                                                                        >  7:257396/1‑147 (MQ=255)
        ggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGATGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccgg                                                                                                                                     <  8:44690/149‑3 (MQ=255)
          cAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGct                                                                                                                                   >  8:271571/1‑149 (MQ=255)
               cGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCCGctgcctc                                                                                                                              >  2:271136/1‑144 (MQ=255)
                ggTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctcact                                                                                                                               >  4:273982/1‑147 (MQ=255)
                             aaCAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTg                                                                                                                >  3:284368/1‑149 (MQ=255)
                                 gcagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGaa                                                                                                                          <  4:44045/135‑1 (MQ=255)
                                 gcagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGaa                                                                                                                          >  3:44045/1‑135 (MQ=255)
                                    ggAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTgg                                                                                                         >  4:248346/1‑149 (MQ=255)
                                                aGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgc                                                                                             <  3:273982/149‑1 (MQ=255)
                                                aGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgc                                                                                             >  2:236089/1‑149 (MQ=255)
                                                 gACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgcg                                                                                            <  1:46541/149‑1 (MQ=255)
                                                  aCCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgcgc                                                                                           <  8:203868/149‑1 (MQ=255)
                                                  aCCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTgcgc                                                                                           >  7:108277/1‑149 (MQ=255)
                                                    cAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACaa                                                                                                   >  7:118699/1‑139 (MQ=255)
                                                    cAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACaa                                                                                                   <  8:118699/139‑1 (MQ=255)
                                                       gctgctCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGCCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCCCCACGCGGCCCCCGCACACACCGAACAAGTAAGTTGCGCTGGCGGCAACCCTGCGCGCaac                                                                                      >  6:240038/1‑147 (MQ=255)
                                                           ctCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGcc                                                                                   <  4:284368/148‑1 (MQ=255)
                                                           ctCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCg                                                                                  <  7:271571/149‑1 (MQ=255)
                                                            tCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCt                                                                                               >  3:307594/1‑135 (MQ=255)
                                                            tCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCt                                                                                               <  4:307594/135‑1 (MQ=255)
                                                                 gtgtGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGtagtag                                                                            >  8:19138/1‑149 (MQ=255)
                                                                 gtgtGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGtagtag                                                                            >  7:193136/1‑149 (MQ=255)
                                                                      gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGc                                                                                                        >  5:116663/1‑116 (MQ=255)
                                                                      gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGc                                                                                                        <  6:116663/116‑1 (MQ=255)
                                                                            ggACGAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGa                                                                 <  5:240038/149‑1 (MQ=255)
                                                                            ggACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGa                                                                 <  8:193136/149‑1 (MQ=255)
                                                                              aCCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGACCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGt                                                               >  6:161852/1‑149 (MQ=255)
                                                                                cAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTca                                                             <  6:380021/149‑1 (MQ=255)
                                                                                      cggATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTa                                                       <  7:19138/149‑1 (MQ=255)
                                                                                             tggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCCCCACACTCAGAAATTGTAATTGGCGCTGGCGACAACGCTGCGCTCAAAGCCGTAGTAGGAGTAGTAGGAGTCACAAGTAGAGAgg                                                 >  4:39472/1‑148 (MQ=255)
                                                                                             tggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAgg                                                 >  1:331746/1‑148 (MQ=255)
                                                                                                  ctctTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGACCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATcc                                           <  5:161852/149‑1 (MQ=255)
                                                                                                   tctTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAAAGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCt                                          >  1:285431/1‑149 (MQ=255)
                                                                                                       catcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTTGCGACAACGCTGCGGACATCGCCGTAGTAGCAGTAGTAGGAGTGACAAGTAGAGAGGTTATCCATGt                                       >  8:33059/1‑148 (MQ=255)
                                                                                                         tcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTca                                    <  1:236089/149‑1 (MQ=255)
                                                                                                                 ccGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtg                            <  3:39472/149‑1 (MQ=255)
                                                                                                                   gATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagagag                                                  <  1:237071/125‑1 (MQ=255)
                                                                                                                   gATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagagag                                                  >  2:237071/1‑125 (MQ=255)
                                                                                                                   gATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtggt                          >  8:85469/1‑149 (MQ=255)
                                                                                                                   gATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtggt                          <  3:248346/149‑1 (MQ=255)
                                                                                                                        ttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCg                      <  7:85469/148‑1 (MQ=255)
                                                                                                                          gtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTCGCGACAACGCTGCGCTCATCGCCGTAGTAGAAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGt                    >  8:187592/1‑148 (MQ=255)
                                                                                                                                   gagctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTg                                        >  5:121221/6‑119 (MQ=255)
                                                                                                                                   gagctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTg                                        <  6:121221/114‑1 (MQ=255)
                                                                                                                                   aaCCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgacc          <  2:331746/149‑1 (MQ=255)
                                                                                                                                      ctggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGt                              >  3:31307/3‑126 (MQ=255)
                                                                                                                                      ctggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGt                              <  4:31307/124‑1 (MQ=255)
                                                                                                                                      ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagag                                                    >  7:290395/1‑104 (MQ=255)
                                                                                                                                      ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagag                                                    <  8:290395/104‑1 (MQ=255)
                                                                                                                                      ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac       <  8:257396/149‑1 (MQ=255)
                                                                                                                                           gCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGcc  >  7:331100/1‑149 (MQ=255)
                                                                                                                                            ccGCAAGGCGGCCCCAGCTCACTCAGAACTTTTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcga            <  1:124432/138‑1 (MQ=255)
                                                                                                                                            ccGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcga            >  2:124432/1‑138 (MQ=255)
                                                                                                                                              gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac       <  6:323712/141‑1 (MQ=255)
                                                                                                                                              gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac       >  5:323712/1‑141 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AGATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCC  >  NC_002947/425204‑425491

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