Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,240,754 (C)5→4 coding (915/1041 nt) galS ← DNA‑binding transcriptional dual regulator GalS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,240,7500C.98.2% 228.4 / ‑3.0 56coding (919/1041 nt)galSDNA‑binding transcriptional dual regulator GalS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (28/27);  minor base A (1/0);  total (29/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 = 1.00e+00

TCACAGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGAATCG  >  NC_000913/2240611‑2240884
                                                                                                                                           |                                                                                                                                      
tcacAGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCtgctgc                                                                                                                                        <  2:555897/140‑1 (MQ=255)
 cacaGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcc                                                                                                                                       <  1:698034/140‑1 (MQ=255)
      aTTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccctgc                                                                                                                                  <  2:615439/140‑4 (MQ=255)
      aTTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccctgc                                                                                                                                  >  1:484060/1‑137 (MQ=255)
        tGAAAGCAGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGGCCCTGCaa                                                                                                                               >  2:396145/1‑140 (MQ=255)
              cGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAg                                                                                                                         >  2:117917/1‑140 (MQ=255)
                gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGtt                                                                                                                       >  1:657033/1‑140 (MQ=255)
                gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGtt                                                                                                                       >  2:490425/1‑140 (MQ=255)
                gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGtt                                                                                                                       >  2:422823/1‑140 (MQ=255)
                gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACACTGGCTGGCCCGAGGATTAATATTGCCTGCTACCCCTGCAAGGCCAGtt                                                                                                                       >  2:255876/1‑140 (MQ=255)
                  tACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCg                                                                                                                     <  2:421002/140‑1 (MQ=255)
                         ccTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTa                                                                                                              >  2:240805/1‑140 (MQ=255)
                          cTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTaa                                                                                                             >  2:645734/1‑140 (MQ=255)
                           tGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAAt                                                                                                            >  1:547485/1‑140 (MQ=255)
                             aTTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGATGATCAATATTGCGTGCCGCCCCTGCAAGGCTAGTTCGGTGGCTAAttt                                                                                                          >  2:305436/1‑140 (MQ=255)
                             aTTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAAttt                                                                                                          <  1:645734/140‑1 (MQ=255)
                                  ttgAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCa                                                                                                     >  2:17695/1‑140 (MQ=255)
                                  ttgAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCa                                                                                                     >  2:49265/1‑140 (MQ=255)
                                    gAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCAtt                                                                                                   >  1:736901/1‑140 (MQ=255)
                                      gTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                                 <  1:490425/140‑1 (MQ=255)
                                          gTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCa                                                                                             <  1:647681/140‑1 (MQ=255)
                                           tGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCaa                                                                                            <  1:147852/140‑1 (MQ=255)
                                               cgccgccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATggg                                                                                        >  2:671324/1‑140 (MQ=255)
                                                gccgccaaATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAACCAGTGGTTGGCACGCGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGa                                                                                       <  2:453868/140‑1 (MQ=255)
                                                 ccgccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGAt                                                                                      >  1:736634/1‑140 (MQ=255)
                                                     cgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTTCAAGGCCAATTCAGAGACTAATTTCTCCATATAAGCAATGGGATAACg                                                                                  >  2:331029/1‑140 (MQ=255)
                                                                   ttGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAAtt                                                                    <  1:17695/140‑1 (MQ=255)
                                                                   ttGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAAtt                                                                    >  2:51228/1‑140 (MQ=255)
                                                                    tGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTg                                                                   <  1:97897/140‑1 (MQ=255)
                                                                      cGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCg                                                                 >  1:37988/1‑140 (MQ=255)
                                                                        aCAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCggg                                                               <  1:597153/140‑1 (MQ=255)
                                                                        aCAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCggg                                                               >  1:403945/1‑140 (MQ=255)
                                                                           gAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCg                                                            <  1:32620/140‑1 (MQ=255)
                                                                           gAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCg                                                            >  2:471763/1‑140 (MQ=255)
                                                                           gAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCg                                                            <  1:240805/140‑1 (MQ=255)
                                                                            aaTGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTAAAGTAATGGGATAACGCACGGTCGTTAATTGCGGGTCgg                                                           >  2:272253/1‑140 (MQ=255)
                                                                                 cGACGCACTAACGTCGGCATAATACAGTGACTGGCACAAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTaa                                                      <  1:331029/140‑1 (MQ=255)
                                                                                 cGACGCACTAACGTCGGCATAAAACGGTGTCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTaa                                                      <  1:317864/140‑1 (MQ=255)
                                                                                      cACTAACGTCGGCATCAAACAGTGGATGGTACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                                 <  1:305436/140‑1 (MQ=255)
                                                                                      cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                                 >  1:651552/1‑140 (MQ=255)
                                                                                       aCTAACGTCGGCATAAAACAGTGGCTGTCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCg                                                <  1:736147/140‑1 (MQ=255)
                                                                                       aCTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCg                                                <  2:403945/140‑1 (MQ=255)
                                                                                        cTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGa                                               <  2:31144/140‑1 (MQ=255)
                                                                                         gaaCGTCGGCATAAAACGGTGGTTGGCACGAAGATGAATTTTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGAt                                              <  2:298195/139‑1 (MQ=255)
                                                                                          aaCGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATg                                             >  2:615983/1‑140 (MQ=255)
                                                                                             gTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGa                                          <  1:587878/140‑1 (MQ=255)
                                                                                                     aaaaCAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGTTCGGTGTAACGGGTGATGGGAATatcatc                                  >  1:298195/1‑140 (MQ=255)
                                                                                                      aaaCAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCg                                 <  1:471763/140‑1 (MQ=255)
                                                                                                       aaCAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGa                                >  2:407565/1‑140 (MQ=255)
                                                                                                            tggctggcACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCg                           <  1:51228/140‑1 (MQ=255)
                                                                                                             ggctggcACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCga                          <  1:558591/140‑1 (MQ=255)
                                                                                                             ggctggcACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCga                          <  1:117917/140‑1 (MQ=255)
                                                                                                                tggcACGAGGATCAATATTGCCTGCTGCCCCTGCAATGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCgatga                       <  1:650536/140‑1 (MQ=255)
                                                                                                                 ggcACGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCgatgat                      >  2:494657/1‑140 (MQ=255)
                                                                                                                     cGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTgag                  <  1:49265/140‑1 (MQ=255)
                                                                                                                     cGAGGATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTgag                  <  2:547485/140‑1 (MQ=255)
                                                                                                                        ggATCAATATTGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTgagaga               <  2:651552/140‑1 (MQ=255)
                                                                                                                                  tGCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGaa     >  2:420114/1‑140 (MQ=255)
                                                                                                                                   gCCTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAAATTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGAAt    >  2:447134/1‑140 (MQ=255)
                                                                                                                                    ccTGCTGCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGAATc   >  2:86670/1‑140 (MQ=255)
                                                                                                                                     ctgctgCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGAATCg  >  2:178851/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                      
TCACAGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGAATCG  >  NC_000913/2240611‑2240884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: