Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,529,044 G→C 10.4% intergenic (+128/‑279) PP_2217 → / → PP_2218 enoyl‑CoA hydratase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,529,0440GC10.4% 171.0 / 11.2 58intergenic (+128/‑279)PP_2217/PP_2218enoyl‑CoA hydratase/transposase
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base C (3/3);  total (25/33)
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.18e-01

GCACAGTTCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGTTACTTTTTGTCGCGTT  >  NC_002947/2528896‑2529186
                                                                                                                                                    |                                                                                                                                               
gCACAGTTCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGa                                                                                                                                                <  1:800709/150‑1 (MQ=255)
  acaGTTCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATg                                                                                                                                              >  2:228131/1‑150 (MQ=255)
   caGTTCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGc                                                                                                                                             <  1:627186/150‑1 (MQ=255)
   caGTTCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGc                                                                                                                                             <  1:810156/150‑1 (MQ=255)
     gTTCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGccc                                                                                                                                           >  1:138539/1‑150 (MQ=255)
      ttCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCg                                                                                                                                          >  2:1041299/1‑150 (MQ=255)
                gTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGCGGCGGGCATGGATAACCGTAAAGCGCTg                                                                                                                                >  2:1296745/1‑150 (MQ=17)
                         aCTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCat                                                                                                                      <  1:1348317/150‑1 (MQ=255)
                           tCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCAGAGCGTAGGGTAAGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatgg                                                                                                                    >  1:674168/1‑150 (MQ=255)
                                    tgccgttgccgttgccgtGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCATGGGGCCt                                                                                                           <  2:1114402/150‑1 (MQ=255)
                                        gttgccgttgccgtGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAg                                                                                                       >  2:587613/1‑150 (MQ=255)
                                         ttgccgttgccgtGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAgc                                                                                                      <  1:587613/150‑1 (MQ=255)
                                         ttgccgttgccgtGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAgc                                                                                                      <  2:138539/150‑1 (MQ=255)
                                         ttgccgttgccgtGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAgc                                                                                                      <  2:858187/150‑1 (MQ=255)
                                             cgttgccgtGGCTTTTGCTTCTAAATGCGCGCAATTATTGGTGTCTTAGAGTTCGGCTTCAGGAGGCCGAGCATTTGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGCCTCAgcgctt                                                                                                  <  1:672846/150‑3 (MQ=9)
                                                     tGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCgg                                                                                          <  1:1326127/150‑1 (MQ=40)
                                                     tGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCgg                                                                                          <  2:66373/150‑1 (MQ=38)
                                                     tGGCTTTTGCTTCTAAGTGCGCGATAGTCCAGGGGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCgg                                                                                          <  2:946280/150‑1 (MQ=255)
                                                       gCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                        <  1:639008/150‑1 (MQ=38)
                                                       gCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                        <  1:507742/150‑1 (MQ=38)
                                                       gCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                        >  2:578874/1‑150 (MQ=38)
                                                         ttttGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGATCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCACGTACTCCCGGGacc                                                                                      >  2:450784/1‑148 (MQ=255)
                                                           ttGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCaa                                                                                    >  2:520111/1‑150 (MQ=35)
                                                              cTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTAAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGAc                                                                                 <  1:1041299/150‑1 (MQ=37)
                                                               ttCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGt                                                                                                                                         <  2:1367028/94‑1 (MQ=39)
                                                               ttCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGt                                                                                                                                         >  1:1367028/1‑94 (MQ=39)
                                                               ttCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatgg                                                                                                                <  1:842634/118‑1 (MQ=37)
                                                               ttCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatgg                                                                                                                >  2:842634/1‑118 (MQ=38)
                                                               ttCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGAcc                                                                                >  2:722995/1‑150 (MQ=38)
                                                                      tGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                                                >  2:1091896/1‑96 (MQ=35)
                                                                      tGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                                                <  1:1091896/96‑1 (MQ=35)
                                                                       gcgcGATAGTTCAGGCGACGCAGATTGCGCCTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAg                                                                        >  1:1254341/1‑150 (MQ=38)
                                                                              aGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCACGGGAGCAAGACCGTAGTGAGTGAAcac                                                                 >  2:1149759/1‑148 (MQ=255)
                                                                                ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCg                                                               <  1:1149759/150‑1 (MQ=38)
                                                                                   aGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCaa                                                                                    >  2:279593/1‑126 (MQ=33)
                                                                                   aGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCaa                                                                                    <  1:279593/126‑1 (MQ=33)
                                                                                       gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                        >  1:1363493/1‑118 (MQ=33)
                                                                                       gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                        <  2:1363493/118‑1 (MQ=33)
                                                                                          gCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                       <  2:278363/69‑1 (MQ=11)
                                                                                          gCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                       >  1:278363/1‑69 (MQ=11)
                                                                                              aTTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAgg                                                 <  1:578874/150‑1 (MQ=38)
                                                                                                      ttCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAgg                                                 >  2:425416/1‑142 (MQ=39)
                                                                                                      ttCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAgg                                                 <  1:425416/142‑1 (MQ=39)
                                                                                                      ttCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgat                                         >  2:424036/1‑150 (MQ=39)
                                                                                                         aggaggCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                      <  1:493818/150‑1 (MQ=38)
                                                                                                           gaggCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATgg                                    >  1:908944/1‑150 (MQ=38)
                                                                                                            aggCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATggg                                   <  2:274840/150‑1 (MQ=38)
                                                                                                              gCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAg                                 <  2:674168/150‑1 (MQ=38)
                                                                                                              gCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAg                                 <  1:520111/150‑1 (MQ=39)
                                                                                                                          gCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGtttttt                     <  1:722995/150‑1 (MQ=35)
                                                                                                                               cGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGTTa                <  2:1354177/150‑1 (MQ=35)
                                                                                                                                    ggaggTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGTTACtttt           >  2:1034655/1‑150 (MQ=32)
                                                                                                                                        ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGa                                  >  2:257411/2‑123 (MQ=12)
                                                                                                                                        ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGa                                  <  1:257411/122‑1 (MQ=12)
                                                                                                                                           aCGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGtt                 <  2:254888/137‑1 (MQ=16)
                                                                                                                                           aCGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGtt                 >  1:254888/1‑137 (MQ=16)
                                                                                                                                             gggCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGTTACTTTTTGTCGCGtt  >  2:1336500/1‑150 (MQ=32)
                                                                                                                                             gggCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGTTACTTTTTGTCGCGtt  <  2:1254341/150‑1 (MQ=32)
                                                                                                                                                    |                                                                                                                                               
GCACAGTTCAAGGACCGTTGATCAGACTCGGTTCGTTGCCGTTGCCGTTGCCGTGGCTTTTGCTTCTAAGTGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCAGCGTTTTTTGGTTACTTTTTGTCGCGTT  >  NC_002947/2528896‑2529186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: