Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,946 C→G 7.1% G78G (GGC→GGG) ‡ moeA → molybdenum::molybdopterin ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,9460CG 7.1% 181.5 / 4.6 56G78G (GGC→GGG) ‡moeAmolybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base C (24/28);  new base G (2/2);  total (26/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 = 7.60e-01

CTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTC  >  NC_002947/2421798‑2422090
                                                                                                                                                     |                                                                                                                                                
cTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGc                                                                                                                                                  <  5:54029/149‑1 (MQ=255)
   cACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAg                                                                                                                                               >  6:333535/1‑149 (MQ=255)
    aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGc                                                                                                                                              >  4:39383/1‑149 (MQ=255)
    aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGc                                                                                                                                              <  8:89155/149‑1 (MQ=255)
     ccGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                             >  6:126558/1‑149 (MQ=255)
           cgcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                        <  8:157596/148‑1 (MQ=255)
            gcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCg                                                                                                                                      <  4:166988/149‑1 (MQ=255)
              gagaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                    <  6:290326/149‑1 (MQ=255)
              gagaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                    <  3:39383/149‑1 (MQ=255)
                taCCGCAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACCGCGCCATGGACGGGTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                                  <  8:269481/148‑1 (MQ=255)
                gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                                  <  5:243257/149‑1 (MQ=255)
                gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                                  <  7:91924/149‑1 (MQ=255)
                     aaaaCCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCa                                                                                                                                                <  8:239567/130‑1 (MQ=255)
                     aaaaCCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCa                                                                                                                                                >  7:239567/1‑130 (MQ=255)
                           tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                    <  5:5052/136‑1 (MQ=255)
                           tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                                    >  6:5052/1‑136 (MQ=255)
                           tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAtt                                                                                                                       >  3:110140/1‑149 (MQ=255)
                             ccGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAtttt                                                                                                                     >  3:64865/1‑149 (MQ=255)
                                   gctgAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCgg                                                                                                               <  5:252924/146‑1 (MQ=255)
                                         gCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATgc                                                                                                         >  2:155077/1‑149 (MQ=255)
                                           cGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTc                                                                                                                                   >  6:287703/1‑121 (MQ=38)
                                           cGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTc                                                                                                                                   <  5:287703/121‑1 (MQ=38)
                                               ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGTGCCCCa                                                                                                   <  2:110557/149‑1 (MQ=255)
                                               ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCa                                                                                                   <  4:64865/149‑1 (MQ=255)
                                                 ccGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGGCGGAACACCCCAATGGACGGTTATGCCTTGCGCATGGCTGCCTGGCAGGCCCAGCCGCTGCCGGTCAGCCAGCGTATTTACGCCGGCCATGCGCCCCAg                                                                                                  >  1:79733/1‑149 (MQ=255)
                                                     gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                             <  1:155077/149‑1 (MQ=255)
                                                         tGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGTCGGCCATGCGCCCCAGCCGCTGCaa                                                                                         >  3:72314/1‑149 (MQ=255)
                                                         tGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCaa                                                                                         >  7:46488/1‑149 (MQ=255)
                                                           gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGTCGGCCATGCGCCCCAGCCGCTGCAAc                                                                                        <  4:72314/148‑1 (MQ=255)
                                                                            ctcgctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGa                                                                      >  2:286994/1‑149 (MQ=255)
                                                                              cgctcgACCTGCCGCCCTGGCCGAACAGCACCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATc                                                                    <  2:79733/149‑1 (MQ=255)
                                                                                ctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                        >  7:147867/1‑79 (MQ=255)
                                                                                ctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                        <  8:147867/79‑1 (MQ=255)
                                                                                ctccACCTCCGGCCCGGGCCCCACCGGGCCTCGGCCGGTGAGGCCTTGCGCCTGGCGGCCTGGCAGGGCCAGCCGCTGCCGCGCAGCCAGCGGATTTTCCCCGGCCATGCGCCCGAGCCGCGGCATCCGGGTCCCTGCGCCCTGATCtt                                                                  <  4:266797/149‑1 (MQ=255)
                                                                                 tcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                  >  8:116815/1‑84 (MQ=21)
                                                                                 tcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                  <  7:116815/84‑1 (MQ=21)
                                                                                   gACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTc                                                                 <  7:17184/147‑1 (MQ=255)
                                                                                   gACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCAc                                                               <  5:26370/149‑1 (MQ=255)
                                                                                            cccTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGcc                                                      >  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                              <  1:286994/149‑1 (MQ=255)
                                                                                                                     ggTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGt                              <  2:225516/149‑1 (MQ=255)
                                                                                                                      gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                             <  5:126558/149‑1 (MQ=255)
                                                                                                                      gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                             <  1:78478/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)
                                                                                                                                                     |                                                                                                                                                
CTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTC  >  NC_002947/2421798‑2422090

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: