Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 4,889,556 G→A intergenic (+2/‑151) hyi → / → glxR hydroxypyruvate isomerase/tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,889,5560GA100.0% 97.7 / NA 29intergenic (+2/‑151)hyi/glxRhydroxypyruvate isomerase/tartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/17);  total (12/17)

ATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTG  >  AE015451/4889261‑4889839
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                           
aTCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGCTCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGCCCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAaacaa                                                                                                                                                                                                                                                                                        >  1:656475/1‑301 (MQ=255)
                gCAGCGACAACCTGTTCCTGCAATCAGCCATCTCCCACATGCAGATCATGGAAGGTGACCTGGCTCGCACAATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTggg                                                                                                                                                                                                                                                                        <  1:642981/301‑1 (MQ=255)
                 cAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTgggg                                                                                                                                                                                                                                                                       >  1:48389/1‑300 (MQ=255)
                    ccaCAACCCGTTCCTGCAATACCACATCTCCCCCATGCAGATCCTGGAAGGTTACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCg                                                                                                                                                                                                                                                                    <  2:20531/299‑1 (MQ=255)
                                       tACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCg                                                                                                                                                                                                                                                 >  1:242750/1‑301 (MQ=255)
                                                   cacaTGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTCGCTACCAGGGCTGGGTTGGCGCGGAGTACAAGCCGCTGACACCCACCGACGCGGGCCTGGGCTGGCTGCACACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGCGAGCGGGCTCGCCCGCTAACACCGGCaaacacggtc                                                                                                                                                                                                                                     >  1:87042/1‑294 (MQ=255)
                                                          aGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCAGCGAAAACCGGCAAAGCCGGTGCCATccc                                                                                                                                                                                                                              >  1:672349/1‑300 (MQ=255)
                                                                     ggtgcacctgctcgaaactgggaggCCAACCTGTAGCTGACAAACCCAATCCGGCGGGGCGACAACCCAGTCCGCCATTAACCAGCACCCGGCGGGATCCAATACCGCTTCCGGTTCGAACACCTGGCCCGCATTGGCTACCAGGGCTGGGTGGGCGCGAAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCt                                                                                                                                                                                                                   <  2:598865/277‑1 (MQ=255)
                                                                                                           gATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAgc                                                                                                                                                                             >  1:574369/1‑301 (MQ=255)
                                                                                                            atcgaCCACATCCAGCTGGCCGACAAGCCCGGCCGCAATGAACCGGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCGGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAgct                                                                                                                                                                            <  2:66853/297‑1 (MQ=255)
                                                                                                                    catccgacgcgCCGACAACCCATGCCGCAAGGAACCAGGCACCGGCGAGCTCAACTACCGCTTTCCGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGGGGGCGCGGAGTACAAGCCGCTTACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTgag                                                                                                                                                                    <  1:578016/291‑1 (MQ=255)
                                                                                                                    aaTCCAGCGGGCCGACAACCCAGGCCGAAATGAACCAGGAACCGGCGAGATCAACTACAGCTTCATGTGAGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTgag                                                                                                                                                                    <  2:811286/300‑1 (MQ=255)
                                                                                                                               ccGAAAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAAtt                                                                                                                                                         <  2:659306/301‑1 (MQ=255)
                                                                                                                                    ccgatctgCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACaaa                                                                                                                                                    <  1:62683/294‑1 (MQ=255)
                                                                                                                                     aCCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACAGACGCTTCCTGTTCCAAAGACTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAAACTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAag                                                                                                                                                   <  2:656245/301‑1 (MQ=255)
                                                                                                                                         cggCCGCACTGAGCCAGCAACCGGCGCGATCAACAACCGCTTCCTGTTCGAACCCCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGt                                                                                                                                               <  2:372905/300‑1 (MQ=255)
                                                                                                                                           gCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGGTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATAGCACAATTACaaaagatcgg                                                                                                                                             >  2:62671/1‑294 (MQ=255)
                                                                                                                                           agttcagacgtgtgctcttccgatctgaTCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTa                                                                                                                                              <  1:707963/275‑1 (MQ=255)
                                                                                                                                              gCAAAGAACCAGGCCACGGCGAGATCAACTACCGCTTCCTGGTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACACGCCGCTGACAACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                                                                                                          <  1:489587/301‑1 (MQ=255)
                                                                                                                                                                     gaTCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAAAGGTaagatcggaaggacgtcgggtagggaa                                                                                                                    >  2:707717/1‑276 (MQ=255)
                                                                                                                                                                               cGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCgg                                                                                                          <  2:870656/301‑1 (MQ=255)
                                                                                                                                                                                  ttCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAACGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGTCACCGGcat                                                                                                       >  1:378902/1‑301 (MQ=255)
                                                                                                                                                                                        ttCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTAGCTGCTTCGCGGGCGAACCCGCTCCCACAGGAGGCCAACGCTTGCTTGAGATATCACAATTAACAAGAGCTAATTTACAATTGCTAAAACCGGTTTCCTCgcacacggcacaatggg                                                                                                 >  2:338531/1‑285 (MQ=255)
                                                                                                                                                                                                  tGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGACATCATGGGCAAGCCCATg                                                                                       >  1:646559/1‑301 (MQ=255)
                                                                                                                                                                                                                              gggCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTcaca                                                           <  1:296636/301‑1 (MQ=255)
                                                                                                                                                                                                                               ggCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAACTCGGTTTCATCGGCACCGGCATCATGGTCCAACCCATGGCTCAGAACCTGCAAAAGGCAGGTcaaca                                                          >  2:296537/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                gCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGcc                         <  1:849347/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ggAAAACCCACACCGCAAGCTGAGAACCAACTCGGGACCCGGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCt         <  2:242663/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ccccaACGCAACCTGGGAAACAACTCGGGACCCCGGGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTg  <  2:48378/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                           
ATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGTTCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTG  >  AE015451/4889261‑4889839

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: