Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,941,695 C→T 7.3% intergenic (+15/+31) PP_2574 → / ← PP_2575 LysR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,941,6950CT 7.3% 156.6 / 2.6 55intergenic (+15/+31)PP_2574/PP_2575LysR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/28);  new base T (2/2);  total (25/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 = 7.65e-01

ATGCTCGACATTGCCGAGCCGGTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCGCGCCGGTCCTGGCGCCGGTCCTGACGCTGGCG  >  NC_002947/2941549‑2941842
                                                                                                                                                  |                                                                                                                                                    
aTGCTCGACATTGCCGAGCCGGTGGCGCTGTCGGGGCTGCCTGCGGCGTGGATGAAGCGAGTCCAGCTGACCAAACCCGCGCACTTGTTAGTGGAGTTTTGTCGGGATGAGCTCGCGAAGTTCCTCAAGTGAGCACGCCATGCCAcccc                                                                                                                                                    >  3:172526/1‑149 (MQ=255)
    tCGACATTGCCGAGCCGGTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCgg                                                                                                                                                  >  2:155321/1‑147 (MQ=255)
      gACATTGCCGAGCCGGTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCa                                                                                                                                                <  1:155321/147‑1 (MQ=255)
      gACATTGCCGAGCCGGTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAg                                                                                                                                               >  1:210363/1‑148 (MQ=255)
             ccGAGCCGGTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAg                                                                                                                                               <  1:236452/141‑1 (MQ=255)
             ccGAGCCGGTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAg                                                                                                                                               >  2:236452/1‑141 (MQ=255)
                     gTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCAACCCGGCAGGAGTAGCAtggcgtgg                                                                                                                               >  8:196595/1‑149 (MQ=255)
                              tCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCt                                                                                                                       >  1:29093/1‑148 (MQ=255)
                                ggggCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAg                                                                                                                   >  2:231326/1‑149 (MQ=255)
                                  ggctggctGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAggg                                                                                                                 <  2:29093/149‑1 (MQ=255)
                                  ggctggctGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAggg                                                                                                                 <  8:119108/149‑1 (MQ=255)
                                      ggctgCGGCGTGGTTGAAGCGAGTCCAGCTGACCACACCGGCCCAGTTGGTCGTGGAGTTTTGTCTGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACACCGGCAGGAGTAGCATGGCGTGTCGATTGCGCAtgggat                                                                                                              >  7:85899/1‑146 (MQ=255)
                                       gctgCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGa                                                                                                            >  2:72479/1‑149 (MQ=255)
                                              cgTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTa                                                                                                     <  5:34665/149‑1 (MQ=255)
                                                tGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGc                                                                                                   <  3:98534/149‑1 (MQ=255)
                                                   ttGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCcgc                                                                                                >  5:234349/1‑149 (MQ=255)
                                                    tGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCcgcg                                                                                               >  6:150131/1‑149 (MQ=255)
                                                                     aCCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACGTgg                                                                                                       <  8:194573/124‑1 (MQ=255)
                                                                     aCCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTgg                                                                                                       >  7:194573/1‑124 (MQ=255)
                                                                     aCCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCCGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGc                                                                              >  4:22817/1‑149 (MQ=255)
                                                                      ccAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCt                                                                             <  6:100335/149‑1 (MQ=255)
                                                                        aaaCCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGt                                                                           <  1:72479/149‑1 (MQ=255)
                                                                        aaaCCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCCGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGt                                                                           <  3:22817/149‑1 (MQ=255)
                                                                          aCCGGCGCAGTTGTTCGGGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAg                                                                         <  7:219382/149‑1 (MQ=255)
                                                                                 cAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTgg                                                                  <  2:210363/149‑1 (MQ=255)
                                                                                  aGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCTACTCCTGCCGGGGTGGCAtg                                                                                                                                     <  3:127585/82‑1 (MQ=25)
                                                                                  aGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCTACTCCTGCCGGGGTGGCAtg                                                                                                                                     >  4:127585/1‑82 (MQ=25)
                                                                                  aGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTgg                                                                  >  1:47264/1‑148 (MQ=255)
                                                                                  aGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGc                                                                 <  8:50253/149‑1 (MQ=255)
                                                                                    ttgttCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCa                                                                <  1:206194/149‑1 (MQ=255)
                                                                                       ttCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATgg                                                             <  7:80312/148‑1 (MQ=255)
                                                                                         cGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTgg                                                                  >  4:212949/1‑141 (MQ=255)
                                                                                         cGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTgg                                                                  <  3:212949/141‑1 (MQ=255)
                                                                                            ggAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAggg                                                                                                                 <  6:235415/91‑1 (MQ=255)
                                                                                            ggAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAggg                                                                                                                 >  5:235415/1‑91 (MQ=255)
                                                                                             gAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTggg                                                      >  2:26416/1‑149 (MQ=255)
                                                                                               gTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGtt                                                    >  4:146429/1‑149 (MQ=255)
                                                                                                ttttGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAg                                                                         >  3:135399/1‑127 (MQ=255)
                                                                                                ttttGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAg                                                                         <  4:135399/127‑1 (MQ=255)
                                                                                                  ttGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGtttt                                                  <  7:196595/148‑1 (MQ=255)
                                                                                                   tGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTg                                                 <  3:146429/148‑1 (MQ=255)
                                                                                                   tGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTg                                                 <  2:47264/148‑1 (MQ=255)
                                                                                                   tGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGAAGGCACTGGCATGGGCCTGGGTTTTgg                                                <  1:231326/149‑1 (MQ=255)
                                                                                                     tCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGcc                                              <  5:150131/149‑1 (MQ=255)
                                                                                                       gggAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTc                                            <  6:234349/149‑1 (MQ=255)
                                                                                                          aaGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTc                                         >  8:88949/1‑149 (MQ=255)
                                                                                                                  gCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGtttt                                                  <  7:107631/132‑1 (MQ=255)
                                                                                                                  gCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGtttt                                                  >  8:107631/1‑132 (MQ=255)
                                                                                                                            tGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCgcgccggtcctg                       >  6:118647/1‑149 (MQ=255)
                                                                                                                                   aGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCgcgccggtcctggcgccgg                >  4:254524/1‑149 (MQ=255)
                                                                                                                                   aGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGCTTTGGCCTCGTCCTGGTCgcgccggtcctggcgccg                 <  7:88949/148‑1 (MQ=255)
                                                                                                                                        gCCATGCCACCCCGGCAGGAGTTGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCgcgccggtcctggcgccggtcctg           <  5:215321/149‑1 (MQ=255)
                                                                                                                                         ccATGCTACTCCTGCCGGGGTGGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTg                                                 <  1:54425/110‑1 (MQ=38)
                                                                                                                                         ccATGCTACTCCTGCCGGGGTGGCATGGCGTGGCGATTGCTAAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCCCTGTAGGCACTGGCATGGGCCAGGCTTTTg                                                 >  2:54425/1‑110 (MQ=25)
                                                                                                                                          cATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCgcgccggtcctggcgccg                 <  2:31286/141‑1 (MQ=255)
                                                                                                                                          cATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCgcgccggtcctggcgccg                 >  1:31286/1‑141 (MQ=255)
                                                                                                                                                 ccccGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCGCGCCGGTCCTGGCGCCGGTCCTGACGCTGgcg  <  2:160161/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATGCTCGACATTGCCGAGCCGGTGGCGCTGTCGGGGCTGGCTGCGGCGTGGTTGAAGCGAGTCCAGCTGACCAAACCGGCGCAGTTGTTCGTGGAGTTTTGTCGGGAAGAGCTGGCGAAGTTCCTGAAGTGAGCACGCCATGCCACCCCGGCAGGAGTAGCATGGCGTGGCGATTGCTCAGGGGTTGACCTGGTAGCCGCGCCCTTGCAGGCAACCGCTGTAGGCACTGGCATGGGCCTGGGTTTTGGCCTCGTCCTGGTCGCGCCGGTCCTGGCGCCGGTCCTGACGCTGGCG  >  NC_002947/2941549‑2941842

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