Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,104,809 C→T 7.0% I345I (ATC→ATT yhhS → carbohydrate efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,104,8090CT 7.0% 165.1 / 2.6 57I345I (ATC→ATTyhhScarbohydrate efflux transporter
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/2);  total (27/30)
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

GCACCCAACGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCCTCTATTTGTACCGCCAAGGGCGTAGGGACAG  >  NC_002947/5104663‑5104955
                                                                                                                                                   |                                                                                                                                                  
gcACCCAACGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGa                                                                                                                                                  >  4:130000/1‑149 (MQ=255)
gcACCCAACGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCAGTGGGCGTGGAGGCAGTGAACCAAGTCTCGGCGGCTAACCGGGGGGAGGCGGTGGGGGCGTATTCGCTGTTCATCGa                                                                                                                                                  >  8:111190/1‑149 (MQ=255)
 cACCCAACGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGAc                                                                                                                                                 <  6:118875/149‑1 (MQ=255)
 cACCCAACGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGAc                                                                                                                                                 >  5:138697/1‑149 (MQ=255)
      aaCGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGt                                                                                                                                             >  1:361197/1‑148 (MQ=255)
       aCGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCa                                                                                                                                           <  8:3495/149‑1 (MQ=255)
         gCCGAACTGGCCTTGGCGCGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGATCCAAGTGTCGGCGGCTAACCGCGGGGCGGCTGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACt                                                                                                                                         >  5:302677/1‑149 (MQ=255)
          ccGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTg                                                                                                                                        <  3:130000/149‑1 (MQ=255)
             aaCTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTgggg                                                                                                                                     <  8:342869/149‑1 (MQ=255)
             aaCTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTAACTgggg                                                                                                                                     >  4:30932/1‑149 (MQ=255)
                cggCGTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGt                                                                                                                                   <  1:181896/147‑1 (MQ=255)
                 ggCCTTGGCGGGGGCTGCATTGAGCGGGGTTGGTTTCTCGCTGGTGTTTCCGGCCTTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGGATTCGCTGTTCATCGCCTTGTCACTGGGGGTGa                                                                                                                                 <  6:312098/149‑1 (MQ=255)
                       ggCGTGGGCTGCATTTAGAGGGTTGGGTTTCTCGCTGGTGTTTGCGGCATTGGGCGGGGATGCATTGAAGCAAGTGTCGGCGGCTAAGCGCGGTGCGGGGGTGGGGTCGTATTCGCTGTTCATCGACTTGTCACTGTGGGTGACCtgg                                                                                                                            <  3:20258/148‑1 (MQ=255)
                       ggCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATTGACTTGTCACTGGGGGTGACCGGGc                                                                                                                           >  4:96354/1‑149 (MQ=255)
                                  cATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTcggt                                                                                                                <  7:17457/149‑1 (MQ=255)
                                            gtttggtttCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTggcgg                                                                                                      >  3:227110/1‑149 (MQ=255)
                                             tttggtttCTCGCTGGTGTTTCCGGCATGGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCCTTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTggcgg                                                                                                      <  1:292229/149‑1 (MQ=255)
                                                ggtttCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGTCTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGcggcgg                                                                                                   <  1:269436/148‑1 (MQ=255)
                                                    tctcGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTc                                                                                              >  1:178176/1‑149 (MQ=255)
                                                    tctcGCTGGTGTTTCAGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCCGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACAGCGGGTGACAGGGCCATTGGTCGGAGCGGTCGCGGCGGGGTTc                                                                                              >  2:115280/1‑149 (MQ=255)
                                                     ctcGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATTGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTc                                                                                              <  3:96354/148‑1 (MQ=255)
                                                            gtgtTTCCGGCATTGGGCGGGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTGATCGTCTTGTCACTGGGGGTGACCGGGCCTTTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGc                                                                                      <  7:111190/149‑1 (MQ=255)
                                                              gtTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGc                                                                                           >  8:100410/1‑142 (MQ=255)
                                                              gtTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGc                                                                                           <  7:100410/142‑1 (MQ=255)
                                                               tttCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTc                                                                                   >  4:366812/1‑149 (MQ=255)
                                                                  ccGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGa                                                                                 <  2:178176/148‑1 (MQ=255)
                                                                  ccGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGAt                                                                                >  7:117349/1‑149 (MQ=255)
                                                                   cGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTcggt                                                                                                                >  8:1611/1‑116 (MQ=255)
                                                                   cGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTcggt                                                                                                                <  7:1611/116‑1 (MQ=255)
                                                                   cGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATg                                                                               >  4:53557/1‑149 (MQ=255)
                                                                      ctttGGGCGTGGAGGCTGTGAACCAAGTGTCGTCGGCTAACCGCGGGGGGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGTGGTGACCGGGCGATTGGTCGGTGCGGTGGGGGCGGGGTTCGGCTTTGCTTCGATtttt                                                                            <  6:151164/147‑1 (MQ=255)
                                                                           gggctGGATGCAGTTAACCAAGTGGCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGt                                                                        <  6:138697/144‑1 (MQ=255)
                                                                           ggCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGtt                                                                       <  4:227110/149‑1 (MQ=255)
                                                                            gCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTc                                                                                              <  5:210849/125‑1 (MQ=255)
                                                                            gCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGCCCATTGGTCGGTGCGGTGGCGGCGGGGTTc                                                                                              >  6:210849/1‑125 (MQ=255)
                                                                                        tGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCg                                                          >  1:261816/1‑149 (MQ=255)
                                                                                               aGTGTCGTCGGCTAACCGGGGGGCGGCGTTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGt                                                   <  8:117349/149‑1 (MQ=255)
                                                                                                       cggcTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTg                                           <  2:261816/149‑1 (MQ=255)
                                                                                                         gcTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGgt                                         <  3:366812/149‑1 (MQ=255)
                                                                                                          cTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGgtg                                        >  6:322466/1‑149 (MQ=255)
                                                                                                            aaCCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCa                                                       <  3:49891/132‑1 (MQ=255)
                                                                                                            aaCCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCa                                                       >  4:49891/1‑132 (MQ=255)
                                                                                                              ccGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGtt                                      >  6:246520/1‑147 (MQ=255)
                                                                                                                   gggCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCctc                               >  8:258110/1‑149 (MQ=255)
                                                                                                                        gcggTGGGGGCGTATTCGCTGTTCATTGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCgg                                               >  4:78616/1‑128 (MQ=255)
                                                                                                                        gcggTGGGGGCGTATTCGCTGTTCATTGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCgg                                               <  3:78616/128‑1 (MQ=255)
                                                                                                                               gggCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTcg                                                                     <  6:42701/99‑1 (MQ=255)
                                                                                                                                 gCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCgcgg                                                                  <  3:130032/100‑1 (MQ=255)
                                                                                                                                 gCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCgcgg                                                                  >  4:130032/1‑100 (MQ=255)
                                                                                                                                     tATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCCTCTATTTGTACCGCCAAgg              >  6:28258/1‑149 (MQ=255)
                                                                                                                                         cGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCCTCTATTTGTACCGCCaa                <  5:185500/143‑1 (MQ=255)
                                                                                                                                         cGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCCTCTATTTGTACCGCCaa                >  6:185500/1‑143 (MQ=255)
                                                                                                                                            tGTTCATCGACTTGTCACTGTGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCCTCTATTTGTACCGCCAAGGGCGTagg       <  3:30932/149‑1 (MQ=255)
                                                                                                                                                cATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCCTCTATTTGTACCGCCAAGGGCGTaggga     >  8:218342/1‑147 (MQ=255)
                                                                                                                                                 aTCGAGTTGTCTCTGGTGGTGTCCGGTCCATTGGTGGGTGCGGTGCCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGTGTTGGTGTTGAGCCTCTATTTGTACCGCCAAGGGCGTagggacag  <  6:302677/149‑1 (MQ=255)
                                                                                                                                                 aTCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGa                                    >  5:106122/1‑115 (MQ=255)
                                                                                                                                                 aTCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGa                                    <  6:106122/115‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCACCCAACGCCGAACTGGCCTTGGCGGGGGCTGCATTGAGCGGGTTTGGTTTCTCGCTGGTGTTTCCGGCATTGGGCGTGGAGGCAGTGAACCAAGTGTCGGCGGCTAACCGCGGGGCGGCGGTGGGGGCGTATTCGCTGTTCATCGACTTGTCACTGGGGGTGACCGGGCCATTGGTCGGTGCGGTGGCGGCGGGGTTCGGCTTTGCTTCGATGTTTCTGTTCGCGGCGGCTGCGGCAGCGTGCGGGTTGGTGTTGAGCCTCTATTTGTACCGCCAAGGGCGTAGGGACAG  >  NC_002947/5104663‑5104955

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: