Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,485,658 Δ5,170 bp [yhfA][pabA] [yhfA], crp, yhfK, argD, yrfJ, [pabA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3485658 3490827 5170 36 [0] [1] 37 [yhfA]–[pabA] [yhfA],crp,yhfK,argD,yrfJ,[pabA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 34856570 (0.000)36 (0.980) 35/474 0.2 100% coding (162/405 nt) yhfA OsmC family protein YhfA
?NC_000913 3490828 = 0 (0.000)coding (2/564 nt) pabA aminodeoxychorismate synthase subunit 2

GAATCGGCAAGCGCCCGGCATAGTGGCGAATAACGTCAAGGGAGATCCCGGCTTCATCTGGCGTACAGGGGCCAGGTGAGATGACAATTTTTTGTGGTTTAAGGGCGTCGATATCCGCCAGCGTCAACGCATCGTTGCGCTTAACCAGCACATCCGCCCCCAGTTCACAAAAGTACTGGTAGAGGTTCCAGGTAAAAGAATCGTAGTTATCTATAAGCAGGATCATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGCCAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATT  >  NC_000913/3490604‑3491074
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
acgagtgagtaatattcaccgctttttccagcatcaaaaacaatgaacaatatttctcgacagagagatcaacagcacgcgcaactgctgcgtctttcaggacgcgaacggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacactcgaccacatcctgccgccctttttgcaggatcATAGCGGATCCGGTACAAAAGAACcgc                                                                                                                                                                                                                              >  1:341415‑M2/225‑251 (MQ=255)
        gtaatattcaccgctttttccagcatcaacgccactgaacaatatttctcggcagagagatcaaccgcacgcgcaaccgctgcgtctttcaggtcgcgaccggtgacgataaaatgcagattaatgtgtgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTAtttt                                                                                                                                                                                                                      <  2:295594‑M2/35‑1 (MQ=255)
           atattcaccgctttttccagcatcaacgccactgaacaatatttctcggcagagagatcaaccgcacgcgcaaccgctgcgtctttcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAAACGCCCTATTTTACt                                                                                                                                                                                                                   >  2:65815‑M2/214‑251 (MQ=255)
                                        cactgaacaatatttctcggcagagagatcaaccgcacgcgcaaccgctgcgtctttcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTg                                                                                                                                                                                      <  2:300011‑M2/67‑1 (MQ=255)
                                         actgaacaatatttctcggcagagagatcaaccgcacgcgcaaccgctgcgtctttcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCCTCGCTTACCACTTTGc                                                                                                                                                                                     >  2:344991‑M2/184‑251 (MQ=255)
                                            gaacaatatttctcggcagagagatcaaccgcacgcgcaaccgctgcgtctgtcaggtcgcgaccggtgacgagaaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcatgatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCTGCTTCGCTTACCACTTTGCTa                                                                                                                                                                                   <  2:172999‑M2/70‑1 (MQ=255)
                                               caatatttctcggcagagagatcaaccgcacgcgcaaccgctgcgtctttcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgccgcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAat                                                                                                                                                                                >  1:207176‑M2/178‑250 (MQ=255)
                                                     ttctcggcagagagatcaaccgcacgcgcaaccgctgcgtctttcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCctct                                                                                                                                                                         >  2:212529‑M2/172‑251 (MQ=255)
                                                           gcagagagatcaaccgcacgcgcaaccgctgcgtctttcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTg                                                                                                                                                                   <  2:42382‑M2/86‑1 (MQ=255)
                                                                              gcgcaaccgctgcgtctttcaggtcgcgaccggggacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccatttttgcaggatcATGGCGGCCCCGGGACAAAAGAACCGCGCTTTTTTACTCAGTTTCCCCGGCTTCGCTTACCACTTTGCTAAATATTATCTGCAGTTCGGTGCGATCCCCCATTGc                                                                                                                                                >  2:188248‑M2/147‑251 (MQ=255)
                                                                                        tgcgtatttcaggtcgcgacgggtgacgataaaatgcagattaatgtgcgtcaacaggcgtggtgcctcttcgcggcgttaagaggtcaattttacttcaccatcgaacacatcctgacgccctttttgcaggatcATGGCGGCCACGGTACAAAAGAACAGCGCTATTTTACTCAGATTCCCCGGCCTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGCGAGATCCACCATTGCCCCACTCTg                                                                                                                                       >  2:366372‑M2/137‑250 (MQ=255)
                                                                                                tcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCCCTCTGATTTGc                                                                                                                                 >  2:16891‑M2/129‑248 (MQ=255)
                                                                                                tcaggtcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGc                                                                                                                                 <  1:16891‑M2/120‑1 (MQ=255)
                                                                                                   cggcgcgaccggtgacgataaaatgcagattaatgtgcgtaaacaggcgttgtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgcgctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCTCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGTTGAGATCACCCAGTGCCCCTCGCTTATTTGCTTGAtt                                                                                                                           <  2:353205‑M2/126‑1 (MQ=255)
                                                                                                           ccggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCct                                                                                                                   <  1:82944‑M2/134‑1 (MQ=255)
                                                                                                            cggtgacgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCAGTGCCCCACTCTGATTTGCCTGATTCCAGGCctc                                                                                                                  <  2:181890‑M2/135‑1 (MQ=255)
                                                                                                                gacgataaaatgcagattaatgtgcgtaaacaggcgaggcgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGACACGGTACACAAGAACCGCGCTATTTTACTCAGATTCACCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTACAGGCCTCtttt                                                                                                              >  1:218942‑M2/113‑251 (MQ=255)
                                                                                                                  cgataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTATGATTTGCCTGATTCCAGGCCTCTTTTTc                                                                                                            >  1:290995‑M2/111‑251 (MQ=255)
                                                                                                                   gataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCg                                                                                                           <  1:200520‑M2/142‑1 (MQ=255)
                                                                                                                    ataaaatgcagattaatgtgcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCg                                                                                                           >  1:95725‑M2/109‑250 (MQ=255)
                                                                                                                         atgcagattaatgtgcggaaacaggcgtggtgcctcttcgcggcattcagaggccaattttactacacaatcgaccacatccagacgccctttttgcaggaacACGGCCGCTACGGCACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGATTAACTTACCACTTTGCTAAATATCATCTGCAGTGAGGTGAGATCCCCCATTACCCCAACCTGATTTGCCTGATTCCAGGCCTCTTTTTCCATAccc                                                                                                     >  1:353979‑M2/104‑251 (MQ=255)
                                                                                                                                       gcgtaaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTGTTCGATTCCCTGCCATCTCAGTTg                                                                                       <  2:212563‑M2/162‑1 (MQ=255)
                                                                                                                                          taaacaggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGat                                                                                     >  1:277537‑M2/87‑250 (MQ=255)
                                                                                                                                              cagtcttggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAACTCAGTTGATATCCg                                                                                <  2:177009‑M2/169‑1 (MQ=255)
                                                                                                                                                ggcgtggtgcctcttcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGACCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTGTTCGATACCCTGCCAGCTCAGTTGATATCCGGc                                                                              <  2:176356‑M2/171‑1 (MQ=255)
                                                                                                                                                              tcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCa                                                                   <  1:127657‑M2/182‑1 (MQ=255)
                                                                                                                                                              tcgcggcgttcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAg                                                                  <  2:98421‑M2/183‑1 (MQ=255)
                                                                                                                                                                 cggcggtcagaggtcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAATCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTTCCTCTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATAGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATTCCCGGCCAGCTCAGTTGATATCCTGAATGAATCTCCAGTGGCt                                                             <  2:96779‑M2/188‑1 (MQ=255)
                                                                                                                                                                              tcaattttacttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCg                                                >  1:223019‑M2/51‑251 (MQ=255)
                                                                                                                                                                                        ttcacaatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAAATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGcc                                         >  1:234304‑M2/41‑248 (MQ=255)
                                                                                                                                                                                             aatcgaccacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGCCAGACCACt                                 >  1:226888‑M2/36‑251 (MQ=255)
                                                                                                                                                                                                 gaccacatcctgacgcccttgttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTACCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGCCAGACCACTTccc                             <  2:167059‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                  accacatcctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGCCAGACCACTTCCCa                            <  2:223019‑M2/221‑1 (MQ=255)
                                                                                                                                                                                                     ccatcctgacgccctttgtgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGAGTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTTCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTGCGATACCCTGCCAGCGCAGTTGATATCCGGCATGACTcgccagttgttcgaagaagatccgagttgccatttcacttcccaccc                         <  2:212822‑M2/224‑24 (MQ=255)
                                                                                                                                                                                                         cctgacgccctttttgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTCATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGCCAGACCACTTCCCACCCGGa                      >  1:170220‑M2/24‑250 (MQ=255)
                                                                                                                                                                                                                       tgcaggatcATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGCCAGACCACTTCCCACCCGGAAGGGATGCAGCACGt       >  2:245686‑M2/10‑251 (MQ=255)
                                                                                                                                                                                                                            gATCATGGCGGCTCCGGTACAACAAAACCGCGCCATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGACTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCCGTTGCTCGAAGAAGATCCGCTGTGCCAGAACACTTCCCACCCCGAAGGGATGCAGCCCGTTGAtt  >  2:354301/1‑251 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
GAATCGGCAAGCGCCCGGCATAGTGGCGAATAACGTCAAGGGAGATCCCGGCTTCATCTGGCGTACAGGGGCCAGGTGAGATGACAATTTTTTGTGGTTTAAGGGCGTCGATATCCGCCAGCGTCAACGCATCGTTGCGCTTAACCAGCACATCCGCCCCCAGTTCACAAAAGTACTGGTAGAGGTTCCAGGTAAAAGAATCGTAGTTATCTATAAGCAGGATCATGGCGGCTCCGGTACAAAAGAACCGCGCTATTTTACTCAGATTCCCCGGCTTCGCTTACCACTTTGCTAAATATCATCTGCAGTGCGGTGAGATCCCCCATTGCCCCACTCTGATTTGCCTGATTCCAGGCCTCTTTTTCGATACCCTGCCAGCTCAGTTGATATCCGGCATGAATCGCCAGTTGCTCGAAGAAGATCCGCTGTGCCAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATT  >  NC_000913/3490604‑3491074

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: