Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,158,172 G→T A175S (GCC→TCC)  yqhG → DUF3828 domain‑containing protein YqhG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,158,1720GT100.0% 116.0 / NA 35A175S (GCC→TCC) yqhGDUF3828 domain‑containing protein YqhG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (18/17);  total (18/17)

TAGAAAGCGTTGCGATAACCCAAACTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACC  >  NC_000913/3157927‑3158399
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
tAGAAAGCGTTGCGATAACCCAAACTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAATCCTTTGTTTCCGa                                                                                                                                                                                                                                 >  2:266543/1‑250 (MQ=255)
    aaGCGATGCGATACCCCAAACTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGCAAATGGTCGTTGGGTCATTGACGATATTGTCCGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAActct                                                                                                                                                                                                                            >  2:415287/1‑251 (MQ=255)
       cGTTGCGATAACCCAAACTGTCGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGCCCACGCAGACACTGAAAATGGTGGCGGAAAATGGTAGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAATCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCCAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCtttt                                                                                                                                                                                                                         >  2:463057/1‑250 (MQ=255)
       cGTTGCGATAACCCAAACTGACGCCGATCCTGCCGATGCCGTTGAGCGCCTTCGTATCTTTAAAGATGATAAAGAAAAGACCAAGCAGACACTGAAAATGGTGGCGGAATATGGTCGTTGGGTCATTGAAGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAACCCTTTGTTTCCGAACTCTTTg                                                                                                                                                                                                                         >  1:448684/1‑251 (MQ=255)
                  cccAAACTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCGTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGTAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCt                                                                                                                                                                                                              <  1:445125/251‑1 (MQ=255)
                    cAAACTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGGTAATAGCGAGAATGAAAAAGCGCTGGGTGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGa                                                                                                                                                                                                            <  1:454420/251‑1 (MQ=255)
                      aaCTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGATCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGAtt                                                                                                                                                                                                          <  2:106254/251‑1 (MQ=255)
                      aaCTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGTATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTGTCCGAACTCTTTGAACATATTGCTGAtt                                                                                                                                                                                                          >  2:352081/1‑251 (MQ=255)
                              ccGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTgg                                                                                                                                                                                                  <  1:58142/251‑1 (MQ=255)
                                          atgccgTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACgtgg                                                                                                                                                                                      >  2:156511/1‑251 (MQ=255)
                                                  gtgCGCTTTCGTATCTTGAAAGATGAAAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGttt                                                                                                                                                                               <  2:85203/250‑1 (MQ=255)
                                                                 tttAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGGGGTTTCCGACTCTTACCGcc                                                                                                                                                                >  1:254745/1‑250 (MQ=255)
                                                                        atgatAAAGAAAAGACCACGCAGACACTGCAAATGGCGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGtt                                                                                                                                                        >  1:156650/1‑251 (MQ=255)
                                                                                 aaaagaCCCCGCAGACACTGACCATGGCGGCGGCAAATGGACGTTGGGTAATTGACGATATTGTCAGCAATCATGGCCGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAACGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAATATCTTGCTGATTATAGCTGGCCGTTGACGTGGGTGGTTTCTGACTCTTACCGCCAGGCGGTTGATGTCTTc                                                                                                                                               >  2:439421/1‑251 (MQ=255)
                                                                                  aaagaCCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCt                                                                                                                                              >  2:367259/1‑251 (MQ=255)
                                                                                     gaCCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCtat                                                                                                                                            >  1:315779/1‑250 (MQ=255)
                                                                                                cacTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATGGCTTATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTGCCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTc                                                                                                                                <  1:210366/251‑1 (MQ=255)
                                                                                                    gAAAATGGTTGCGGAAAATGTTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGATGCCTTTGGTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAg                                                                                                                             <  1:210581/250‑1 (MQ=255)
                                                                                                                  aaaaTGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATccc                                                                                                              <  2:208058/251‑1 (MQ=255)
                                                                                                                    aatgggctttGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGACCTCTTTGAACAGATTGCTGATTATAGCTGGCCTTGGACGTGGGTGGTGTCCGACTCTTACCGCCAGGCTGTTAATTCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGa                                                                                                            <  1:24406/243‑1 (MQ=255)
                                                                                                                                  cATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCATTCTATAAAACCACCTTCAAGACGGACAATAATCCCGATGAAGATATGCaaa                                                                                              <  1:98762/251‑1 (MQ=255)
                                                                                                                                            atTGTCAGCAATCATGACAACGTGTTAGTAGCCGTTAATAGCGAGAATGAAGCAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGATTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGAGGTTTCCGACTCTTACCGCCAGGAGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGTTATGCAAATAGAACGGcc                                                                                    >  2:293121/1‑250 (MQ=255)
                                                                                                                                                                 gTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGAt                                                               >  2:402883/1‑251 (MQ=255)
                                                                                                                                                                      aCAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGtt                                                          <  1:402883/251‑1 (MQ=255)
                                                                                                                                                                                aaTAGCGAGAAGGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAag                                                 <  1:60284/250‑1 (MQ=255)
                                                                                                                                                                                   aGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTc                                             <  1:181557/251‑1 (MQ=255)
                                                                                                                                                                                    gCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCg                                            >  1:33046/1‑251 (MQ=255)
                                                                                                                                                                                       agaATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCt                                          >  1:347002/1‑250 (MQ=255)
                                                                                                                                                                                        gaATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATCTTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTCTAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATCTGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTAt                                        >  2:320687/1‑251 (MQ=255)
                                                                                                                                                                                            aaaaaaaCGCGGGCCGCTTTAGCTGCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTc                                    <  2:164043/250‑1 (MQ=255)
                                                                                                                                                                                                aaaCGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGCAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTTTTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCAcgc                                >  1:297354/1‑251 (MQ=255)
                                                                                                                                                                                                              ttAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTc                   <  1:266543/250‑1 (MQ=255)
                                                                                                                                                                                                                        tGCAAAAAGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTCCTGATTATAGCAAGCCGTGGACGTGGGTGGTTTACGACTCTTACCGCCAGGCGGTTAATGACTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTCACGAAAATCCGATCTGTTTTGGCGACGAGTCGCTATTTTCACGCGTTGATGCACTTCGAGTACTgggg        >  1:46049/1‑249 (MQ=255)
                                                                                                                                                                                                                              aaGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAAcc  <  1:458563/251‑1 (MQ=255)
                                                                                                                                                                                                                              aaGAACAGCCGGAAGCCTTTGTTTCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAAcc  <  1:367259/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
TAGAAAGCGTTGCGATAACCCAAACTGACGCCGATCATGCCGATGCCGTTGTGCGCTTTCGTATCTTTAAAGATGATAAAGAAAAGACCACGCAGACACTGAAAATGGTGGCGGAAAATGGTCGTTGGGTCATTGACGATATTGTCAGCAATCATGGCAGCGTCTTACAAGCAGTTAATAGCGAGAATGAAAAAACGCTGGCCGCTTTAGCTTCGTTGCAAAAAGAACAGCCGGAAGCCTTTGTTGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACC  >  NC_000913/3157927‑3158399

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: