Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,675,534 T→C 10.5% D136G (GAC→GGC) ‡ yceJ ← iron‑regulated cytochrome b561

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,675,5340TC10.5% 148.8 / 7.1 57D136G (GAC→GGC) ‡yceJiron‑regulated cytochrome b561
Reads supporting (aligned to +/- strand):  ref base T (29/22);  new base C (4/2);  total (33/24)
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.84e-01

TCAAGCTTTGCGGCCCAGCATGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGC  >  NC_002947/5675389‑5675681
                                                                                                                                                 |                                                                                                                                                   
tCAAGCTTTGCGGCCCAGCATGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTcag                                                                                                                                                  >  7:270670/1‑149 (MQ=255)
  aaGCTTTGCGGCCCAGCATGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTcagg                                                                                                                                                 >  8:166657/1‑148 (MQ=255)
   aGCTTTGCGGCCCAGCATGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGt                                                                                                                                                     >  8:33712/1‑143 (MQ=255)
   aGCTTTGCGGCCCAGCATGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGt                                                                                                                                                     <  7:33712/143‑1 (MQ=255)
                   aTGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGc                                                                                                                               >  6:18686/1‑149 (MQ=255)
                    tGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGcc                                                                                                                              <  7:132160/149‑1 (MQ=255)
                    tGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGcc                                                                                                                              <  3:158051/149‑1 (MQ=255)
                      cGAGTCAGGGTCGCGTCACGGTCGATGAAATGATCCTTGAGCCCCGCAAGGGCATGCAGCACGGCAAAAATCACCCGCCCCCAAGCCAGCCACAAATCAACCCCCCCCGCCACAACTGCCTGGCCAGGCAGCTCGCTGATCACTGCCgg                                                                                                                            >  6:37997/1‑149 (MQ=255)
                          tCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCAc                                                                                                                         >  7:164195/1‑148 (MQ=255)
                          tCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGACcaggcagatgtg                                                                                                                                         <  5:334635/132‑6 (MQ=255)
                          tCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGACcaggcagatgtg                                                                                                                                         >  6:334635/1‑127 (MQ=255)
                            aGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTc                                                                                                                      <  5:18686/149‑1 (MQ=255)
                            aGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTc                                                                                                                      <  8:164195/149‑1 (MQ=255)
                                tCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTcagg                                                                                                                                                 >  8:193750/1‑118 (MQ=255)
                                tCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTcagg                                                                                                                                                 <  7:193750/118‑1 (MQ=255)
                                    gTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGc                                                                                                              <  7:166657/149‑1 (MQ=255)
                                        cGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGaa                                                                                                          <  8:270670/149‑1 (MQ=255)
                                         ggTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGc                                                                                                              <  3:103473/144‑1 (MQ=255)
                                         ggTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGc                                                                                                              >  4:103473/1‑144 (MQ=255)
                                                     gatgCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGcc                                                                                              >  6:284252/1‑148 (MQ=255)
                                                        gCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAAt                                                                                                 <  5:30723/142‑1 (MQ=255)
                                                        gCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAAt                                                                                                 >  6:30723/1‑142 (MQ=255)
                                                               cgcgGCAAGGGGATGCAGCACGGCAAAAATCACCAGCCCCCACGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCTCTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTcgg                                                                                   <  5:37997/149‑1 (MQ=255)
                                                                gcgGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGACcaggcagat                                                                                                                                            <  4:40341/91‑3 (MQ=255)
                                                                gcgGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGACcaggcagat                                                                                                                                            >  3:40341/1‑89 (MQ=255)
                                                                  ggCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTcggcgg                                                                                <  3:274387/149‑1 (MQ=255)
                                                                          cATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGGCGGCGGTGGAAATc                                                                        <  5:284252/149‑1 (MQ=255)
                                                                            tGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAg                                                                      <  6:349684/149‑1 (MQ=255)
                                                                            tGCAGCACGGCAAAAATCACAAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAAAAGGCCGAACACCGGAATGCCGACGCCGACGGCGGTGGAAAACAg                                                                      >  7:278565/1‑149 (MQ=255)
                                                                                          aaTCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAAcc                                                         >  6:13116/1‑148 (MQ=255)
                                                                                                cAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGGCGTcggcattccg                                                                            <  6:57748/123‑7 (MQ=255)
                                                                                                cAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGGCGTcggcattccg                                                                            >  5:57748/1‑117 (MQ=255)
                                                                                                     cccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCaaa                                              >  4:338004/1‑148 (MQ=255)
                                                                                                     cccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCaaa                                              >  1:5546/1‑148 (MQ=255)
                                                                                                     cccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAAc                                             <  3:338004/149‑1 (MQ=255)
                                                                                                      ccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACa                                            >  3:287681/1‑149 (MQ=255)
                                                                                                        aaGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGc                                          <  2:5546/149‑1 (MQ=255)
                                                                                                          gccagccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGcc                                         >  2:328978/1‑148 (MQ=255)
                                                                                                               ccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAgg                                    >  6:82732/1‑148 (MQ=255)
                                                                                                               ccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAgg                                    >  4:21914/1‑148 (MQ=255)
                                                                                                               ccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAgg                                    >  7:38835/1‑148 (MQ=255)
                                                                                                               ccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAgg                                    >  7:33584/1‑148 (MQ=255)
                                                                                                                   aaaTGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAg                               >  2:180589/1‑149 (MQ=255)
                                                                                                                    aaTGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAgg                              >  3:338730/1‑149 (MQ=255)
                                                                                                                      tGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGc                            >  1:232892/1‑149 (MQ=255)
                                                                                                                      tGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGc                            <  5:13116/149‑1 (MQ=255)
                                                                                                                       gAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGc                            >  2:314956/1‑148 (MQ=255)
                                                                                                                           cacaCCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACAGCAAACAGCCCAAGGTACAGGGCCAGGt                       >  7:285070/1‑149 (MQ=255)
                                                                                                                               cccGCGAAATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGcc                   <  7:264647/149‑1 (MQ=255)
                                                                                                                                      gaccTGCCTGACCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATTACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTg            >  2:173697/4‑149 (MQ=255)
                                                                                                                                      caTCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTg            <  1:314956/149‑1 (MQ=255)
                                                                                                                                       aTCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTgg           >  4:51788/1‑149 (MQ=255)
                                                                                                                                       aTCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTgg           >  7:1872/1‑149 (MQ=255)
                                                                                                                                           gCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACCGCCCAAGGTACAGGGCCAGGAGGCCCAGCTTGgccgc       >  6:20777/1‑149 (MQ=255)
                                                                                                                                           gCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGgccgc       >  4:63856/1‑149 (MQ=255)
                                                                                                                                           gCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGgccgc       >  4:89988/1‑149 (MQ=255)
                                                                                                                                            ccTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGgc          <  4:287681/145‑1 (MQ=255)
                                                                                                                                             cTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCa     <  3:89988/149‑1 (MQ=255)
                                                                                                                                                gTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGc  <  2:232892/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCAAGCTTTGCGGCCCAGCATGCGAGTCAGGGTCGCGTCACGGTCGATGAAATGATGCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGC  >  NC_002947/5675389‑5675681

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: