Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,223,113 C→T 7.1% intergenic (‑420/+140) PP_2823 ← / ← PP_2824 methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,223,1130CT 7.1% 157.2 / 2.4 56intergenic (‑420/+140)PP_2823/PP_2824methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/23);  new base T (2/2);  total (31/25)
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.17e-01

GCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGA  >  NC_002947/3222976‑3223258
                                                                                                                                         |                                                                                                                                                 
gcCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGTGCTATCGGGGc                                                                                                                                        >  4:214917/1‑149 (MQ=255)
gcCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGc                                                                                                                                        <  7:57211/149‑1 (MQ=255)
gcCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGc                                                                                                                                        <  3:325459/149‑1 (MQ=255)
gcCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTGATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACCTGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGTGCTATCGGGGc                                                                                                                                        >  4:219224/1‑149 (MQ=255)
       tGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGc                                                                                                                                 >  5:353467/1‑149 (MQ=255)
         cTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCtt                                                                                                                               >  8:47951/1‑149 (MQ=255)
         cTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCtt                                                                                                                               <  5:58356/149‑1 (MQ=255)
         cTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCtt                                                                                                                               <  3:350891/149‑1 (MQ=255)
          tAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTg                                                                                                                              >  2:258249/1‑149 (MQ=255)
          tAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTg                                                                                                                              >  3:112907/1‑149 (MQ=255)
              gCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGtt                                                                                                                          >  2:242833/1‑149 (MQ=255)
               cccATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTCTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTg                                                                                                                         <  4:17059/149‑1 (MQ=255)
                       cGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGc                                                                                                                  <  7:47951/148‑1 (MQ=255)
                               ccACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCgg                                                                                                         <  2:123664/149‑1 (MQ=255)
                                      gTTGTGTTTTATTTGTTATATCACGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGTGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGggcggc                                                                                                   <  3:219224/148‑1 (MQ=255)
                                      gTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGTGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGggcggc                                                                                                   <  3:214917/148‑1 (MQ=255)
                                        tgtgTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGt                                                                                                >  7:76417/1‑149 (MQ=255)
                                            ttttATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTgcg                                                                                            >  6:218351/1‑149 (MQ=255)
                                                aTTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGcga                                                                                         <  7:13249/148‑1 (MQ=255)
                                                 tttGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGcgacg                                                                                       <  6:353467/149‑1 (MQ=255)
                                                 tttGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGcgac                                                                                        >  3:201524/1‑148 (MQ=255)
                                                    gTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCt                                                                                                                                      >  8:185637/1‑99 (MQ=255)
                                                    gTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCt                                                                                                                                      <  7:185637/99‑1 (MQ=255)
                                                    gTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGcgacga                                                                                      >  1:307143/1‑147 (MQ=255)
                                                            aaGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAggcgg                                                                            >  1:1105/1‑149 (MQ=255)
                                                                          aTACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCg                                                              >  8:227599/1‑149 (MQ=255)
                                                                              aGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCCGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGttt                                                          >  5:291966/1‑149 (MQ=255)
                                                                                     ttGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTggg                                                    <  4:112907/148‑1 (MQ=255)
                                                                                     ttGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTggg                                                    <  1:258249/148‑1 (MQ=255)
                                                                                         ttACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGa                                               >  6:123365/1‑149 (MQ=255)
                                                                                         ttACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGa                                               >  2:79798/1‑149 (MQ=255)
                                                                                           aCATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACg                                             >  5:14876/1‑149 (MQ=255)
                                                                                           aCATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACg                                             >  4:107202/1‑149 (MQ=255)
                                                                                           aCATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACg                                             >  2:74442/1‑149 (MQ=255)
                                                                                               gCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGa                                                                 >  4:268134/1‑125 (MQ=255)
                                                                                               gCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGa                                                                 <  3:268134/125‑1 (MQ=255)
                                                                                               gCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCtt                                         >  2:114081/1‑149 (MQ=255)
                                                                                               gCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCtt                                         >  7:290840/1‑149 (MQ=255)
                                                                                                 aaGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTcc                                       <  4:77858/149‑1 (MQ=255)
                                                                                                  aGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTc                                        >  6:242466/1‑147 (MQ=255)
                                                                                                   gACGCAATGGCGCCTGAGCCTGAAAAAAGCTGGGAAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGAGCGGCAAGGCGCGGGGCGGAAGTTGCGCTACTACGTCAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAg                                     <  1:74442/149‑1 (MQ=255)
                                                                                                   gACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAg                                     >  4:65529/1‑149 (MQ=255)
                                                                                                        aaTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAg                                >  5:359351/1‑149 (MQ=255)
                                                                                                        aaTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAg                                >  8:340071/1‑149 (MQ=255)
                                                                                                         aTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAgg                               >  3:264682/1‑149 (MQ=255)
                                                                                                          tGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGgc                              <  6:359351/149‑1 (MQ=255)
                                                                                                            gcgcCAGAGCCTGCAAAAAGCACGGCAGGCGCGATCGGGGCCTGGTGCTAGGGCTGGATCGGCAAGGCGCGGGGCGGCAGCTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCTGGAAGgcgt                            <  8:123710/149‑2 (MQ=255)
                                                                                                                       tGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCg                 >  7:88747/1‑149 (MQ=255)
                                                                                                                        gCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGAGCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                <  5:218351/149‑1 (MQ=255)
                                                                                                                         cAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCTAGGGAAGGCGCCGGGTCTACCg                 <  7:227599/147‑1 (MQ=255)
                                                                                                                          aaaaaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                >  4:188902/1‑147 (MQ=255)
                                                                                                                          aaaaaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGa              <  2:307143/149‑1 (MQ=255)
                                                                                                                             aaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGc           >  7:66694/1‑149 (MQ=255)
                                                                                                                                    gcaggCGCGATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGGTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAg    <  1:242833/149‑1 (MQ=255)
                                                                                                                                     caggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAgg   <  7:340071/149‑1 (MQ=255)
                                                                                                                                      aggCGCTTTCGGGGCCTGGTGCTTGGGTTGTATCGGCAAGGCGCGGTGCGGCTGTTGCGCTACGACGGAAGGCGGCTAAGACACGAGCTGTGTGGTGTGGGTGACGCCTTCCAGTGATGTCGCCGGGTCTACCGGGAGGCGTCTTAGGa  <  8:290840/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
GCCTTGGTGCTAAAGCCCATTTTCGCTCCTGCCACTTGGTTGTGTTTTATTTGTTATATCAAGTGTTTTGTAAAATACAGGAGTTTTGCTTACATGCAAGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGA  >  NC_002947/3222976‑3223258

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: