Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,355,891 T→C 7.1% L124P (CTG→CCG)  glpC → anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,355,8910TC 7.1% 166.6 / 3.4 56L124P (CTG→CCG) glpCanaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster
Reads supporting (aligned to +/- strand):  ref base T (21/31);  new base C (2/2);  total (23/33)
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.76e-01

TATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAAC  >  NC_000913/2355745‑2356036
                                                                                                                                                  |                                                                                                                                                 
tatCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAgctgc                                                                                                                                                 >  2:269180/1‑149 (MQ=255)
    cAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGa                                                                                                                                             <  6:66075/149‑1 (MQ=255)
    cAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGa                                                                                                                                             <  2:98847/149‑1 (MQ=255)
      gcgcgcgcGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATg                                                                                                                                           <  4:144137/149‑1 (MQ=255)
       cgcgcgctGGCGAAATATGACCACACGCGCCCGTCGCTGCGTAATTTTTTGTGGAGTCATACCGACCTTCGGGGTAGCGTTTCCACGCCGGTCTCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATgc                                                                                                                                          <  2:197741/149‑1 (MQ=255)
        gcgcgcgGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATgcg                                                                                                                                         <  1:300761/149‑1 (MQ=255)
                   aaaTATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTaaaaa                                                                                                                               >  6:206464/1‑148 (MQ=255)
                   aaaTATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAAt                                                                                                                              <  5:301489/149‑1 (MQ=255)
                   aaaTATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAAt                                                                                                                              <  1:269180/149‑1 (MQ=255)
                        tGACACCCCGCTCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATAGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGGGGCGTTAAAAATCGat                                                                                                                          <  7:246680/148‑1 (MQ=255)
                          acacCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGatca                                                                                                                        <  3:215674/148‑1 (MQ=255)
                           caccccGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCCTACCGAACTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGatcatc                                                                                                                      <  1:97170/149‑1 (MQ=255)
                            accacGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATcg                                                                                                                     >  2:226992/1‑149 (MQ=255)
                             ccacGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATcgc                                                                                                                    >  2:307398/1‑149 (MQ=255)
                               acGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATcgccg                                                                                                                  >  7:265339/1‑149 (MQ=255)
                                  cgcCCGGCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCAc                                                                                                               <  1:276319/149‑1 (MQ=255)
                                       gTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACTCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTAc                                                                                                          <  6:43902/149‑1 (MQ=255)
                                            tGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGaaa                                                                                                     <  7:163837/149‑1 (MQ=255)
                                            tGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGaaa                                                                                                     >  3:186620/1‑149 (MQ=255)
                                            agCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGaaa                                                                                                     >  1:175746/2‑149 (MQ=255)
                                              cGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGaaa                                                                                                     <  4:186620/147‑1 (MQ=255)
                                                      tgtgTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGGCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTc                                                                                            <  4:254930/148‑1 (MQ=255)
                                                          ttGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCAc                                                                                       >  8:104367/1‑149 (MQ=255)
                                                          ttGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCAc                                                                                       >  4:87745/1‑149 (MQ=255)
                                                          ttGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCAc                                                                                       >  8:280502/1‑149 (MQ=255)
                                                                 aTACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGt                                                                                <  4:154181/149‑1 (MQ=255)
                                                                     cGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCt                                                                            >  1:273948/1‑149 (MQ=255)
                                                                      gACCTGATGGGTAGCGTTTCCACGCCGTTCTCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTg                                                                           <  5:84579/149‑1 (MQ=255)
                                                                        ccTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGt                                                                         <  8:183862/149‑1 (MQ=255)
                                                                          tGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTAt                                                                       >  6:191329/1‑149 (MQ=255)
                                                                            aTGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTAt                                                                       <  5:206464/147‑1 (MQ=255)
                                                                               ggTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAg                                                                  >  8:38032/1‑149 (MQ=255)
                                                                                   gCGTTTCCACGCCGTTCGCACCACTCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCTGCTGCTTGATGCGGCGGTAAAACTCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTg                                                              <  7:208145/149‑1 (MQ=255)
                                                                                   gCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTg                                                              <  7:208144/149‑1 (MQ=255)
                                                                                    cGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTgg                                                             <  2:243474/149‑1 (MQ=255)
                                                                                     gTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAACCCGGTGCGGCCGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTCCCGAAATACTCCTTCGGCGCGTTCCGTCGCTGGTATCGCCGCGTggc                                                            >  5:336798/1‑149 (MQ=255)
                                                                                        tCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTggcgg                                                          <  4:28929/148‑1 (MQ=255)
                                                                                               cGTTCGCACCAATCGTCAACACCGCTACCTCGCTGATACCGTTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCAAGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaa                                                  <  2:175746/149‑1 (MQ=255)
                                                                                                gTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaag                                                 >  6:99000/1‑149 (MQ=255)
                                                                                                gTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaag                                                 <  1:60540/149‑1 (MQ=255)
                                                                                                    gCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCCGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGcac                                              >  6:131088/1‑148 (MQ=255)
                                                                                                     cACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACaa                                            <  3:87745/149‑1 (MQ=255)
                                                                                                              tCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGAc                                   >  4:75346/1‑149 (MQ=255)
                                                                                                                 acacCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAg                                >  8:24933/1‑149 (MQ=255)
                                                                                                                     cGCTACCTCGCTGAAACCGGTGCGGCAGCCGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCg                            >  2:139105/1‑149 (MQ=255)
                                                                                                                      gCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCAAGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGc                           >  6:130610/1‑149 (MQ=255)
                                                                                                                             cGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGctttcttt                    >  5:257776/1‑149 (MQ=255)
                                                                                                                              gCTGAAACCGGTGCGGCAGCCGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCtttctttc                   <  1:139105/149‑1 (MQ=255)
                                                                                                                                  aaaCCGGTGCGGCAGCCGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCAc                 <  5:131088/147‑1 (MQ=255)
                                                                                                                                    aCCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGgct             <  7:280502/149‑1 (MQ=255)
                                                                                                                                    aCCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGgct             <  7:38032/149‑1 (MQ=255)
                                                                                                                                         tGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGGCGCTTTCTTTCACGgctgctgc        <  1:282728/149‑3 (MQ=255)
                                                                                                                                           cGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCg       <  7:24933/148‑1 (MQ=255)
                                                                                                                                            ggCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGt      >  7:176906/1‑148 (MQ=255)
                                                                                                                                             gCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGtt     >  7:165664/1‑148 (MQ=255)
                                                                                                                                               aGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTg                                                                           >  7:206872/1‑76 (MQ=255)
                                                                                                                                               aGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTg                                                                           <  8:206872/76‑1 (MQ=255)
                                                                                                                                                gctgctTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGACGCTTTCTTTCACGGCTGCTTCGTTAAc  >  3:104414/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAAC  >  NC_000913/2355745‑2356036

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: