Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,926,096 G→A *219* (TAG→TAA yobB → C‑N hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,0960GA100.0% 225.1 / NA 69*219* (TAG→TAAyobBC‑N hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/30);  total (39/30)

AAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTG  >  NC_000913/1925964‑1926230
                                                                                                                                    |                                                                                                                                      
aaaTGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCtttt                                                                                                                                  <  2:61802/139‑1 (MQ=255)
          ggTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGa                                                                                                                        <  2:1034723/139‑1 (MQ=255)
              aCAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCAt                                                                                                                    >  1:111241/1‑139 (MQ=255)
                 gTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGcc                                                                                                                 >  1:659221/1‑139 (MQ=255)
                    cACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATg                                                                                                              >  2:966000/1‑139 (MQ=255)
                    cACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATg                                                                                                              <  2:111241/139‑1 (MQ=255)
                          gggATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTgcgc                                                                                                        <  1:941029/139‑1 (MQ=255)
                               gAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTAt                                                                                                   <  1:966000/139‑1 (MQ=255)
                                 atatGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCg                                                                                                 <  1:767792/139‑1 (MQ=255)
                                    tGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGcc                                                                                                                 >  1:691985/1‑120 (MQ=255)
                                    tGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGcc                                                                                                                 <  2:691985/120‑1 (MQ=255)
                                      gtcgtcTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGAAGCATGCCATGTTGCGCATTATCGATACa                                                                                            >  2:996641/1‑139 (MQ=255)
                                       tcgtcTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAg                                                                                           >  1:1126804/1‑139 (MQ=255)
                                          tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATc                                                                                                  >  2:504667/1‑129 (MQ=255)
                                          tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATc                                                                                                  <  1:504667/129‑1 (MQ=255)
                                          tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGaaa                                                                                        >  2:1321234/1‑139 (MQ=255)
                                          tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGaaa                                                                                        >  2:669251/1‑139 (MQ=255)
                                          tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGaaa                                                                                        >  2:623129/1‑139 (MQ=255)
                                          tctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGaaa                                                                                        >  1:445622/1‑139 (MQ=255)
                                           ctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTgcgc                                                                                                        >  2:939300/1‑122 (MQ=255)
                                           ctcATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTgcgc                                                                                                        <  1:939300/122‑1 (MQ=255)
                                              aTCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGa                                                                                                                        >  2:1128090/1‑103 (MQ=255)
                                              aTCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGa                                                                                                                        <  1:1128090/103‑1 (MQ=255)
                                              aTCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTgc                                                                                                          <  2:1447514/117‑1 (MQ=255)
                                              aTCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTgc                                                                                                          >  1:1447514/1‑117 (MQ=255)
                                                 gTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCgg                                                                                 <  1:669251/139‑1 (MQ=255)
                                                     gcgcCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTg                                                                             >  1:1044444/1‑139 (MQ=255)
                                                     gcgcCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATGATCGATACAGAAACCTGCGGTTTg                                                                             >  2:820153/1‑139 (MQ=255)
                                                           aTCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCg                                                                                  <  2:1471325/128‑1 (MQ=255)
                                                           aTCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCg                                                                                  >  1:1471325/1‑128 (MQ=255)
                                                             cGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggagg                                                                     <  2:445622/139‑1 (MQ=255)
                                                                       tattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGcc                                                                                                                 <  1:728401/85‑1 (MQ=255)
                                                                       tattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGcc                                                                                                                 >  2:728401/1‑85 (MQ=255)
                                                                       tattGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTgag                                                           >  2:110481/1‑139 (MQ=255)
                                                                         ttGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGAt                                                         <  1:1321234/139‑1 (MQ=255)
                                                                            ttAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGc                                                      >  1:863690/1‑139 (MQ=255)
                                                                            ttAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGc                                                      >  1:1190042/1‑139 (MQ=255)
                                                                             tAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGcc                                                     <  2:1044444/139‑1 (MQ=255)
                                                                              aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGcc                                                                                                                            >  1:608573/1‑67 (MQ=255)
                                                                              aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGcc                                                                                                                            <  2:608573/67‑1 (MQ=255)
                                                                              aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGc                                                                                                                      >  2:1018095/1‑73 (MQ=255)
                                                                              aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGc                                                                                                                      <  1:1018095/73‑1 (MQ=255)
                                                                              aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCCGGGAGGGATCGTTGAGATTGCct                                                    >  1:1291425/1‑139 (MQ=255)
                                                                              aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCct                                                    <  2:1190042/139‑1 (MQ=255)
                                                                              aGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCct                                                    >  2:18723/1‑139 (MQ=255)
                                                                                tCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGcc                                                                                                                 >  1:938640/1‑76 (MQ=255)
                                                                                tCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCct                                                    <  1:585523/137‑1 (MQ=255)
                                                                                tCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCct                                                    >  2:585523/1‑137 (MQ=255)
                                                                                acGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGcc                                                                                                                 <  2:938640/75‑1 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCtttt                                                                                                                                  >  1:344356/1‑56 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCtttt                                                                                                                                  <  2:344356/56‑1 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGcc                                                                                                                            <  2:956025/62‑1 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGcc                                                                                                                            >  1:956025/1‑62 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATc                                                                                                  <  2:188450/88‑1 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATc                                                                                                  >  1:188450/1‑88 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGtt                                               >  2:1445939/1‑139 (MQ=255)
                                                                                   gTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGtt                                               >  1:556260/1‑139 (MQ=255)
                                                                                    tCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCct                                                    <  2:416700/133‑1 (MQ=255)
                                                                                    tCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCct                                                    >  1:416700/1‑133 (MQ=255)
                                                                                                     gTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggag                                                                      <  2:274847/98‑1 (MQ=255)
                                                                                                     gTTGGCAAGGCGATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggag                                                                      >  1:274847/1‑98 (MQ=255)
                                                                                                                gATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGAt                                                         <  2:326887/100‑1 (MQ=255)
                                                                                                                gATATCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGAt                                                         >  1:326887/1‑100 (MQ=255)
                                                                                                                   atCATTCCATTACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAAGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCAt               <  1:623129/139‑1 (MQ=255)
                                                                                                                            ttACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCt      >  1:100162/1‑139 (MQ=255)
                                                                                                                            ttACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCt      >  1:1003471/1‑139 (MQ=255)
                                                                                                                             tACGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTg     >  2:1048953/1‑139 (MQ=255)
                                                                                                                               cGCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGt   <  1:1115881/139‑1 (MQ=255)
                                                                                                                                gCTAAGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTg  >  1:519076/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                      
AAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGATCGTGGTTCATTATTGTTAGTCGGTCAGCGTTCATCACAGGGTTGGCAAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTG  >  NC_000913/1925964‑1926230

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