Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,631,689 Δ1 bp intergenic (+251/‑579) emrB → / → RPPX_20645 multidrug resistance protein B/oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,631,6890G.98.3% 190.2 / ‑2.0 58intergenic (+251/‑579)emrB/RPPX_20645multidrug resistance protein B/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (27/30);  minor base A (1/0);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTG  >  CP009974/4631543‑4631712
                                                                                                                                                  |                       
gcgcAGGGCGTGCGAAGGGATGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACGCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCg                    <  2:5772399/151‑1 (MQ=11)
 cgcaGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGt                   >  1:2492548/4‑151 (MQ=255)
  gcaGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCg                    >  1:531680/3‑149 (MQ=255)
  gcaGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCg                    <  2:531680/147‑1 (MQ=255)
  gcaGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGtt                  >  2:2616366/3‑151 (MQ=255)
  gTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGtt                  >  1:2555100/1‑151 (MQ=255)
  gTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGtt                  <  1:5023914/151‑1 (MQ=255)
   tAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGttt                 <  1:7091219/151‑1 (MQ=255)
   tAGGGCTGGCGAAGGGAGGGGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGttt                 <  2:3059547/151‑1 (MQ=255)
   caGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGttt                 <  1:4522464/150‑1 (MQ=255)
    aGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGtttt                >  2:698933/1‑151 (MQ=255)
    aGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGtttt                <  1:7007882/151‑1 (MQ=255)
     gggCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGttttt               >  2:2053917/1‑151 (MQ=255)
      ggCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGtttttt              >  1:3430205/1‑151 (MQ=255)
       ggttGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGAGGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTg             <  1:3474702/149‑1 (MQ=255)
       gCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTg             <  2:4815568/151‑1 (MQ=255)
       gCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTg             <  1:4422556/151‑1 (MQ=255)
        cTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTgg            >  2:6214200/1‑151 (MQ=255)
         ttGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCTCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGtttttgtgt           >  1:943113/1‑147 (MQ=255)
         ttGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGt           <  2:5725811/151‑1 (MQ=16)
         ttGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGt           <  1:6393419/151‑1 (MQ=16)
         ttGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGt           <  1:5898639/151‑1 (MQ=16)
         ttGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGt           >  1:618186/1‑151 (MQ=255)
         ttGCGAAGGGAGGGGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGt           <  2:3750766/151‑1 (MQ=255)
          tGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGtt          >  1:1134505/1‑151 (MQ=255)
          tGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGtt          <  2:7066548/151‑1 (MQ=16)
           gCGAAGGGAGGTGACGGGCATGGATGCCCGTCTAGCGCTGAGGCTCCATGGATGGGGCCTGCAGCGCGTATTCCTGGGAGCATGACCGTAGCGAGGCGATCCCGGAGCGTAGCGCCGGGGCCGGATGATAGGAGCCAGCGTCTTCTGGTTa         >  1:3194712/1‑151 (MQ=11)
           gCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTa         <  1:6890355/151‑1 (MQ=255)
           gCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTa         <  1:4647701/151‑1 (MQ=255)
           gCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTa         <  1:4362561/151‑1 (MQ=255)
           gCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGACAGGAGCCAGCGTTTTTTggctc         >  1:5197030/1‑148 (MQ=255)
           gCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGAAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCACGCGATTTTTGGTTa         >  1:3192298/1‑151 (MQ=11)
            ggAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTAc        <  2:2442308/150‑1 (MQ=255)
            cGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTAc        >  2:1373240/1‑151 (MQ=19)
            cGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTAc        <  1:5618608/151‑1 (MQ=16)
            cGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGGTTTTTGGTTAc        >  2:7178050/1‑151 (MQ=255)
             gAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACt       >  1:1835126/1‑151 (MQ=16)
             gAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACt       >  2:5502339/1‑151 (MQ=16)
               agtgAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACttt     <  2:2203612/148‑1 (MQ=2)
               aGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACttt     >  2:3613542/1‑151 (MQ=9)
               aGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACttt     >  1:3439788/1‑151 (MQ=9)
               aGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACttt     >  2:5109184/1‑151 (MQ=9)
                gggAGGTGACGGGCGTGGATGCCCGTCGAGCGCTGAGGCGCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACtttt    <  2:6669936/151‑1 (MQ=11)
                gggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACtttt    <  2:738338/151‑1 (MQ=9)
                gggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACtttt    >  1:3674499/1‑151 (MQ=9)
                gggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACtttt    <  1:2561382/151‑1 (MQ=9)
                gggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACtttt    <  1:1415841/151‑1 (MQ=9)
                 ggAGGTGAGGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACttttt   <  1:1392473/151‑1 (MQ=9)
                 ggAGGTGACGGGCATGGTTGCCCGTCAAGCGCTGAGGCCCGATGGATGGGGCCTGCAGCGCGTACTCGCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACttttt   <  2:5419133/151‑1 (MQ=11)
                 ggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACttttt   >  1:6672267/1‑151 (MQ=9)
                  gAGGTGAGGGGCATGGATGCCCGGCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  <  2:1010968/151‑1 (MQ=255)
                  gAGGTGACGGGCATGGATGCCCGTCAAGGGCTGAGGTCCCATGGATGGGGCCGGCAGCGCGTGCTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  <  2:4535086/151‑1 (MQ=11)
                  gAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  >  2:3944631/1‑151 (MQ=2)
                  gAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  >  2:5369700/1‑151 (MQ=2)
                  gAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  >  2:3282716/1‑151 (MQ=2)
                  gAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  >  1:3477739/1‑151 (MQ=2)
                  gAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  >  1:132919/1‑151 (MQ=2)
                  gAGGTGACGGGCATGGAGGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTg  <  2:3706968/151‑1 (MQ=2)
                                                                                                                                                  |                       
GCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGACCGTAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTG  >  CP009974/4631543‑4631712

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: