Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,137,695 C→T 12.4% intergenic (+80/‑210) PP_2753 → / → PP_2754 ABC transporter ATP‑binding protein/OprD family outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,137,6950CT12.4% 137.0 / 10.6 57intergenic (+80/‑210)PP_2753/PP_2754ABC transporter ATP‑binding protein/OprD family outer membrane porin
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base T (3/4);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTT  >  NC_002947/3137549‑3137842
                                                                                                                                                  |                                                                                                                                                   
gggCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAgg                                                                                                                                                  <  2:1883997/150‑1 (MQ=255)
          gcCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCTGGGCCGCGttgaa                                                                                                                                        <  1:1725364/150‑4 (MQ=255)
             gAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCCGGCCGCGTTTTTTTc                                                                                                                                     <  1:1807821/150‑1 (MQ=255)
             gAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTc                                                                                                                                     <  1:1459102/150‑1 (MQ=255)
              aGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCTGGGCCGCGttgaatgga                                                                                                                                    <  2:1757691/150‑8 (MQ=255)
               gCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGCCGGCTGTGTAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGc                                                                                                                                   <  1:230447/150‑1 (MQ=255)
                         ggggAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGGTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGaaa                                                                                                                         <  2:114126/149‑1 (MQ=255)
                            gAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAgcg                                                                                                                      >  2:257425/1‑150 (MQ=255)
                            gAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAgcg                                                                                                                      >  1:1461678/1‑150 (MQ=255)
                             aaGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAgcgc                                                                                                                     <  1:1132118/150‑1 (MQ=255)
                              aGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCt                                                                                                                    >  1:1509004/1‑150 (MQ=255)
                                    tCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGgcg                                                                                                              >  2:658657/1‑150 (MQ=255)
                                      aGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGcc                                                                                                            >  1:1826189/1‑150 (MQ=255)
                                       gCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCAAGAGCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCg                                                                                                           >  1:959012/1‑150 (MQ=255)
                                         aTTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCTACTTGTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGc                                                                                                         >  1:99349/1‑150 (MQ=255)
                                         aTTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCAGCGGTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGc                                                                                                         >  1:99350/1‑150 (MQ=255)
                                          tttCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGcc                                                                                                        >  2:14494/1‑149 (MQ=255)
                                           ttCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCaa                                                                                                       >  1:396770/1‑150 (MQ=255)
                                                ccTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAAcagca                                                                                                  <  2:396770/150‑1 (MQ=255)
                                                  tGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGt                                                                                                <  1:658657/150‑1 (MQ=255)
                                                        gcAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCa                                                                                          <  2:128292/150‑1 (MQ=255)
                                                                           tGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGAATGGGTGTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACCGAAGTTGTGCTTATTTTTTCTCAGACGGAg                                                                       >  1:1368208/1‑150 (MQ=255)
                                                                                aCGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGcatca                                                                  >  2:1641847/1‑150 (MQ=255)
                                                                                   gCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAg                                                               >  1:1067333/1‑150 (MQ=255)
                                                                                      gtgGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAgtgt                                                            <  1:460046/150‑1 (MQ=255)
                                                                                        ggAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTg                                                          <  1:1372859/150‑1 (MQ=255)
                                                                                         gAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAgtgt                                                            <  2:1675527/147‑1 (MQ=255)
                                                                                         gAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAgtgt                                                            >  1:1675527/1‑147 (MQ=255)
                                                                                              gcgcAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAg                                                    >  2:1763587/1‑150 (MQ=255)
                                                                                                   gTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAgttgtt                                               <  1:1360127/150‑1 (MQ=255)
                                                                                                   gTCATACCGCCTTTTACTGAGTAGGTAATGTTCAACTTGAGGCGGCCCAGGCCGCGTGTTTTTCGCCCGCACGAAAGGGCGGACCGTTCTGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCACAGTGTGGAGTAAtttgtt                                               >  1:1383479/1‑150 (MQ=255)
                                                                                                      aTACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCATATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCaa                                            <  1:333275/150‑1 (MQ=255)
                                                                                                       tACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTAGGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCGTGATAGGGTTGAGGAAGTCGTTCAAt                                           >  1:124015/1‑150 (MQ=255)
                                                                                                        aCCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGACACAGCAGTTGTGCAGCATGCTTCGCATACTGAGCATCATATTTTTTAGGAAGTTGTTCAAtt                                          <  2:1454448/150‑1 (MQ=255)
                                                                                                          cGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAgttgt                                                <  1:1274573/142‑1 (MQ=255)
                                                                                                          cGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAgttgt                                                >  2:1274573/1‑142 (MQ=255)
                                                                                                          cGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGc                                        >  1:1543571/1‑150 (MQ=255)
                                                                                                                  aCTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGa                                    >  1:608092/1‑146 (MQ=255)
                                                                                                                  aCTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGa                                    <  2:608092/146‑1 (MQ=255)
                                                                                                                      aGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAAt                            <  2:1496989/150‑1 (MQ=255)
                                                                                                                         aGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGACTAACAATGGt                         <  2:1461678/150‑1 (MQ=255)
                                                                                                                          ggTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACCATGGtt                        >  1:613492/1‑150 (MQ=255)
                                                                                                                           gTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCGACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAgt                                                              <  2:1262922/111‑1 (MQ=255)
                                                                                                                           gTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAgt                                                              >  1:1262922/1‑111 (MQ=255)
                                                                                                                           gTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGttt                       >  2:891954/1‑150 (MQ=255)
                                                                                                                             gATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGttttt                     >  1:533363/1‑150 (MQ=255)
                                                                                                                                 ttCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCtgtg                                                                                                                 <  1:409618/54‑1 (MQ=255)
                                                                                                                                 ttCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCtgtg                                                                                                                 >  2:409618/1‑54 (MQ=255)
                                                                                                                                  tCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGtttttttttt                >  2:1595124/1‑149 (MQ=255)
                                                                                                                                   ccATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTttgt               <  2:1826189/150‑1 (MQ=255)
                                                                                                                                     aTTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCt                                                                                   <  2:1547192/80‑1 (MQ=255)
                                                                                                                                     aTTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCt                                                                                   >  1:1547192/1‑80 (MQ=255)
                                                                                                                                         aaCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTtgt                                                                                             <  1:483753/66‑1 (MQ=39)
                                                                                                                                         aaCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTtgt                                                                                             >  2:483753/1‑66 (MQ=39)
                                                                                                                                         aaCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTttgtt              <  1:244508/145‑1 (MQ=255)
                                                                                                                                         aaCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTttgtt              >  2:244508/1‑145 (MQ=255)
                                                                                                                                         aaCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCAt         >  1:1807508/1‑150 (MQ=255)
                                                                                                                                                cccAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGgttgtt  <  1:1763587/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTT  >  NC_002947/3137549‑3137842

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