Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,950 C→G 7.7% P80A (CCG→GCG) ‡ moeA → molybdenum::molybdopterin ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,9500CG 7.7% 162.6 / 4.5 52P80A (CCG→GCG) ‡moeAmolybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base G (2/2);  total (23/29)
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.23e-01

ACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATGG  >  NC_002947/2421802‑2422098
                                                                                                                                                     |                                                                                                                                                    
aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGc                                                                                                                                                      <  8:89155/149‑1 (MQ=255)
aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGc                                                                                                                                                      >  4:39383/1‑149 (MQ=255)
 ccGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                                     >  6:126558/1‑149 (MQ=255)
       cgcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                <  8:157596/148‑1 (MQ=255)
        gcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCg                                                                                                                                              <  4:166988/149‑1 (MQ=255)
          gagaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                            <  6:290326/149‑1 (MQ=255)
          gagaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                            <  3:39383/149‑1 (MQ=255)
            taCCGCAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACCGCGCCATGGACGGGTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                                          <  8:269481/148‑1 (MQ=255)
            gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                                          <  7:91924/149‑1 (MQ=255)
            gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                                          <  5:243257/149‑1 (MQ=255)
                       tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                            >  6:5052/1‑136 (MQ=255)
                       tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                            <  5:5052/136‑1 (MQ=255)
                       tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAtt                                                                                                                               >  3:110140/1‑149 (MQ=255)
                         ccGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAtttt                                                                                                                             >  3:64865/1‑149 (MQ=255)
                               gctgAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCgg                                                                                                                       <  5:252924/146‑1 (MQ=255)
                                     gCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATgc                                                                                                                 >  2:155077/1‑149 (MQ=255)
                                       cGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTc                                                                                                                                           >  6:287703/1‑121 (MQ=38)
                                       cGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTc                                                                                                                                           <  5:287703/121‑1 (MQ=38)
                                           ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGTGCCCCa                                                                                                           <  2:110557/149‑1 (MQ=255)
                                           ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCa                                                                                                           <  4:64865/149‑1 (MQ=255)
                                             ccGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGGCGGAACACCCCAATGGACGGTTATGCCTTGCGCATGGCTGCCTGGCAGGCCCAGCCGCTGCCGGTCAGCCAGCGTATTTACGCCGGCCATGCGCCCCAg                                                                                                          >  1:79733/1‑149 (MQ=255)
                                                 gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                     <  1:155077/149‑1 (MQ=255)
                                                     tGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGTCGGCCATGCGCCCCAGCCGCTGCaa                                                                                                 >  3:72314/1‑149 (MQ=255)
                                                     tGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCaa                                                                                                 >  7:46488/1‑149 (MQ=255)
                                                       gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGTCGGCCATGCGCCCCAGCCGCTGCAAc                                                                                                <  4:72314/148‑1 (MQ=255)
                                                                        ctcgctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGa                                                                              >  2:286994/1‑149 (MQ=255)
                                                                          cgctcgACCTGCCGCCCT‑GGCCGAACAGCACCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATc                                                                            <  2:79733/149‑1 (MQ=255)
                                                                            ctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                <  8:147867/79‑1 (MQ=255)
                                                                            ctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                >  7:147867/1‑79 (MQ=255)
                                                                            ctccACCTCCGGCCCG‑GGCCCCACCGGGCCTCGGCCGGTGAGGCCTTGCGCCTGGCGGCCTGGCAGGGCCAGCCGCTGCCGCGCAGCCAGCGGATTTTCCCCGGCCATGCGCCCGAGCCGCGGCATCCGGGTCCCTGCGCCCTGATCtt                                                                          <  4:266797/149‑1 (MQ=255)
                                                                             tcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                          >  8:116815/1‑84 (MQ=21)
                                                                             tcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                          <  7:116815/84‑1 (MQ=21)
                                                                               gACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTc                                                                         <  7:17184/147‑1 (MQ=255)
                                                                               gACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCAc                                                                       <  5:26370/149‑1 (MQ=255)
                                                                                        cccT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGcc                                                              >  3:248065/1‑149 (MQ=255)
                                                                                               ccGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGcc                                                        >  4:36050/1‑149 (MQ=255)
                                                                                                cGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCg                                                       >  2:232683/1‑149 (MQ=255)
                                                                                                cGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCACTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCg                                                       <  5:333535/149‑1 (MQ=255)
                                                                                                        gcCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGc                                               <  4:248065/149‑1 (MQ=255)
                                                                                                                 ggTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGt                                      <  2:225516/149‑1 (MQ=255)
                                                                                                                 ggTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGt                                      <  1:286994/149‑1 (MQ=255)
                                                                                                                  gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                                     <  1:78478/149‑1 (MQ=255)
                                                                                                                  gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                                     <  5:126558/149‑1 (MQ=255)
                                                                                                                            gcgcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCa                             <  4:110140/147‑1 (MQ=255)
                                                                                                                              gcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGa                          >  5:110252/1‑148 (MQ=255)
                                                                                                                                   gCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAacac                    >  7:3558/1‑149 (MQ=255)
                                                                                                                                     tGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc                   >  3:212625/1‑148 (MQ=255)
                                                                                                                                     tGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc                   >  2:77328/1‑148 (MQ=255)
                                                                                                                                       aCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCg                  >  4:158398/1‑147 (MQ=255)
                                                                                                                                         tGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTTTTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAg               >  3:208538/1‑148 (MQ=255)
                                                                                                                                         tGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAg               >  6:246213/1‑148 (MQ=255)
                                                                                                                                           gCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTAc            <  8:131603/149‑1 (MQ=255)
                                                                                                                                             aGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTc          >  4:294720/1‑149 (MQ=255)
                                                                                                                                                   aGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGAt    >  7:208063/1‑149 (MQ=255)
                                                                                                                                                    gCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATg   <  8:208063/149‑1 (MQ=255)
                                                                                                                                                     ccGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATgg  >  1:237409/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATGG  >  NC_002947/2421802‑2422098

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