Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 450,996 A→G G205G (GGT→GGC cyoA ← cytochrome bo3 ubiquinol oxidase subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913450,9960AG100.0% 214.9 / NA 59G205G (GGT→GGCcyoAcytochrome bo3 ubiquinol oxidase subunit 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/32);  total (27/32)

GCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATACCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTGAAGTACACCG  >  NC_000913/450862‑451132
                                                                                                                                      |                                                                                                                                        
gCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTc                                                                                                                                     >  2:177372/1‑140 (MQ=255)
      aTGTCAGACATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGt                                                                                                                               <  2:48603/140‑1 (MQ=255)
      aTGTCAGACATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGt                                                                                                                               <  2:868921/140‑1 (MQ=255)
         tCAGACATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGcc                                                                                                                            <  1:237918/140‑1 (MQ=255)
            gACATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCggg                                                                                                                         >  1:194561/1‑140 (MQ=255)
              cATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGtt                                                                                                                       <  2:903410/140‑1 (MQ=255)
              cATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGtt                                                                                                                       >  2:467810/1‑140 (MQ=255)
                 ggTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCAt                                                                                                                                   <  2:236867/125‑1 (MQ=255)
                 ggTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCAt                                                                                                                                   >  1:236867/1‑125 (MQ=255)
                  gtgtTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGtt                                                                                                                   >  2:910307/1‑140 (MQ=255)
                  gtgtTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGtt                                                                                                                   >  1:581465/1‑140 (MQ=255)
                  gtgtTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGtt                                                                                                                   >  2:68936/1‑140 (MQ=255)
                   tgtTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTg                                                                                                                  <  1:786381/140‑1 (MQ=255)
                   tgtTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTg                                                                                                                  >  2:528633/1‑140 (MQ=255)
                    gtTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGTGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTgg                                                                                                                 <  2:441906/140‑1 (MQ=255)
                       cGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGa                                                                                                              >  2:293371/1‑140 (MQ=255)
                         gCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATc                                                                                                            <  2:581465/140‑1 (MQ=255)
                           gACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAg                                                                                                          <  2:960112/140‑1 (MQ=255)
                            aCTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGa                                                                                                         >  1:278320/1‑140 (MQ=255)
                               gCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGc                                                                                                      >  2:1030103/1‑140 (MQ=255)
                                  tCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAgg                                                                                                   <  1:1030103/140‑1 (MQ=255)
                                    gCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCg                                                                                                 <  1:467810/140‑1 (MQ=255)
                                    gCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCg                                                                                                 <  1:1300/140‑1 (MQ=255)
                                    ccctttGCGACCCACTGGTCGAATGCGGCGCGGTTCGGTGTTGCAAGAGCTTTGAACTTCATGCCCGAGAAGCCCGGGCCGCTGTAGCTGGCGGTGATGCCGTCAGAAGTGCCCGGTTCGTTGGCGATCAGATGTAGGCg                                                                                                 <  2:737353/137‑1 (MQ=255)
                                       tttGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGt                                                                                              <  2:921842/140‑1 (MQ=255)
                                                  cTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCg                                                                                                 <  1:732064/126‑1 (MQ=255)
                                                  cTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCg                                                                                                 >  2:732064/1‑126 (MQ=255)
                                                    ggTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGcc                                                                                 <  2:1032287/140‑1 (MQ=255)
                                                     gTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCa                                                                                <  1:859406/140‑1 (MQ=255)
                                                           tGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCAt                                                                          >  1:1021304/1‑140 (MQ=255)
                                                               gcgcGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAAt                                                                      >  2:138678/1‑140 (MQ=255)
                                                                           gtTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAg                                                          <  1:293371/140‑1 (MQ=255)
                                                                              gCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGAcg                                                       <  1:138678/140‑1 (MQ=255)
                                                                               cAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGAcgc                                                      >  2:222090/1‑140 (MQ=255)
                                                                               cAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGAcgc                                                      >  2:552166/1‑140 (MQ=255)
                                                                               cAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGAcg                                                       >  1:508573/1‑139 (MQ=255)
                                                                               cAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGAcg                                                       <  2:508573/139‑1 (MQ=255)
                                                                                       ttGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGAcgcg                                                     <  2:766723/133‑1 (MQ=255)
                                                                                       ttGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGAcgcg                                                     >  1:766723/1‑133 (MQ=255)
                                                                                        tGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATgaag                                             >  2:340258/1‑140 (MQ=255)
                                                                                        tGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATgaag                                             >  1:461712/1‑140 (MQ=255)
                                                                                         gAACTTCATGCCTGAGAAGCCCGGGACGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATgaaga                                            >  1:1046851/1‑140 (MQ=255)
                                                                                               cATGCCTGAGAAGCCCGGGACGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGt                                      <  2:1046851/140‑1 (MQ=255)
                                                                                                          aGCCCGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAg                           <  1:528633/140‑1 (MQ=255)
                                                                                                            cccGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCAc                               >  1:571438/1‑134 (MQ=255)
                                                                                                            cccGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCAc                               <  2:571438/134‑1 (MQ=255)
                                                                                                              cGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTgg                       >  1:627774/1‑140 (MQ=255)
                                                                                                              cGGGCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTgg                       >  1:80229/1‑140 (MQ=255)
                                                                                                                ggCCGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAg                     >  2:634246/1‑140 (MQ=255)
                                                                                                                  ccGCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGt                   <  1:222090/140‑1 (MQ=255)
                                                                                                                    gCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCa                 <  2:461712/140‑1 (MQ=255)
                                                                                                                    gCTGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCa                 <  2:924670/140‑1 (MQ=255)
                                                                                                                      tGTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACt               <  2:80229/140‑1 (MQ=255)
                                                                                                                       gTAGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACtt              <  1:910307/140‑1 (MQ=255)
                                                                                                                         aGCTGGCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTg            <  2:278320/140‑1 (MQ=255)
                                                                                                                             ggCGGAGATGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTGAAGt        <  2:1021304/140‑1 (MQ=255)
                                                                                                                                  agaTGCCGTCATAANTGCCGGGTTCGTNGGNGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTGAAGTACAcc   >  1:474753/1‑140 (MQ=255)
                                                                                                                                  agaTGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTGAAGt        >  2:860678/1‑135 (MQ=255)
                                                                                                                                  agaTGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTGAAGt        <  1:860678/135‑1 (MQ=255)
                                                                                                                                   gaTGCCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTGAAGTACACCg  <  2:998000/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
GCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCGCTTTTGCGACCCACTGGTCGAATGCGGCGCGATCCGGTGTTGCAATAGCTTTGAACTTCATGCCTGAGAAGCCCGGGCCGCTGTAGCTGGCGGAGATACCGTCATAAGTGCCGGGTTCGTTGGCGATCAGATGCAGGCGAGTCTGCATACCGGCCATGGCATAAATCTGGCTACCCAGACGCGGAATGAAGAAGGAGTTCATCACGGAGTTGGAGGTCACTTTGAAGTACACCG  >  NC_000913/450862‑451132

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: