Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,549 G→A 11.1% G138D (GGC→GAC)  PP_4611 → ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,233,5490GA11.1% 141.8 / 6.6 54G138D (GGC→GAC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base G (29/19);  new base A (4/2);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTT  >  NC_002947/5233401‑5233695
                                                                                                                                                      |                                                                                                                                                  
cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACgg                                                                                                                                                    >  6:348010/1‑149 (MQ=255)
   cAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCt                                                                                                                                                 >  8:185284/1‑149 (MQ=255)
      gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTggg                                                                                                                                              <  7:185284/149‑1 (MQ=255)
      gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTggg                                                                                                                                              <  3:173842/149‑1 (MQ=255)
       cgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGc                                                                                                                                             <  3:221854/149‑1 (MQ=255)
        cctcgcacggcgcgacctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGc                                                                                                                                             >  1:364120‑M2/38‑148 (MQ=255)
          tcgcacggcgcgacctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGcc                                                                                                                                            >  5:244031‑M2/36‑147 (MQ=255)
            GCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGgccagca                                                                                                                                        >  2:12039‑M1/1‑146 (MQ=255)
                 ggTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                        <  4:151845/144‑4 (MQ=255)
                 ggTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                        >  3:151845/1‑141 (MQ=255)
                   tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                                 >  3:115201/1‑149 (MQ=255)
                   gcgacctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                                 <  2:364120‑M2/123‑1 (MQ=255)
                       ggAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                             >  6:119135/1‑149 (MQ=255)
                       ggAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                             >  4:267536/1‑149 (MQ=255)
                       cctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                             <  6:244031‑M2/127‑1 (MQ=255)
                                cgcgCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                                    <  5:119135/149‑1 (MQ=255)
                                   gcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTGGGAGCCCTTGGTAAAAACCGCCATAAATCTGTCCCACCTCGACCGCCAGGGGCAAAGCGAAATAGCTAGCGAACTgtg                                                                                                                   >  1:294070/1‑147 (MQ=255)
                                          gagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCt                                                                                                          >  3:52603‑M2/4‑149 (MQ=255)
                                                 aGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACg                                                                                                   >  4:129183/1‑149 (MQ=255)
                                                      gCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCaa                                                                                              <  4:52603/149‑1 (MQ=255)
                                                      gCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCaa                                                                                              <  3:267536/149‑1 (MQ=255)
                                                            atcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                                 >  4:321032/1‑140 (MQ=255)
                                                            atcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                                 <  3:321032/140‑1 (MQ=255)
                                                                   agaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGTTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                                 <  4:115201/149‑1 (MQ=255)
                                                                                          ttGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                          >  4:350440/1‑149 (MQ=255)
                                                                                              tGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcgcg                                                      >  3:348176/1‑149 (MQ=255)
                                                                                                 aCGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtg                                                   >  1:362084/1‑149 (MQ=255)
                                                                                                  cGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                                  <  3:129183/149‑1 (MQ=255)
                                                                                                     ccTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcgc                                                                  <  1:283006/130‑1 (MQ=255)
                                                                                                     ccTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcgc                                                                  >  2:283006/1‑130 (MQ=255)
                                                                                                      cTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAAc                                              >  2:242089/1‑149 (MQ=255)
                                                                                                            a‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                                         >  5:52794/1‑149 (MQ=255)
                                                                                                            a‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                                         >  1:14713/1‑149 (MQ=255)
                                                                                                                    gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcc                                  >  2:298919/1‑149 (MQ=255)
                                                                                                                     tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGGAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTccc                                 >  2:97994/1‑149 (MQ=255)
                                                                                                                     tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTccc                                 >  1:345164/1‑149 (MQ=255)
                                                                                                                       aaaCGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcccc                                <  5:348010/148‑1 (MQ=255)
                                                                                                                             gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggtg                         >  7:199845/1‑149 (MQ=255)
                                                                                                                               cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc                       >  7:239425/1‑149 (MQ=255)
                                                                                                                               cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc                       >  7:336729/1‑149 (MQ=255)
                                                                                                                                ttCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGTGCATTTCGCTGGCCACGGTCAAGCGCt                                                                                         <  2:61354/82‑1 (MQ=255)
                                                                                                                                ttCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCt                                                                                         >  1:61354/1‑82 (MQ=255)
                                                                                                                                 tCCTTCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCCGCCCTGCGCTGCTACTTGGCGCTGTGAACGCCTCGTTTGCCCCGCTGGTGGCg                      >  2:233072/1‑148 (MQ=255)
                                                                                                                                  ccTGCTGGCCCGGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                          <  5:246570/111‑1 (MQ=255)
                                                                                                                                  ccTGCTGGCCCGGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                          >  6:246570/1‑111 (MQ=255)
                                                                                                                                  ccTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg                      <  3:350440/147‑1 (MQ=255)
                                                                                                                                     cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa                   >  5:26685/3‑147 (MQ=255)
                                                                                                                                       tgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAgg               >  3:341295/1‑149 (MQ=255)
                                                                                                                                         gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGc              <  2:294070/148‑1 (MQ=255)
                                                                                                                                          cccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCg            >  6:92283/1‑149 (MQ=255)
                                                                                                                                               gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCgg                     >  7:300725/1‑135 (MQ=255)
                                                                                                                                               gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCgg                     <  8:300725/135‑1 (MQ=255)
                                                                                                                                                  gACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTgg    >  7:110934/1‑149 (MQ=255)
                                                                                                                                                   aCGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGt   <  4:348176/149‑1 (MQ=255)
                                                                                                                                                   aCGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGt   <  2:345164/149‑1 (MQ=255)
                                                                                                                                                    gggCGTGGGCCAAGGGGTATTGGCTCGCTAAGTGGGAATTTCGCTGGCCAGGGGCAAGGGCTACCTCAAAACGGCTGCAATGGGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGtt  <  1:97994/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTT  >  NC_002947/5233401‑5233695

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: ‑