Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,170,111 G→C pseudogene (57/475 nt) gatR ← DNA‑binding transcriptional repressor GatR, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,170,1110GC100.0% 237.9 / NA 61pseudogene (57/475 nt)gatRDNA‑binding transcriptional repressor GatR, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (27/34);  total (27/34)

GCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  >  NC_000913/2169977‑2170245
                                                                                                                                      |                                                                                                                                      
gCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTc                                                                                                                                   >  2:665721/1‑140 (MQ=255)
  tCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCg                                                                                                                                      >  1:665698/1‑135 (MQ=255)
  tCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCg                                                                                                                                      >  2:120328/1‑135 (MQ=255)
  tCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCAt                                                                                                                                 >  2:545552/1‑140 (MQ=255)
  tCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCAt                                                                                                                                 >  1:16245/1‑140 (MQ=255)
  tCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCCAGCCACCACAGGCCACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTTTTGGTAATCACCGGGATAATTTTGCCGGTCAt                                                                                                                                 >  2:741969/1‑140 (MQ=255)
     gCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTaa                                                                                                                              <  2:431248/140‑1 (MQ=255)
           gAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTc                                                                                                                        >  2:686885/1‑140 (MQ=255)
            aaCCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCg                                                                                                                       <  1:784688/140‑1 (MQ=255)
            aaCCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCATGATATTTTTGCCGGTCATTAATCCTTCg                                                                                                                       >  2:655646/1‑140 (MQ=255)
                 tGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGAt                                                                                                                  >  2:146434/1‑140 (MQ=255)
                 tGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGAt                                                                                                                  >  2:73568/1‑140 (MQ=255)
                 tGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGAt                                                                                                                  >  2:738983/1‑140 (MQ=255)
                 tGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGAt                                                                                                                  <  1:94723/140‑1 (MQ=255)
                        aGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCAtt                                                                                                           >  2:510533/1‑140 (MQ=255)
                          cgcgTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGt                                                                                                         <  1:665721/140‑1 (MQ=255)
                                    ggCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGTAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCg                                                                                               <  1:240527/140‑1 (MQ=255)
                                    ggCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCg                                                                                               <  1:1085319/140‑1 (MQ=255)
                                       ggACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTg                                                                                            <  1:1043256/140‑1 (MQ=255)
                                         aCCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCTTTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTc                                                                                          >  2:350998/1‑140 (MQ=255)
                                                caccacAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATa                                                                                   <  1:686885/140‑1 (MQ=255)
                                                caccacAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTCGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATa                                                                                   <  2:918857/140‑1 (MQ=255)
                                                 accacAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATaa                                                                                  <  2:16245/140‑1 (MQ=255)
                                                     caGACCGCCAGAGTAATGTATTTATTTTCGCAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGCTCTTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGt                                                                              <  2:552863/140‑1 (MQ=255)
                                                       gacgacCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCt                                                                            >  1:352402/1‑140 (MQ=255)
                                                              agagTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCTTTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTa                                                                     <  1:350998/140‑1 (MQ=255)
                                                              agagTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTa                                                                     <  1:510533/140‑1 (MQ=255)
                                                                agTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACgg                                                                               <  1:981124/128‑1 (MQ=255)
                                                                agTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACgg                                                                               >  2:981124/1‑128 (MQ=255)
                                                                agTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTAcc                                                                   >  2:325439/1‑140 (MQ=255)
                                                                agTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCATGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTAcc                                                                   <  1:655646/140‑1 (MQ=255)
                                                                   aaTGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAg                                                                <  1:545552/140‑1 (MQ=255)
                                                                   aaTGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAg                                                                <  1:738983/140‑1 (MQ=255)
                                                                    aTGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGc                                                               >  2:852266/1‑140 (MQ=255)
                                                                      gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCaa                                                             >  2:1015821/1‑140 (MQ=255)
                                                                          ttATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGt                                                         >  2:186158/1‑140 (MQ=255)
                                                                              tttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGac                                                     <  1:795705/140‑1 (MQ=255)
                                                                               ttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGaca                                                    <  1:1009181/140‑1 (MQ=255)
                                                                                  ggAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGac                                                     >  2:228093/1‑136 (MQ=255)
                                                                                  ggAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGac                                                     <  1:228093/136‑1 (MQ=255)
                                                                                      aGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGACt                                             <  2:588998/140‑1 (MQ=255)
                                                                                       gggCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTa                                            <  1:673632/140‑1 (MQ=255)
                                                                                        ggCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTaa                                           >  2:564712/1‑140 (MQ=255)
                                                                                        ggCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTaa                                           <  1:1015821/140‑1 (MQ=255)
                                                                                                 tGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTTATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTaaa                                  <  1:146434/140‑1 (MQ=255)
                                                                                                    cGGGAGACTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAActc                               <  2:352402/140‑1 (MQ=255)
                                                                                                         gaCTGTTGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagt                          <  1:564712/140‑1 (MQ=255)
                                                                                                              ttGGTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGGAAACTCTCAGTCagag                     >  1:502500/1‑140 (MQ=255)
                                                                                                                ggTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTaa                                   <  2:1103306/124‑1 (MQ=255)
                                                                                                                ggTAATCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTaa                                   >  1:1103306/1‑124 (MQ=255)
                                                                                                                    aTCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACt               <  1:852266/140‑1 (MQ=255)
                                                                                                                    aTCACCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACt               >  2:538839/1‑140 (MQ=255)
                                                                                                                       aCCGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac            <  1:538839/140‑1 (MQ=255)
                                                                                                                        ccGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcaca           >  1:621401/1‑140 (MQ=255)
                                                                                                                        ccGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcaca           >  1:112051/1‑140 (MQ=255)
                                                                                                                        ccGTGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCCGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcaca           >  2:1071764/1‑140 (MQ=255)
                                                                                                                           tGATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGa        >  2:847259/1‑140 (MQ=255)
                                                                                                                            gATATTTTTGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc       <  2:333653/140‑1 (MQ=255)
                                                                                                                                 ttttGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa  <  2:621401/140‑1 (MQ=255)
                                                                                                                                 ttttGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa  <  1:325439/140‑1 (MQ=255)
                                                                                                                                 ttttGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa  <  1:1049921/140‑1 (MQ=255)
                                                                                                                                 ttttGCCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCCGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa  <  1:1071764/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAA  >  NC_000913/2169977‑2170245

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