Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,486,819 C→G 7.0% intergenic (+52/‑132) spuC‑I → / → PP_2181 polyamine:pyruvate transaminase/Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,486,8190CG 7.0% 177.5 / 3.7 57intergenic (+52/‑132)spuC‑I/PP_2181polyamine:pyruvate transaminase/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base G (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GCCACCGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGA  >  NC_002947/2486677‑2486964
                                                                                                                                              |                                                                                                                                                 
gcCACCGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGcgcg                                                                                                                                             <  5:227/149‑1 (MQ=255)
      gCTGTTCATGAGTATCGAAGAGAGCGACGTTCTGGGGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGATGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGATc                                                                                                                                       <  8:146536/149‑1 (MQ=255)
           tGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGGCCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCGAGCACCGCGCAGGCGGTGCTTGATct                                                                                                                                  <  3:259896/149‑2 (MQ=38)
               cAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTgc                                                                                                                              >  6:400229/1‑149 (MQ=255)
                    tCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTgccgcc                                                                                                                          <  4:375795/148‑1 (MQ=255)
                     cGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCAt                                                                                                                        >  8:353431/1‑149 (MQ=255)
                         gagaTCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTa                                                                                                                    <  5:15985/149‑1 (MQ=255)
                              cgacgaACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACAttt                                                                                                               >  5:396651/1‑149 (MQ=255)
                               gacgaACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTc                                                                                                              <  3:329636/149‑1 (MQ=255)
                                acgaACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCa                                                                                                             <  3:197996/149‑1 (MQ=255)
                                 cgaACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCAc                                                                                                            <  7:380855/149‑1 (MQ=255)
                                        gtggAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGCTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAgg                                                                                                      >  4:89602/1‑148 (MQ=255)
                                              aaaGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTc                                                                                               >  7:387251/1‑149 (MQ=255)
                                                   ccGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCGAGCACCGCGCAGGCGGTGCTTGATct                                                                                                                                  <  2:390520/109‑2 (MQ=38)
                                                   ccGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCa                                                                                           >  3:379909/1‑148 (MQ=255)
                                                   ccGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACAATTCACGTCAGGCACCTTCCGCa                                                                                           >  5:188352/1‑148 (MQ=255)
                                                      cAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGc                                                                                       <  6:35510/149‑1 (MQ=255)
                                                           gCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCaaaat                                                                                  >  4:271070/1‑145 (MQ=255)
                                                                gacctgacctATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCttt                                                                             <  4:379909/149‑1 (MQ=255)
                                                                   ctgacctATGAGGCTGTTCGGTAACAGCGGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCttt                                                                             <  2:46669/146‑1 (MQ=255)
                                                                   ctgacctATGAGGCTGTTCGGTAACAGCGGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCttt                                                                             >  1:46669/1‑146 (MQ=255)
                                                                     gaactATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTATCACCACCTGTGAGATCAAGCACCGCATGCGCGGTGCACGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATATTCCTTTGTGa                                                                         >  6:270401/4‑148 (MQ=255)
                                                                      acctATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCAAAATGCTTtgcgg                                                                         <  3:271070/147‑4 (MQ=255)
                                                                         tATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCCAAGCATTTTGCTTTGTGAAATg                                                                     >  5:319912/1‑148 (MQ=255)
                                                                          aTGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCGCCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGaa                                                                   >  3:30911/1‑149 (MQ=255)
                                                                                 tGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTa                                                            >  5:183481/1‑149 (MQ=255)
                                                                                 tGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTa                                                            >  2:106601/1‑149 (MQ=255)
                                                                                     cGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACttt                                                        >  2:437187/1‑149 (MQ=255)
                                                                                       gTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTc                                                      >  8:89865/1‑149 (MQ=255)
                                                                                       gTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTc                                                      >  4:199550/1‑149 (MQ=255)
                                                                                          aCAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTgg                                                   <  8:387251/149‑1 (MQ=255)
                                                                                            agcagAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTgg                                                   <  1:109670/147‑1 (MQ=255)
                                                                                            agcagAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTgg                                                   >  2:109670/1‑147 (MQ=255)
                                                                                                     ccTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGAACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTc                                        <  4:35726/149‑1 (MQ=255)
                                                                                                      cTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTcc                                       <  5:400229/149‑1 (MQ=255)
                                                                                                         gTTAGCTGTTCTCAGCACCTGTGAGATCAAGCGCCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGt                                    <  4:30911/149‑1 (MQ=255)
                                                                                                                 ttCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGt                            >  3:162102/1‑149 (MQ=255)
                                                                                                                  tctcAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTc                           >  5:105270/1‑149 (MQ=255)
                                                                                                                    tcAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTcgc                         <  3:345722/149‑1 (MQ=255)
                                                                                                                    tcAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTcgc                         >  2:367531/1‑149 (MQ=255)
                                                                                                                      aGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGa                       >  3:275719/1‑149 (MQ=255)
                                                                                                                         aCCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGAtt                     >  8:284208/1‑148 (MQ=255)
                                                                                                                                 gaTCGAGCACCGCGCAGGCGGTGCTTGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATg             >  3:49874/1‑148 (MQ=255)
                                                                                                                                 gaTCGAGCACCGCGCAGGCGGTGCTTGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATg             >  1:218696/1‑148 (MQ=255)
                                                                                                                                 gaTCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGa            <  8:25595/149‑1 (MQ=255)
                                                                                                                                   tCAAGCACCGCCTGCGCTGTGCTCTATCGTTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTTTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGaaa          <  3:89602/149‑1 (MQ=255)
                                                                                                                                   tCAACCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGaaa          <  7:362240/149‑1 (MQ=255)
                                                                                                                                    cAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAg         >  1:341352/1‑149 (MQ=255)
                                                                                                                                     aaGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCt                                                     <  7:412246/104‑1 (MQ=255)
                                                                                                                                     aaGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCt                                                     >  8:412246/1‑104 (MQ=255)
                                                                                                                                     aaGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGt        >  3:417524/1‑149 (MQ=255)
                                                                                                                                      aGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTg       >  1:401267/1‑149 (MQ=255)
                                                                                                                                      aGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTg       <  7:263278/149‑1 (MQ=255)
                                                                                                                                       gCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAAGCCTGTACTTTTCTGGGCATCCATTTCCCGGCTGAACGTCGCGATTCCTCCAGGAGAGTGc      >  8:157176/1‑149 (MQ=255)
                                                                                                                                         acggCCTGCTCGGTGGTCGATCTGTGCCGCCAACGTAAATTTCTCGTCAGGAACCTTCCGCAATGCATTTTGCTTTGTGAATTGAATCCTGTACTTTTCTTGGCATCCATTTCCCGTCTGAAAGTCGCTATTCAACCATGAAAGTGCAc    <  5:270401/146‑1 (MQ=255)
                                                                                                                                         aCCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATACATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCAc    >  7:379247/1‑149 (MQ=255)
                                                                                                                                           cGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACga  <  7:89865/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCCACCGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGACCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGA  >  NC_002947/2486677‑2486964

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