Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,392,927:1 +A 10.7% intergenic (+91/‑95) PP_1214 → / → ruvC transcriptional regulator/crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,392,9271.A10.7% 201.3 / 18.3 56intergenic (+91/‑95)PP_1214/ruvCtranscriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base . (17/33);  new base A (2/4);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGG  >  NC_002947/1392779‑1393076
                                                                                                                                                     |                                                                                                                                                     
gCAGAACGTCTATTCCAATGCCCAGATTTCCGACGACATCAGTGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGccc                                                                                                                                                      <  2:1378853/150‑1 (MQ=255)
   gAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCAtt                                                                                                                                                   >  1:2249326/1‑150 (MQ=255)
        tCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaaa                                                                                                                                              <  2:1050888/150‑1 (MQ=255)
           aTTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAAtgt                                                                                                                                           >  1:1151251/1‑150 (MQ=255)
           aTTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAAtg                                                                                                                                            >  1:451406/1‑150 (MQ=37)
                 aTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCa                                                                                                                                     <  2:1151251/150‑1 (MQ=255)
                       aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCt                                                                                                                               <  2:2175158/150‑1 (MQ=255)
                             atgcgcgagccggggtcgatacctagaatcagagtcataaagcccgtcgaaaaaatgatgcatctctacaaactgcGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGa                                                                                                                         >  1:2194299‑M1/77‑150 (MQ=255)
                               gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTAAAAATCTTGTTCATTACCTTTAGGAGc                                                                                                                       <  2:2040650/150‑1 (MQ=38)
                                acgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCg                                                                                                                      <  1:1169325/150‑1 (MQ=255)
                                acgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCg                                                                                                                      >  2:1847056/1‑150 (MQ=255)
                                acgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCg                                                                                                                      <  2:451406/150‑1 (MQ=18)
                                 cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                     <  1:2415181/150‑1 (MQ=255)
                                  gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCACCTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAATCCCGCTCCTACAGGTACTGCACCACTTTTGAAATGGGCGTAATACCTGTAGGAGCggg                                                                                                                    <  1:1260162/150‑1 (MQ=17)
                                  gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                                    <  2:2249326/150‑1 (MQ=255)
                                  gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                                    >  2:1733917/1‑150 (MQ=255)
                                  gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                                    <  1:2325547/150‑1 (MQ=255)
                                   agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGt                                                                                                                   <  2:1399202/150‑1 (MQ=255)
                                   agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGt                                                                                                                   <  1:1692242/150‑1 (MQ=255)
                                   agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGt                                                                                                                   <  1:1847056/150‑1 (MQ=255)
                                    gtctaTGGAGACCCTCGGCTAAGGCGACTTGGTTCCTTTTTGCAATTCTTAGCGTTTGATAGCGCCTGTGCCTGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGtt                                                                                                                  <  2:22975/146‑1 (MQ=25)
                                       cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAc                                                                                                               <  1:1733917/150‑1 (MQ=255)
                                        aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCgg                                                                                                                     <  2:637471/143‑1 (MQ=17)
                                        aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCgg                                                                                                                     >  1:637471/1‑143 (MQ=17)
                                           gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCcgc                                                                                                           <  1:411020/150‑1 (MQ=255)
                                               atacctagaatcagagtcataaagcccgtcgaaaaaatgatgcatctctacaaactgcGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAg                                                                                                                        >  2:2426848‑M1/59‑133 (MQ=255)
                                               atacctagaatcagagtcataaagcccgtcgaaaaaatgatgcatctctacaaactgcGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAg                                                                                                                        <  1:2426848‑M1/75‑1 (MQ=255)
                                                                            gCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaaa                                                                                                                                              <  1:2156526/82‑1 (MQ=255)
                                                                            gCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaaa                                                                                                                                              >  2:2156526/1‑82 (MQ=255)
                                                                                       gCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCACATTTGAAAtgggc                                                                                                                                         <  1:928508/76‑4 (MQ=21)
                                                                                       gCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCACATTTGAAAtgggc                                                                                                                                         >  2:928508/1‑73 (MQ=21)
                                                                                                             GCCCTTTCGCGGGTAAACCCGCTCCTACCGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGCGCCGTTTGTAGAGATGCATCATTTTTTCGACGGGGTTTATGTCTCTGattcctg                                         <  2:1783160‑M1/150‑4 (MQ=255)
                                                                                                                CTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTa                                      <  1:2346463‑M1/150‑1 (MQ=255)
                                                                                                                     GCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAc                                 <  2:342578‑M1/150‑1 (MQ=255)
                                                                                                                       GGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                   <  1:862517‑M1/82‑1 (MQ=255)
                                                                                                                       GGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                   >  2:862517‑M1/1‑82 (MQ=255)
                                                                                                                       GGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                               <  1:2064363‑M1/150‑1 (MQ=255)
                                                                                                                              CCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                        <  2:2190912‑M1/150‑1 (MQ=255)
                                                                                                                               CCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                       <  2:938042‑M1/150‑1 (MQ=255)
                                                                                                                               CCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                       >  1:342578‑M1/1‑150 (MQ=255)
                                                                                                                               CCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                       <  2:1341663‑M1/150‑1 (MQ=255)
                                                                                                                                CGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgcg                      <  2:2174782‑M1/150‑1 (MQ=255)
                                                                                                                                   TCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGAc                                                                >  1:887431‑M1/1‑105 (MQ=255)
                                                                                                                                   TCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGAc                                                                <  2:887431‑M1/105‑1 (MQ=255)
                                                                                                                                     CTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                 <  1:1419204‑M1/150‑1 (MQ=255)
                                                                                                                                     CTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                 <  2:2297205‑M1/150‑1 (MQ=255)
                                                                                                                                      TACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                <  1:779910‑M1/150‑1 (MQ=255)
                                                                                                                                            TATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTa          <  2:1563373‑M1/150‑1 (MQ=255)
                                                                                                                                             ATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTAt         >  2:2064363‑M1/1‑150 (MQ=255)
                                                                                                                                              ttACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                          >  2:1567959/1‑116 (MQ=255)
                                                                                                                                              ttACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                          <  1:1567959/116‑1 (MQ=255)
                                                                                                                                                aCGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTa                                                                                     <  2:519258/71‑1 (MQ=255)
                                                                                                                                                aCGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTa                                                                                     >  1:519258/1‑71 (MQ=255)
                                                                                                                                                  gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGa                                                     <  2:402055/101‑1 (MQ=255)
                                                                                                                                                  gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGa                                                     >  1:402055/1‑101 (MQ=255)
                                                                                                                                                  gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt    >  2:2279609/1‑150 (MQ=255)
                                                                                                                                                  gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt    >  1:210967/1‑150 (MQ=255)
                                                                                                                                                   cccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg   <  1:1563373/150‑1 (MQ=255)
                                                                                                                                                   cccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg   >  1:713798/1‑150 (MQ=255)
                                                                                                                                                    ccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg  >  1:95279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGG  >  NC_002947/1392779‑1393076

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