Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,082,281 T→A 8.8% *318L (TAA→TTA)  PP_1861 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,082,2810TA 8.8% 161.8 / 5.6 57*318L (TAA→TTA) PP_1861LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base A (2/3);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

TCAGCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCTTCAGG  >  NC_002947/2082137‑2082426
                                                                                                                                                |                                                                                                                                                 
tcaGCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATc                                                                                                                                               <  7:4992/149‑1 (MQ=255)
  aGCACCGAGGGGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGATTAATTcgc                                                                                                                                             <  4:177332/149‑1 (MQ=255)
   gCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATcgcg                                                                                                                                            >  6:210519/1‑149 (MQ=255)
   gCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATcgcg                                                                                                                                            <  8:221703/149‑1 (MQ=255)
     aaCGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGATTAATTCGCGCa                                                                                                                                          <  7:32473/149‑1 (MQ=255)
          ggTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTc                                                                                                                                      >  5:57397/1‑148 (MQ=255)
          ggTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCa                                                                                                                                     <  5:210519/149‑1 (MQ=255)
           gTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCAc                                                                                                                                    <  5:33696/149‑1 (MQ=255)
            tGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCATTGCGCGATTAATTCGCGCactga                                                                                                                                      >  2:58901/1‑142 (MQ=255)
            tGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCATTGCGCGATTAATTCGCGCactga                                                                                                                                      <  1:58901/146‑5 (MQ=255)
                   cAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCaa                                                                                                                             <  4:169178/148‑1 (MQ=255)
                     gggCGAAACAGACCTGGTTGTCAAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGc                                                                                                                          >  2:106649/1‑149 (MQ=255)
                      ggCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAgga                                                                                                                          >  4:18212/1‑147 (MQ=255)
                      ggCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGc                                                                                                                          >  8:49060/1‑148 (MQ=255)
                       gCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGc                                                                                                                          >  1:97387/1‑147 (MQ=255)
                        cGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATACACTGCGCAAGGCCAg                                                                                                                       >  2:228966/1‑149 (MQ=255)
                           aaCAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTc                                                                                                                    >  3:191814/1‑149 (MQ=255)
                                aCCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCAAgg                                                                                                               >  6:89923/1‑149 (MQ=255)
                                  cTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGt                                                                                                             <  5:89923/149‑1 (MQ=255)
                                  cTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGAACGCTGATCATAAGTGCGCGAATTAATCGCGCAATTAAATGAGCAAGGCCAGGTCCCAGCGt                                                                                                             >  7:19675/1‑149 (MQ=255)
                                   tGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTg                                                                                                            >  5:252458/1‑149 (MQ=255)
                                     gTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAg                                                                                                                       <  1:88016/136‑1 (MQ=255)
                                     gTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAg                                                                                                                       >  2:88016/1‑136 (MQ=255)
                                          ccAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCgg                                                                                                     >  8:73520/1‑149 (MQ=255)
                                           cagcagcagCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCggg                                                                                                    <  1:106649/149‑1 (MQ=255)
                                           cagcagcagCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCggg                                                                                                    <  6:57397/149‑1 (MQ=255)
                                                  agcTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGc                                                                                             <  6:138152/149‑1 (MQ=255)
                                                        tCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGctct                                                                                       >  2:319157/1‑149 (MQ=255)
                                                            aggaggCTTGTGTGTCGGCGTTCATGCAGGTCCTTGAAACTCATAAATGGAAAATTCTGCGACGCTGATCATCAGTTCGCGAATTAATCGCGCAATTCACTGCGCAATGCCAGGTCCCAGGGTGTAATCGGGCTTAAGCGGCTCTTGat                                                                                   <  2:80946/149‑2 (MQ=255)
                                                            aggaggCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAg                                                                                   >  8:215323/1‑149 (MQ=255)
                                                                   ttGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTc                                                                             >  8:18719/1‑148 (MQ=255)
                                                                       gtgtCGGCGTTCATGCAGCTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATAGCGCAATTCACTGCGCAAGGACAGGTCCCAGGGTGGAATCGGGCTGAAGCGCCACTTCAGGAATTCGagc                                                                         >  2:72941/1‑148 (MQ=255)
                                                                               gTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACgg                                                                >  7:299240/1‑149 (MQ=255)
                                                                                  cATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGa                                                                                    <  7:87927/126‑1 (MQ=255)
                                                                                  cATGCAGCTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGa                                                                                    >  8:87927/1‑126 (MQ=255)
                                                                                       aGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACAGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCga                                                          >  3:304312/1‑147 (MQ=255)
                                                                                                  aGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCAt                                             <  1:319157/149‑1 (MQ=255)
                                                                                                      atcaATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTg                                                                                     >  7:156637/1‑105 (MQ=255)
                                                                                                      atcaATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTg                                                                                     <  8:156637/105‑1 (MQ=255)
                                                                                                      atcaATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGCCTCCGAGAGTTCGTTTCATGTGc                                         >  5:190660/1‑149 (MQ=255)
                                                                                                           tGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGc                                    >  2:182623/1‑149 (MQ=255)
                                                                                                              aaaaTTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCa                                       >  5:271469/1‑143 (MQ=255)
                                                                                                              aaaaTTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCa                                       <  6:271469/143‑1 (MQ=255)
                                                                                                              aaaaTTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCg                                   >  3:293813/1‑147 (MQ=255)
                                                                                                               aaaTTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAg                                <  7:18719/149‑1 (MQ=255)
                                                                                                                  ttCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCAt                                             <  3:250150/133‑1 (MQ=255)
                                                                                                                  ttCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCAt                                             >  4:250150/1‑133 (MQ=255)
                                                                                                                       gCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGt                           <  7:73520/146‑1 (MQ=255)
                                                                                                                        cACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAAt                       <  4:191814/149‑1 (MQ=255)
                                                                                                                        cACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAAt                       >  4:109228/1‑149 (MQ=255)
                                                                                                                              gATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCt                 >  5:254122/1‑149 (MQ=255)
                                                                                                                               atcatcAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTg                <  1:182623/149‑1 (MQ=255)
                                                                                                                                    aattGCGCGATTAATTCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGtctc           >  3:290264/4‑149 (MQ=255)
                                                                                                                                       tGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCCGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAATCGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGc        <  2:97387/149‑1 (MQ=255)
                                                                                                                                         cgcgAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCtt      <  4:293813/149‑1 (MQ=255)
                                                                                                                                           cgATTTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCGAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCTTCa    <  1:72941/149‑1 (MQ=255)
                                                                                                                                             aaTTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCTTCAgg  <  7:215323/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TCAGCAACGAGGTGGAGTGCAGGGCGAAACAGACCTGGTTGTCCAGCAGCAGCTCGTCACAGGAGGCTTGGGTGTCGGCGTTCATGCAGGTCCTTGAAAGTCATCAATGGAAAATTCTGGCACGCTGATCATCAGTGCGCGAATTAATCGCGCAATTCACTGCGCAAGGCCAGGTCCCAGGGTGGAATCGGGCTGAAGCGGCTCTTGAGGAATTCGAGCAGCAAACGGCTGCGAGAGTTCGTTTCATGTGCCAGGCGTAGCGCGTAAATGCCGCTGGTCTCAGCTTCAGG  >  NC_002947/2082137‑2082426

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: