Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,088,010 2 bp→TT intergenic (+150/‑273) metK → / → galP methionine adenosyltransferase/galactose:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,088,0100CT97.2% 105.0 / ‑6.7 36intergenic (+150/‑274)metK/galPmethionine adenosyltransferase/galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/18);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,088,0110CT100.0% 109.7 / NA 36intergenic (+151/‑273)metK/galPmethionine adenosyltransferase/galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/18);  total (18/18)

ACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTC  >  NC_000913/3087768‑3088252
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
aCCGCAGCATCCGGTCCTTTTGGGCTTGAACATTCCCCGGGGGAAAAACCCGAAACAGCGCAGT‑GGCTGCGCTATGCGGCCGGTCTGAAGTATGTTTTCTTCACCTGCGTTCAAAGGCCATCCACCCGCTGGCCTTCTTCTTTTGGAGAGGTGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAATCTTCCTGATCACGTTGCACCGCAGATTATCATCACAACTGAAATTGA‑TTaca                                                                                                                                                                                                                                            <  2:524846/251‑1 (MQ=255)
            ggTCACTTTGGTCGTGAACATTTCCCGGGGGAAAAACCCGACAAAGCGATGC‑TGCTGTGCGATGCTGTCGGTCTGAATTAATCTTTCTCCACCGGCGTTCAAAGGCCAGCCACGCGCTGGCCTTTTTCTTTGGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCacaaca                                                                                                                                                                                                                                <  2:416758/251‑1 (MQ=255)
                     gggcGTGAACGTTCCCCGGGGGAAAAAACCGAAAAAGC‑CTGCCTGTCGCGCGATGCTGCCGTGCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCCCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACaaatat                                                                                                                                                                                                                       <  2:262894/248‑4 (MQ=255)
                      gTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCCTCACAACTGCAATTGA‑TTACACCAACCACAACAAACAAAAAtt                                                                                                                                                                                                                      >  2:365430/1‑251 (MQ=255)
                        cGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTg                                                                                                                                                                                                                    >  1:341442/1‑251 (MQ=255)
                           aaaCATTTCCCGTGGGAAAAAACCGACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTaa                                                                                                                                                                                                                 <  1:751496/250‑1 (MQ=255)
                                   cccGTGGGAAAAAACCGACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCCCCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTc                                                                                                                                                                                                          <  1:330651/250‑1 (MQ=255)
                                    ccggggaaaaaaaCCGACATAGCGCACC‑TGTGGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCGGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCat                                                                                                                                                                                                        <  2:130663/244‑1 (MQ=255)
                                                 ccGACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGt                                                                                                                                                                                           >  1:231788/1‑251 (MQ=255)
                                                   gACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGttt                                                                                                                                                                                         >  1:516429/1‑251 (MQ=255)
                                                    aCAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGtttt                                                                                                                                                                                        >  2:412893/1‑251 (MQ=255)
                                                                 tgctgcGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTAc                                                                                                                                                                            >  2:713186/1‑251 (MQ=255)
                                                                             ctgacggtttaatgaaatattgcttcACCTGAGTACAAAAGCCAGCCACCCGCTGGCCTTTTTCTTTTGGATAGGCGTACACGACGCCTCCGGCAAAAAAAACAACCGCACAATAACATCATTCTTCATGATCACGTTTCACAGCAAATTATCATCAAAAAT‑CAACTTACTTAAACCAACCACAAAAGACAAAGAATTGTAATATTTTCATATTATTATTCGGTTTTAACAGTGGTTACAATTCTTTTCa                                                                                                                                                                 >  2:572488/23‑250 (MQ=255)
                                                                                          gTAATCTTTCTTCTCCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTTGATAGGCGCTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAAtt                                                                                                                                                   <  2:182499/251‑1 (MQ=255)
                                                                                              tctttcttCACATGCGTACAAAGGCCAGCCTCGCGCTGGCCTGTTTTTTTTGGATGGGTGTTCACGCCTCATCCGGCAAAAAAACCGCCCGCACAATAACATCATGCTTCCTGATCACGTTTCTCCTCAGATTAACATCACAACTGAAATTGA‑TTACACCATCCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACAGTTCTTTTCAGTAAAGTCTTAATTGCa                                                                                                                                                <  2:516429/250‑1 (MQ=255)
                                                                                                                   aGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCCTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTatga                                                                                                                          <  2:533731/251‑1 (MQ=255)
                                                                                                                     accAGCCTCGCGCTGGCGTTTTTCTTTTGGACCGGCGTCCACGCCGCATCCGACAAAAAAACCGCCCGCACCATAACATCCTTCTTCCTGATAACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTATAATATTTTCATATTATTATTCGGGTTTCACCGTTGTTACCTTCATTTTCAGTAAAGTCTTAATTGCAGCTAACAGCGTTTAATCttttatg                                                                                                                        >  2:379522/2‑245 (MQ=255)
                                                                                                                          ccTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGatat                                                                                                                    >  2:278171/1‑250 (MQ=255)
                                                                                                                                     ccTTTTTCTTTTGGATAGGAGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAattatt                                                                                                        <  1:311783/251‑1 (MQ=255)
                                                                                                                                            cTTTTGGATAGGCGTTCACGCCGCATCCGGCCAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGc                                                                                                 >  2:324665/1‑251 (MQ=255)
                                                                                                                                                tGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAATCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACcaca                                                                                                                                                                                                                                   <  2:239045/117‑1 (MQ=255)
                                                                                                                                                         gatCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTACTTGCAGATAACAGCGTTTAATCTATGATGATCTAACTCAATTATTTTCATGCACTTAAATCATa                                                                                     >  1:873708/3‑250 (MQ=255)
                                                                                                                                                          ttCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAAc                                                                                   >  1:739015/1‑251 (MQ=255)
                                                                                                                                                               gccgcATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTa                                                                                 >  1:566449/1‑248 (MQ=255)
                                                                                                                                                                  gcATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTGTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATaa                                                                           <  1:112206/251‑1 (MQ=255)
                                                                                                                                                                             aaaaaCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCACCCACAACAGACAAAGATTTGTTATATTTTCATAATATTATCCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTAATTGCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTa                                                                <  2:354074/251‑1 (MQ=255)
                                                                                                                                                                                aaCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGc                                                             <  2:644005/251‑1 (MQ=255)
                                                                                                                                                                                 aCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCg                                                            <  2:231788/251‑1 (MQ=255)
                                                                                                                                                                                                 catcatTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGAt                                            <  1:760601/251‑1 (MQ=255)
                                                                                                                                                                                                             tGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATc                                >  1:479501/1‑251 (MQ=255)
                                                                                                                                                                                                                 cACGTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCtttt                            >  2:66598/1‑251 (MQ=255)
                                                                                                                                                                                                                    gTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACCGACAACGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTCTTTTCATGCACTTAAATCATAACTACGATCAATATTAGTGTAAGCGATTACACTGATTTGATTTGCTTCACATCTTTTTAc                         >  2:453283/1‑251 (MQ=255)
                                                                                                                                                                                                                    gTTTCACCGCAGATTATCATCACAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACTCTGATGTGATTTGCTTCACATCTTTTTAc                         <  1:850408/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tcatcaCAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATc         >  1:701198/1‑251 (MQ=255)
                                                                                                                                                                                                                                     catcaCAACTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTAAGTCGTACTCCCctctct        >  2:776217/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aaCTGAAATTGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTc  <  1:713186/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
ACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGC‑TGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGA‑TTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTC  >  NC_000913/3087768‑3088252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑