Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,824,517 Δ1 bp coding (421/576 nt) ves ← HutD family protein Ves

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,824,5140C.100.0% 236.5 / NA 57coding (424/576 nt)vesHutD family protein Ves
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (28/29);  total (28/29)

TCCGGCGAATGACCAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATCGACATCTGCCCTGCCGT  >  NC_000913/1824375‑1824646
                                                                                                                                           |                                                                                                                                    
tCCGGCGAATGACCAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGc                                                                                                                                      <  1:929894/140‑1 (MQ=255)
tCCGGCGAATGACCAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGc                                                                                                                                      >  2:929894/1‑140 (MQ=255)
    gCGAATGACCAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGcccgt                                                                                                                                  >  1:1504067/1‑138 (MQ=255)
         tGACCAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGcccgtta                                                                                                                                <  2:644042/137‑5 (MQ=255)
         tGACCAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGcccgtta                                                                                                                                >  1:644042/1‑133 (MQ=255)
            ccAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATcaca                                                                                                                         <  2:500975/140‑1 (MQ=255)
            ccAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATcaca                                                                                                                         <  2:1254893/140‑1 (MQ=255)
              aGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACaaa                                                                                                                       <  1:497513/140‑1 (MQ=255)
                  aGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACACGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAAcacc                                                                                                                   >  2:1100283/1‑140 (MQ=255)
                    ccAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAAcaccac                                                                                                                 <  2:1504067/140‑1 (MQ=255)
                      aGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCAccc                                                                                                               >  2:1086510/1‑140 (MQ=255)
                             ttCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCcgcg                                                                                                        <  1:209909/140‑1 (MQ=255)
                             ttCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCcgcg                                                                                                        <  1:816404/140‑1 (MQ=255)
                             ttCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCcgcg                                                                                                        >  1:1332938/1‑140 (MQ=255)
                                aCTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGa                                                                                                     <  1:527145/140‑1 (MQ=255)
                                 cTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCcgccgc                                                                                                         >  1:1447302/1‑135 (MQ=255)
                                 cTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCcgccgc                                                                                                         <  2:1447302/135‑1 (MQ=255)
                                  tGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGTCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGc                                                                                                   <  2:1545284/140‑1 (MQ=255)
                                   gAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGcc                                                                                                  <  2:295280/140‑1 (MQ=255)
                                   gAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGcc                                                                                                  <  1:1086510/140‑1 (MQ=255)
                                    aacaacAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCa                                                                                                 <  1:1470286/140‑1 (MQ=255)
                                    aacaacAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCa                                                                                                 >  2:557087/1‑140 (MQ=255)
                                     acaacaGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCaa                                                                                                <  1:725360/140‑1 (MQ=255)
                                        acaGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCTACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAgg                                                                                             <  1:81615/140‑1 (MQ=255)
                                         caGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGgt                                                                                            >  1:834344/1‑140 (MQ=255)
                                                ccTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAac                                                                                     <  2:1332938/140‑1 (MQ=255)
                                                  ttGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAacgt                                                                                   <  2:125946/140‑1 (MQ=255)
                                                   tGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGgta                                                                                           >  1:1447033/1‑131 (MQ=255)
                                                   tGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGCTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGgta                                                                                           <  2:1447033/131‑1 (MQ=255)
                                                                  cGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTg                                                                   >  1:527249/1‑140 (MQ=255)
                                                                       ggTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACttt                                                              <  2:834344/140‑1 (MQ=255)
                                                                                 ccATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACttt                                                              <  1:1277484/130‑1 (MQ=255)
                                                                                 ccATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACttt                                                              >  2:1277484/1‑130 (MQ=255)
                                                                                 ccATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAg                                                    >  2:341003/1‑140 (MQ=255)
                                                                                 ccATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAg                                                    >  2:313638/1‑140 (MQ=255)
                                                                                     cGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTc                                                             <  2:393122/127‑1 (MQ=255)
                                                                                     cGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTc                                                             >  1:393122/1‑127 (MQ=255)
                                                                                     cGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAg                                                    >  2:1326729/1‑136 (MQ=255)
                                                                                     cGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAg                                                    <  1:1326729/136‑1 (MQ=255)
                                                                                           aGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGc                                          >  1:1589530/1‑140 (MQ=255)
                                                                                              aCGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCgtg                                       <  1:557087/140‑1 (MQ=255)
                                                                                                 cgcCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCa                                    >  1:510615/1‑140 (MQ=255)
                                                                                                   cccTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCatg                                  <  1:341003/140‑1 (MQ=255)
                                                                                                     cTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCatgat                                >  1:195085/1‑140 (MQ=255)
                                                                                                       ggTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGt                              >  2:422074/1‑140 (MQ=255)
                                                                                                                tGAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCa                     >  1:293691/1‑140 (MQ=255)
                                                                                                                 gAGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCAt                    <  2:1589530/140‑1 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCat                                   >  2:1295302/1‑124 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCat                                   <  1:1295302/124‑1 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGaaa                         >  2:564443/1‑134 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGaaa                         <  1:564443/134‑1 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATc                   >  2:938044/1‑140 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATc                   >  1:940370/1‑140 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATc                   >  1:1277520/1‑140 (MQ=255)
                                                                                                                  aGCAATTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATc                   >  1:1156659/1‑140 (MQ=255)
                                                                                                                     aaTTTATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATCGAc                >  2:766096/1‑140 (MQ=255)
                                                                                                                       tttATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATc                   >  2:855905/1‑135 (MQ=255)
                                                                                                                       tttATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATc                   <  1:855905/135‑1 (MQ=255)
                                                                                                                       tttATCGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATCGACAt              >  2:1383387/1‑140 (MQ=255)
                                                                                                                            cGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATCGACATctgcc         <  1:766096/140‑1 (MQ=255)
                                                                                                                            cGCCCAGTTGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATCGACATctgcc         >  1:1291723/1‑140 (MQ=255)
                                                                                                                                   ttGCCACGCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATCGACATCTGCCCTGCCGt  <  1:1100283/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TCCGGCGAATGACCAGCGAGCCAGTTGATTTCACTGAACAACAGTTTACCTTGCGGTTGTAAAAGCCGCAGGGTGTGCCTGCCATCGAACCAGCACGCGCCCTGGTCGGTAGTGAGCAATTTATCGCCCAGTTGCCACGCCCCGTTAATCACAAACACCACCCCGCCGCGTGAGCCAAAGGTAGTAAACGTACGTTCGGCAATTCTGACTTTCGCCTTACAGACATCCAGGCGTGTCATGATGTTGAAATCCATCGACATCTGCCCTGCCGT  >  NC_000913/1824375‑1824646

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