Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 977,876 T→A 8.4% intergenic (+39/‑141) cysE → / → iscR serine acetyltransferase/DNA‑binding transcriptional regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947977,8760TA 8.4% 169.5 / 5.2 57intergenic (+39/‑141)cysE/iscRserine acetyltransferase/DNA‑binding transcriptional regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base A (2/3);  total (27/30)
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.57e-01

GCGCTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCGACTG  >  NC_002947/977727‑978025
                                                                                                                                                     |                                                                                                                                                     
gcgcTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCgtgt                                                                                                                                                       <  1:1102264/150‑1 (MQ=255)
       ccAAGATGGGCAGTGACTATTGTGCCAAGGAAATGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACg                                                                                                                                                <  2:1123071/150‑1 (MQ=255)
        cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGc                                                                                                                                               <  1:461755/150‑1 (MQ=255)
        cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGAGGGGCGTGAGCACACGc                                                                                                                                               >  1:338689/1‑150 (MQ=255)
         aaGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGcc                                                                                                                                              >  1:1085077/1‑150 (MQ=255)
             tGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACGccccgc                                                                                                                                          <  2:30288/150‑3 (MQ=255)
              gggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGccc                                                                                                                                             <  2:499258/146‑1 (MQ=255)
              gggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGccc                                                                                                                                             >  1:499258/1‑146 (MQ=255)
               ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACGccccgcat                                                                                                                                        <  1:666447/150‑5 (MQ=255)
               ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGAGGGGCGTGAGCACACGccccgcat                                                                                                                                        <  1:1145204/150‑5 (MQ=255)
                                cAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACATTCGcgccc                                                                                                                       >  1:899644/1‑150 (MQ=255)
                                 aaGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcccg                                                                                                                      >  1:892683/1‑150 (MQ=255)
                                   ggAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGc                                                                                                                                    >  2:889620/1‑134 (MQ=255)
                                   ggAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGc                                                                                                                                    <  1:889620/134‑1 (MQ=255)
                                   ggAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcccgcc                                                                                                                    >  2:617769/1‑150 (MQ=255)
                                       cTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTccc                                                                                                                <  1:954336/150‑1 (MQ=255)
                                             gCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGc                                                                                                          >  1:545948/1‑150 (MQ=255)
                                                ttGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCAACACTGAAGTGATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAAt                                                                                                       >  2:41051/1‑150 (MQ=255)
                                                    cGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCg                                                                                                   >  1:125102/1‑150 (MQ=255)
                                                     ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGa                                                                                                  >  1:124216/1‑150 (MQ=255)
                                                            gacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGc                                                                                           <  2:125102/150‑1 (MQ=255)
                                                             acTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGGCCGCCTCCCTGATGCAATCCCGACTAAAGCa                                                                                          <  2:740366/150‑1 (MQ=255)
                                                              cTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATg                                                                                                           >  2:277839/1‑132 (MQ=255)
                                                              cTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATg                                                                                                           <  1:277839/132‑1 (MQ=255)
                                                               ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGt                                                                                        <  2:1085077/150‑1 (MQ=255)
                                                               ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGt                                                                                        >  2:124188/1‑150 (MQ=255)
                                                                 cAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAg                                                                                      <  1:145186/150‑1 (MQ=255)
                                                                   gCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGt                                                                                    >  1:213889/1‑150 (MQ=255)
                                                                        gCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTCGGTCTAAt                                                                               >  2:581090/1‑150 (MQ=255)
                                                                              cAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa                                                                         <  2:387280/150‑1 (MQ=255)
                                                                                         gcgcGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTc                                                              <  2:213889/150‑1 (MQ=255)
                                                                                          cgcgACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCg                                                             <  2:124216/150‑1 (MQ=255)
                                                                                          cgcgACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGCGTCTTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCg                                                             <  2:635566/150‑1 (MQ=255)
                                                                                                acTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATccc                                                                                                    <  1:748250/105‑1 (MQ=255)
                                                                                                acTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATccc                                                                                                    >  2:748250/1‑105 (MQ=255)
                                                                                                acTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATc                                                       <  1:594008/150‑1 (MQ=255)
                                                                                                  tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTaaa                                                                                             >  2:231428/1‑110 (MQ=255)
                                                                                                  tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTaaa                                                                                             <  1:231428/110‑1 (MQ=255)
                                                                                                  tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTaaa                                                                                             <  1:16963/110‑1 (MQ=255)
                                                                                                  tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTaaa                                                                                             >  2:16963/1‑110 (MQ=255)
                                                                                                  tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGc                                                     >  1:154167/1‑150 (MQ=255)
                                                                                                   cAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCa                                                    <  1:568028/150‑1 (MQ=255)
                                                                                                      ccGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAAt                                                                                                       <  1:714327/96‑1 (MQ=255)
                                                                                                      ccGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAAt                                                                                                       >  2:714327/1‑96 (MQ=255)
                                                                                                        gCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTc                                                              >  2:894383/1‑135 (MQ=255)
                                                                                                        gCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTc                                                              <  1:894383/135‑1 (MQ=255)
                                                                                                          aTTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGc                                             <  1:124188/150‑1 (MQ=255)
                                                                                                               tCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGcacaca                                        >  2:1141546/1‑150 (MQ=255)
                                                                                                                   gccgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCt                                    <  2:892683/150‑1 (MQ=255)
                                                                                                                       cACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGaaa                                <  1:1141546/150‑1 (MQ=255)
                                                                                                                              gTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGtt                                                                           >  2:553987/1‑100 (MQ=255)
                                                                                                                              gTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGtt                                                                           <  1:553987/100‑1 (MQ=255)
                                                                                                                               tCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTg                        >  2:426653/1‑150 (MQ=255)
                                                                                                                                         aTGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATa              <  1:617769/150‑1 (MQ=255)
                                                                                                                                          tGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCCCACATCCTGAAACTCCCTTGGTACCCAATAg             >  2:395711/1‑150 (MQ=255)
                                                                                                                                          tGCGGGGAGTGAGCTAAAGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAg             >  2:897513/1‑150 (MQ=255)
                                                                                                                                              gggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCa          <  1:234484/149‑1 (MQ=255)
                                                                                                                                              gggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCa          >  2:234484/1‑149 (MQ=255)
                                                                                                                                                     tgCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACGCCCTTGGTACCCAATAGCCATGCgactg  <  1:1125152/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCGACTG  >  NC_002947/977727‑978025

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: ‑