Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,919,557 A→G 7.1% L397P (CTG→CCG)  gudX ← glucarate dehydratase‑related protein, substrate unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,919,5570AG 7.1% 147.0 / 3.5 56L397P (CTG→CCG) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base A (23/29);  new base G (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

TGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCG  >  NC_000913/2919410‑2919704
                                                                                                                                                   |                                                                                                                                                   
tGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAg                                                                                                                                                    >  1:1033960/1‑149 (MQ=255)
  aCGGCCGAAAACGGGCCGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTGATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCCCGCCCAGAc                                                                                                                                                  >  2:1013351/1‑149 (MQ=255)
         aaaaCGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGgcgc                                                                                                                                           <  1:394532/149‑1 (MQ=255)
           aaCGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                                                                                         <  1:781667/149‑1 (MQ=255)
           aaCGGGACGTTTACGGGCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCGGTCTGGgcgtgg                                                                                                                                         <  1:638338/149‑4 (MQ=255)
              gggACGTTTACGGGCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAg                                                                                                                                      <  2:579885/149‑1 (MQ=255)
                           tCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCACCCCCCAGACCGGGCGCATCAGGAACGGCAAtatt                                                                                                                         >  2:1058382/1‑149 (MQ=255)
                             aaaGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCAAGACCGGGCGCATCAGGAAGGGCATTTTTTc                                                                                                                       >  2:293095/1‑149 (MQ=255)
                               aGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCg                                                                                                                     >  2:387952/1‑149 (MQ=255)
                                  tCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCAGACCGGGCGCATCAGGAACGGAAAAATTGCCGttcc                                                                                                                 <  2:713800/149‑3 (MQ=255)
                                       ccGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAAt                                                                                                             <  2:136946/149‑1 (MQ=255)
                                          gggATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAAtctc                                                                                                          <  2:5359/149‑1 (MQ=255)
                                             aTCAGGTACTGCCTCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAg                                                                                                       <  2:273865/149‑1 (MQ=255)
                                               cagcgaCTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGATCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCg                                                                                                     <  2:654145/144‑1 (MQ=255)
                                                 ggTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGATGTTTATAGGCCTCCTGTGCCTTTTGTACCTGTTCCCAGACCAGTTCCACGCCCAGCCCGGGCGCATCAGGAACGGCAAGTTTTCCGATGATAATCTCCAGCGGa                                                                                                   >  1:1141609/1‑149 (MQ=255)
                                                      tGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                                                                                              <  2:1118163/149‑1 (MQ=255)
                                                          tCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTc                                                                                          >  1:83431/1‑149 (MQ=255)
                                                          tCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTc                                                                                          <  1:370514/149‑1 (MQ=255)
                                                           cGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                                                                                         <  1:269920/149‑1 (MQ=255)
                                                             gACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAgg                                                                                       <  1:268863/149‑1 (MQ=255)
                                                              aCCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGc                                                                                      >  1:359476/1‑149 (MQ=255)
                                                              aCCTGCGTCGTTACGCGCACCGCCAGGCAGACGTGTATAGGCCTGATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCACAGACCGGGCGCAGCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGc                                                                                      <  2:222482/149‑1 (MQ=255)
                                                                 tGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGAc                                                                                   <  1:986109/149‑1 (MQ=255)
                                                                  gCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACa                                                                                  >  1:582563/1‑149 (MQ=255)
                                                                        ttACGCGCCCCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTGTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATGAGGAACCGCAATTTTTCCGTTTTTAATCTCCAGCGGCGTTTGGGTCAGGCGTCACTCGcc                                                                            >  2:878460/1‑149 (MQ=255)
                                                                        ttACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCCCGCCCAGCCCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGcc                                                                            >  1:9522/1‑149 (MQ=255)
                                                                        ttACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGcc                                                                            <  2:1034440/149‑1 (MQ=255)
                                                                        ttACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGcc                                                                            <  2:899282/149‑1 (MQ=255)
                                                                                  cGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAGCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCaa                                                                  >  1:1084996/1‑149 (MQ=255)
                                                                                    ccAGGCAGACGTTTATAGGCCTCTTGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCGTGCCAAAt                                                                <  1:198475/149‑1 (MQ=255)
                                                                                     cAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATc                                                               <  1:671768/149‑1 (MQ=255)
                                                                                          aGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTg                                                          <  1:491555/149‑1 (MQ=255)
                                                                                                  ataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGa                                                  <  1:1147831/149‑1 (MQ=255)
                                                                                                    aGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATa                                                >  2:162314/1‑149 (MQ=255)
                                                                                                             gtgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGATCGGCAATTTTTCCGTTTTTAATGTCGAGCGGATTTTGGGTCAGGCTACAATCGCCCTGCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGa                                       <  2:395725/149‑1 (MQ=255)
                                                                                                                  ttttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTAcc                                  <  1:999427/149‑1 (MQ=255)
                                                                                                                   tttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCCCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTAccc                                 >  1:911507/1‑149 (MQ=255)
                                                                                                                            gTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGATCGGCAATTTTTCAGTTTTTAATCTCCCCCGCATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCCGTGGGTATCGATAGCGCTAGGATTACCCGGTgccgcc                        >  1:745035/1‑149 (MQ=255)
                                                                                                                              tCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACATTCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCcgccgc                      <  1:1058382/149‑1 (MQ=255)
                                                                                                                               cccAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCcgcg                     >  2:917990/1‑149 (MQ=255)
                                                                                                                               cccAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCcgcg                     >  2:556693/1‑149 (MQ=255)
                                                                                                                                ccagtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCgcgc                    >  2:18949/1‑149 (MQ=255)
                                                                                                                                  agtccagtTCCACGCCCAGACCGGGCGAATCAGGAACGGCAATTTTTCCGTTGTTAATATCAAGCGGATTTTGGGTCAGGCGACAATCGCCCGCCAGCCAAATCCAGGGGGTAGCGATAGCGGTTGGATTTCCCGGTGCCGCCGCGccc                  >  2:76769/1‑149 (MQ=255)
                                                                                                                                     ccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCcaca               >  2:179394/1‑149 (MQ=255)
                                                                                                                                      cagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACAt              <  2:185787/149‑1 (MQ=255)
                                                                                                                                      cagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACAt              <  1:387952/149‑1 (MQ=255)
                                                                                                                                      cagagccACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGGCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGGATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACAt              <  1:1013351/144‑1 (MQ=255)
                                                                                                                                       agtTCCACGCCCAGACAGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATg             >  1:561032/1‑149 (MQ=255)
                                                                                                                                        gtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg            <  2:156999/149‑1 (MQ=255)
                                                                                                                                         tgatgcGCCCGGTCTGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATggg           >  2:735108/6‑149 (MQ=255)
                                                                                                                                         tgatgcGCCCGGTCTGGGCGCAGCAGGAACGGCAATTTTTCCGTGTTTAATCTCCAGCGGATTTTGGGTCAGGCGCCAATCGCCCTCCTGCCATCTCCAGTGGGGATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCGCACATggg           >  1:697381/6‑149 (MQ=255)
                                                                                                                                           ccACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAAACGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTa         <  1:273258/149‑1 (MQ=255)
                                                                                                                                              cGCCCGGTCTGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAAc      <  1:735108/149‑1 (MQ=255)
                                                                                                                                                cccAGACCGGGCGCATCAGGAAAGGCAATTTTTCCGTTTTTAATCTCCAGAGGATTTTGGGTCAGTCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATCTCGGTAGGATTACCCGGTGCCGCCGCGCCCACAGGGGTATACAt    >  2:1087984/1‑149 (MQ=255)
                                                                                                                                                 ccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATc   >  1:255786/1‑149 (MQ=255)
                                                                                                                                                  cAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCg  <  2:1030923/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCG  >  NC_000913/2919410‑2919704

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: