Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,777,977 A→T intergenic (+195/+169) ypjF → / ← ypjA CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,9770AT100.0% 462.9 / NA 115intergenic (+195/+169)ypjF/ypjACP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (48/67);  total (48/67)

Only 100 of 115 total aligned reads displayed.

AGGCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGG  >  NC_000913/2777744‑2778206
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
aGGCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCcca                                                                                                                                                                                                                                  <  1:1386755/239‑1 (MQ=255)
               ttGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   <  1:86554/239‑1 (MQ=255)
               ttGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   <  1:1126205/239‑1 (MQ=255)
               ttGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   <  1:1368620/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAATCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACCCAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:1319547/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCCCCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:719200/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:496579/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  1:9359/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:587107/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  1:587931/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:918623/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:95584/239‑1 (MQ=255)
                tGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCGCCCCTCCCCCGAAAAAACCTGACATTTGCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:461726/239‑1 (MQ=255)
                  cTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTCCACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTgg                                                                                                                                                                                                                    >  1:498355/1‑235 (MQ=255)
                  cTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACgg                                                                                                                                                                                                                >  2:157862/1‑239 (MQ=255)
                            aTGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGaca                                                                                                                                                                                                      >  1:232061/1‑239 (MQ=255)
                                gAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGacaaca                                                                                                                                                                                                   <  1:60223/238‑1 (MQ=255)
                                  aaTAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACAcc                                                                                                                                                                                                 <  1:157862/238‑1 (MQ=255)
                                                  cAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACAACAAAAACAGGAAGTTa                                                                                                                                                                               >  1:507050/1‑239 (MQ=255)
                                                     agCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGa                                                                                                                                                                             >  1:505682/1‑239 (MQ=255)
                                                     agCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGa                                                                                                                                                                             >  1:201752/1‑239 (MQ=255)
                                                     agCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGa                                                                                                                                                                             >  1:1281603/1‑239 (MQ=255)
                                                      gCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGa                                                                                                                                                                                    <  1:272134/231‑1 (MQ=255)
                                                       cAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGac                                                                                                                                                                                                       >  1:1156966/1‑211 (MQ=255)
                                                       cAGACTGAAGGAAAGCAAAATGCTAAGCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGac                                                                                                                                                                                                       <  2:1156966/211‑1 (MQ=255)
                                                            tGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCa                                                                                                                                                                      >  2:108518/1‑239 (MQ=255)
                                                               aGGAAAGCAAAATGCTAATCTCATAGACGAAGAGTCTCCCGCTGTAACCTCCCCGCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGtct                                                                                                                                                                        <  2:888608/234‑1 (MQ=255)
                                                               aGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:560322/192‑1 (MQ=255)
                                                               aGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGt                                                                                                                                                                          <  1:1427541/232‑1 (MQ=255)
                                                                    aGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  1:1184888/187‑1 (MQ=255)
                                                                             cTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAAt                                                                                                                                                                                                                             >  2:560511/1‑167 (MQ=255)
                                                                             cTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAAt                                                                                                                                                                                                                             <  1:560511/167‑1 (MQ=255)
                                                                               aaTCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACCCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACg                                                                                                                                                   <  1:216760/239‑1 (MQ=255)
                                                                                     aTAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGACTCAGCAAAACACCGACCAGACGGTgagg                                                                                                                                             >  1:1119680/1‑239 (MQ=255)
                                                                                      tAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGCgagga                                                                                                                                            >  2:158928/1‑239 (MQ=255)
                                                                                       aGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgaggag                                                                                                                                           <  2:959671/239‑1 (MQ=255)
                                                                                            aaGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  >  2:855149/1‑163 (MQ=255)
                                                                                            aaGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  1:855149/163‑1 (MQ=255)
                                                                                             agagACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATaa                                                                                                                                     <  1:802892/239‑1 (MQ=255)
                                                                                                         tGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTgg                                                                                                                                                                                                                    <  1:699290/148‑1 (MQ=255)
                                                                                                                 cccctccccCGAAAAAACCTTACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCcg                                                                                                                 <  1:357764/239‑1 (MQ=255)
                                                                                                                  ccctccccCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCcgc                                                                                                                <  2:959767/239‑1 (MQ=255)
                                                                                                                     tccccCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCgg                                                                                                              >  1:846573/1‑238 (MQ=255)
                                                                                                                            aaaaaaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGGGAGGAGACATAAAAGGAAAAGCAAAGAAGCCAGGGCTCCTGGt                                                                                                      >  2:548091/1‑239 (MQ=255)
                                                                                                                                 aCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   >  1:343036/1‑125 (MQ=255)
                                                                                                                                 aCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   <  2:343036/125‑1 (MQ=255)
                                                                                                                                 aCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   >  2:24997/1‑125 (MQ=255)
                                                                                                                                 aCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACaaa                                                                                                                                                                                           <  2:82032/149‑1 (MQ=255)
                                                                                                                                 aCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACaaa                                                                                                                                                                                           >  1:82032/1‑149 (MQ=255)
                                                                                                                                     gACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGaaa                                                                                             >  1:160792/1‑239 (MQ=255)
                                                                                                                                          tttCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCca                                                                                        <  1:417522/239‑1 (MQ=255)
                                                                                                                                                   ggACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCa                                                                                                                                                                   <  1:976592/155‑1 (MQ=255)
                                                                                                                                                        aacaaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  2:351849/103‑1 (MQ=255)
                                                                                                                                                        aacaaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  >  1:351849/1‑103 (MQ=255)
                                                                                                                                                        aacaaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCttt                                                                          >  2:14136/1‑239 (MQ=255)
                                                                                                                                                         acaaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCca                                                                                        >  1:735671/1‑224 (MQ=255)
                                                                                                                                                         acaaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCca                                                                                        <  2:735671/224‑1 (MQ=255)
                                                                                                                                                          caaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTgg                                                                                                                                                                                                                    >  2:450435/1‑99 (MQ=255)
                                                                                                                                                          caaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTgg                                                                                                                                                                                                                    <  1:450435/99‑1 (MQ=255)
                                                                                                                                                                gACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGa                                                                  <  2:507050/238‑1 (MQ=255)
                                                                                                                                                                    aaaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgagga                                                                                                                                            <  1:543029/161‑1 (MQ=255)
                                                                                                                                                                     aaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  1:813931/90‑1 (MQ=255)
                                                                                                                                                                     aaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  >  2:813931/1‑90 (MQ=255)
                                                                                                                                                                     aaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGa                                                                                               <  2:49646/205‑1 (MQ=255)
                                                                                                                                                                     aaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGa                                                                                               >  1:49646/1‑205 (MQ=255)
                                                                                                                                                                             aTTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGAc                                                     <  2:160792/239‑1 (MQ=255)
                                                                                                                                                                                                tcaatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGcc                                                                                                                                                                                                                                    <  2:145156/45‑1 (MQ=255)
                                                                                                                                                                                                tcaatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   <  2:943342/62‑1 (MQ=255)
                                                                                                                                                                                                tcaatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   >  1:943342/1‑62 (MQ=255)
                                                                                                                                                                                                tcaatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGa                                                                                                                                                                                    >  1:318032/1‑93 (MQ=255)
                                                                                                                                                                                                tcaatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGa                                                                                                                                                                                    <  2:318032/93‑1 (MQ=255)
                                                                                                                                                                                                  aatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAgg                                                                                                                                                                                     >  2:1274592/1‑90 (MQ=255)
                                                                                                                                                                                                  aatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAgg                                                                                                                                                                                     <  1:1274592/90‑1 (MQ=255)
                                                                                                                                                                                                    tttAACAACAAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGCACGACCAGTTCACACCAGCTACCGCGt                               >  1:1068516/1‑238 (MQ=255)
                                                                                                                                                                                                        ccaccaaGTTACCCACCCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCAc                          <  2:498355/235‑1 (MQ=255)
                                                                                                                                                                                                        acaacaAGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCAc                          <  1:548091/239‑1 (MQ=255)
                                                                                                                                                                                                            caaGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  >  2:490049/1‑51 (MQ=255)
                                                                                                                                                                                                            caaGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGc                                                                                                             <  2:926554/152‑1 (MQ=255)
                                                                                                                                                                                                            caaGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGc                                                                                                             >  1:926554/1‑152 (MQ=255)
                                                                                                                                                                                                            caaGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGAGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCCCCCGGTAACATGAAAGCCCACAGATGTGGGCTTATTCGTTGATGGTCAGAACGACAAGATCAAACCACCTACCGCTTTCCACgggg                      >  2:606200/1‑239 (MQ=255)
                                                                                                                                                                                                              aGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  1:1330556/49‑1 (MQ=255)
                                                                                                                                                                                                              aGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  >  2:1330556/1‑49 (MQ=255)
                                                                                                                                                                                                              aGTTACACAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAt                    >  2:592890/1‑239 (MQ=255)
                                                                                                                                                                                                                  acacAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGa                                                                                                                                                     <  2:563281/106‑1 (MQ=255)
                                                                                                                                                                                                                  acacAACCAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGa                                                                                                                                                     >  1:563281/1‑106 (MQ=255)
                                                                                                                                                                                                                        ccAATTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCg           >  2:531792/1‑238 (MQ=255)
                                                                                                                                                                                                                         cAATTCAGACTCCGCCTGCCCACCAATCATGATTGGAc                                                                                                                                                                                                                  <  1:510160/38‑1 (MQ=255)
                                                                                                                                                                                                                          aaTTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTc                  <  2:1281930/229‑1 (MQ=255)
                                                                                                                                                                                                                           aTTCAGACTCCGCCTGCCCACCAATCATGATTGGa                                                                                                                                                                                                                   <  1:166538/35‑1 (MQ=255)
                                                                                                                                                                                                                           aTTCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGCCCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGc         >  2:316154/1‑237 (MQ=255)
                                                                                                                                                                                                                            ttCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGCCCAGACGGTGAGGAGATATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACt      <  2:1119680/239‑1 (MQ=255)
                                                                                                                                                                                                                             tCATACTCCGCCTGCCCACCAATCATGATTGGACTGGGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCGCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTa     <  2:201752/239‑1 (MQ=255)
                                                                                                                                                                                                                             tCAGACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTa     <  1:14136/239‑1 (MQ=255)
                                                                                                                                                                                                                                gACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGa                                                                                                                                                     <  1:1302031/92‑1 (MQ=255)
                                                                                                                                                                                                                                gACTCCGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGa                                                                                                                                                     >  2:1302031/1‑92 (MQ=255)
                                                                                                                                                                                                                                     cGCCTGCCCACCAATCATGATTGGACGGTGTAAgg                                                                                                                                                                                                         >  2:166538/1‑35 (MQ=255)
                                                                                                                                                                                                                                     cGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGt                                                                                                                                                                          <  2:942418/66‑1 (MQ=255)
                                                                                                                                                                                                                                     cGCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGt                                                                                                                                                                          >  1:942418/1‑66 (MQ=255)
                                                                                                                                                                                                                                      gCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACt      <  1:1379740/229‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCTGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACt      >  2:1379740/1‑229 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
AGGCGCTCTACCGGCTTGCTAAAAACTAATGTGAAATAAACGCTTAAATACAGAGCAGACTGAAGGAAAGCAAAATGCTAATCTCATAGACGAAGAGACTCCCGCTGTAACCTCCCCTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGG  >  NC_000913/2777744‑2778206

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