Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 798,662 C→A 8.8% intergenic (+17/+10) fklB‑I → / ← PP_0685 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947798,6620CA 8.8% 152.9 / 4.6 53intergenic (+17/+10)fklB‑I/PP_0685FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base A (2/3);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

CAGCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTTCG  >  NC_002947/798514‑798810
                                                                                                                                                    |                                                                                                                                                     
cAGCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAggccggcc                                                                                                                                                       >  6:167851/1‑149 (MQ=255)
  gCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCt                                                                                                                                                     <  8:285062/149‑1 (MQ=255)
     gatgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGcc                                                                                                                                                  >  1:350669/1‑149 (MQ=255)
      atgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGcc                                                                                                                                                  <  2:359390/148‑1 (MQ=255)
         aaCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGc                                                                                                                                               <  6:28410/148‑1 (MQ=255)
         aaCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGggccggcc                                                                                                                                              <  4:25233/149‑1 (MQ=255)
                gcagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTccag                                                                                                                                        <  1:74207/147‑4 (MQ=255)
                ccagcaAATGGCGCCTGTCCGGGCCAAGCGAGCGGGCTAACGGGCCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCGGGTGTTCGACGTCGAGCTGCTCTAGGTTCTGTAATGCCCTGGAGGCCGGCCCTCCCGGCCTCAATTc                                                                                                                                      <  8:314560/148‑1 (MQ=255)
                 cagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTc                                                                                                                                      >  6:315143/1‑148 (MQ=255)
                  agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                    <  6:291430/149‑1 (MQ=255)
                    caAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAg                                                                                                                                   >  5:252382/1‑148 (MQ=255)
                     aaaTGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTATAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTc                                                                                                                                 >  6:40105/1‑149 (MQ=255)
                      aaTGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCGCAAGTCCAGTCt                                                                                                                                >  1:265177/1‑149 (MQ=255)
                                 tACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAg                                                                                                                      <  5:167851/148‑1 (MQ=255)
                                 tACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCATGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAgg                                                                                                                     >  7:323226/1‑149 (MQ=255)
                                   cGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTATTGCCCTGGCGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAggg                                                                                                                    <  8:28936/148‑1 (MQ=255)
                                              aGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                    >  6:134073/1‑121 (MQ=255)
                                              aGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                    <  5:134073/121‑1 (MQ=255)
                                                 tGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCgcgc                                                                                                      >  8:256519/1‑148 (MQ=255)
                                                   gCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCcgc                                                                                                        >  3:93428/1‑144 (MQ=255)
                                                   gCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCcgc                                                                                                        <  4:93428/144‑1 (MQ=255)
                                                     ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcc                                                                                                                                          >  8:233394/1‑107 (MQ=255)
                                                     ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcc                                                                                                                                          <  7:233394/108‑2 (MQ=255)
                                                       aCGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAg                                                                                                  >  3:193550/1‑146 (MQ=255)
                                                         ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaa                                                                                             >  2:24516/1‑149 (MQ=255)
                                                         ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaa                                                                                             >  7:142276/1‑149 (MQ=255)
                                                                         gcagcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTTCCGGACTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGCTTCCg                                                                              >  7:234890/1‑148 (MQ=255)
                                                                           agcgTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCATCTCGACCTTCTGTAATACCATGGAGGCCGGCACTGCCGGCCTCAATTCAAGTATGTCCCCCCAGCGCGCAACGCCCCCGCATAGCAGCACAGCAACAGATACCGTAc                                                                           >  2:97652/1‑149 (MQ=255)
                                                                           agcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGTCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTAc                                                                           <  5:130971/149‑1 (MQ=255)
                                                                               tCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGt                                                                       >  8:143456/1‑149 (MQ=255)
                                                                                 ccgccgCACAGCGTTCTGGTGTTCGACGTCTCGCTGCTCGAGGTTCTGAAATGCCCTGGCGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTc                                                                     <  2:265177/149‑1 (MQ=255)
                                                                                 ccgccgCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCa                                                                                                                       >  3:10571/1‑99 (MQ=255)
                                                                                 ccgcagCACAGCGTTCTTGGGGTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGtt                                                                      <  6:252382/148‑1 (MQ=255)
                                                                                   gccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCt                                                                   <  7:256519/149‑1 (MQ=255)
                                                                                    ccgcACAGCGTTCTGGTGTTCGACGTCTAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCtt                                                                  <  5:315143/149‑1 (MQ=255)
                                                                                    ccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCtt                                                                  >  6:259976/1‑149 (MQ=255)
                                                                                       cacaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAg                                                               >  3:211917/1‑149 (MQ=255)
                                                                                             gTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGc                                                              <  2:2336/144‑1 (MQ=255)
                                                                                             gTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGc                                                              >  1:2336/1‑144 (MQ=255)
                                                                                                   gtgtTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGtt                                                    <  5:259976/148‑1 (MQ=255)
                                                                                                          ccggcgtgctgctCCCCGTTCTGAATTGCCCTGGCGGCCGGCCTGCCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGGCCGCGCATAGCAGACCAGAAACCGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                            <  7:302977/142‑1 (MQ=255)
                                                                                                          aCGTCGTCCTGCTCGACGTTCTGTAAGGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                            <  2:113358/149‑1 (MQ=255)
                                                                                                                     ctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCt                                 <  4:211917/149‑1 (MQ=255)
                                                                                                                        gACGTTCTATAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTc                                <  5:40105/147‑1 (MQ=255)
                                                                                                                           gTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTc                           <  2:80651/149‑1 (MQ=255)
                                                                                                                           gTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTc                           >  1:80651/1‑149 (MQ=255)
                                                                                                                                 tAATGCCCTGGAGGCCGGCCCTGCCGGCCTTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGt                      >  1:214949/1‑149 (MQ=255)
                                                                                                                                 tAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTGCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACGCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                        <  8:234890/146‑1 (MQ=255)
                                                                                                                                 tAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGTACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTc                     >  3:23767/1‑149 (MQ=255)
                                                                                                                                    tGCCCTGGAGGCCGCCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCAAAGCAGAACAGAAACAGATTCCGTACCAGTTCCATGAGCACACCCGGTTCACTCGAATAAAGCCCGCTCATGTCCAGGTCAcc                  >  4:248387/1‑149 (MQ=255)
                                                                                                                                      cccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGTACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCAccc                 <  4:23767/148‑1 (MQ=255)
                                                                                                                                       ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCAccc                 <  1:61193/147‑1 (MQ=255)
                                                                                                                                       ccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCCCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCAccc                 >  2:61193/1‑147 (MQ=255)
                                                                                                                                        cTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTAAGc                                      >  6:27775/1‑125 (MQ=255)
                                                                                                                                                  gccCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGtt    >  3:77694/1‑149 (MQ=255)
                                                                                                                                                    ccTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTtcg  >  1:317557/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTTCG  >  NC_002947/798514‑798810

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: