Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,410 (T)7→6 100% coding (1280/1332 nt) gtrS → serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4040T.91.7% 265.6 / 17.4 72coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base . (34/32);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 5.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGC  >  NC_000913/2470306‑2470497
                                                                                                  |                                                                                             
aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtt                                                                                              <  1:1981011/100‑1 (MQ=255)
 ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttt                                                                                             >  1:234961/1‑100 (MQ=255)
  tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                            >  1:2251195/1‑100 (MQ=255)
   ttATTTGTCAGTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTAGTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                           <  1:1261871/100‑1 (MQ=255)
   ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                           <  2:1397980/100‑1 (MQ=255)
    tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                          <  2:2265944/100‑1 (MQ=255)
    tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                          >  2:1719569/1‑100 (MQ=255)
    tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                          <  1:668627/100‑1 (MQ=255)
      tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                        <  2:2022909/100‑3 (MQ=255)
       gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                       >  1:844539/1‑97 (MQ=255)
       gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                       >  1:827388/1‑97 (MQ=255)
        tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                      <  1:1856293/100‑5 (MQ=255)
         gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                    >  2:956068/1‑100 (MQ=255)
           cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                  >  1:2021516/1‑100 (MQ=255)
           cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                  <  2:1317691/100‑1 (MQ=255)
            aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTtat                                                                                  >  1:372636/1‑100 (MQ=255)
             tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                <  2:2296175/100‑1 (MQ=255)
              ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTTCTTATACATAACAGTCCTTTTTTAATTTATaaa                                                                               <  1:2478351/100‑1 (MQ=255)
              ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                               >  2:158490/1‑100 (MQ=255)
               ttgcGGGGAAGGATATAGAGATTTTAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaaa                                                                              <  2:353722/97‑1 (MQ=255)
                   ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                          >  2:1728401/1‑100 (MQ=255)
                     ggAAGGATATAGATATTTAAGCTATCATTTATTTGGTATGTATTTTTCTTTCAACCTTTACTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                        <  1:1726699/100‑1 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                        <  2:2321925/100‑1 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                        >  2:2129226/1‑100 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                        <  2:1623170/100‑1 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                        >  2:1154667/1‑100 (MQ=255)
                       aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaa                                                                      <  2:1574322/100‑1 (MQ=255)
                       aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaa                                                                      <  2:712633/100‑1 (MQ=255)
                        aGGTTATAGAGTTTTAAGCAAGTATTTATTTGGAATGTATTTCTCGTTAGATCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                     <  1:1818875/100‑1 (MQ=255)
                         ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACGTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                    >  2:1816540/1‑100 (MQ=255)
                         ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                    <  2:111066/100‑1 (MQ=255)
                         ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                    <  2:631591/100‑1 (MQ=255)
                            tattGATATTTAAGCAAGCAGTTATTTGGAATTTATTTTTCTTTAGACCTTTGCTTATATATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                 <  1:360149/100‑1 (MQ=255)
                            tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                 >  1:1592371/1‑100 (MQ=255)
                              taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                               <  1:1277896/100‑1 (MQ=255)
                              taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACATTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                               >  1:2343940/1‑100 (MQ=255)
                              nagagATTTAAGCAAGCATTTACTTGGAATGTATTATTCGTTCGACCTTTGCTTATACAAAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                               >  1:487501/2‑100 (MQ=255)
                                gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                             >  1:1839363/1‑100 (MQ=255)
                                 agaTTTAAGCAAGAATTTATTTTGAATGTATTTTTCGTTCGACCTTTGCTTATATATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                            <  2:1434488/100‑1 (MQ=255)
                                   aTTTTAGCAAGCATTTATTTGGAATGTATTTTTCTTTCAACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                          <  1:825237/100‑1 (MQ=255)
                                    tttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCATTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACa                                                         >  1:1408143/1‑100 (MQ=255)
                                         gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                    >  2:524138/1‑100 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                   >  2:2364588/1‑100 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                   >  2:1523780/1‑100 (MQ=255)
                                           aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                   <  1:1958562/99‑1 (MQ=255)
                                             gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                >  2:680017/1‑100 (MQ=255)
                                              caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                <  1:1098641/99‑1 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                          <  2:234963/100‑1 (MQ=255)
                                                     ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                        <  2:1478758/100‑1 (MQ=255)
                                                     ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                        <  1:2503223/100‑1 (MQ=255)
                                                       ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                      >  2:1463671/1‑100 (MQ=255)
                                                          aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                   <  1:1854919/100‑1 (MQ=255)
                                                            gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                 <  1:557000/100‑1 (MQ=255)
                                                            gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                 >  1:2308324/1‑100 (MQ=255)
                                                              aTTTTTCTTTTTACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                               >  1:805990/1‑100 (MQ=255)
                                                              aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                               >  2:872025/1‑100 (MQ=255)
                                                               tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                              <  2:1384849/100‑1 (MQ=255)
                                                                  ttcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                            >  1:2482035/1‑100 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                        <  2:2435351/100‑1 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                        <  2:2220145/100‑1 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                        >  2:10254/1‑100 (MQ=255)
                                                                       tcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                      >  2:212094/1‑100 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTTTTTTTTTTTATAATACAAAGAAAGCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                      >  1:1677254/1‑99 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                     >  2:2480297/1‑100 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                     >  2:133346/1‑100 (MQ=255)
                                                                         gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                    >  2:1042813/1‑100 (MQ=255)
                                                                           ccTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                  >  1:754265/1‑100 (MQ=255)
                                                                           ccTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCAGTAAGTTTAAATTGCGCGCCaa                   >  1:1749720/1‑99 (MQ=255)
                                                                              ttGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt               <  2:86957/100‑1 (MQ=255)
                                                                                gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg             >  2:566960/1‑100 (MQ=255)
                                                                                gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg             >  2:1683008/1‑100 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTTTTTTTTATAATAATTCAAAGAAATCAAGATAATAGCGTTGTCAAGGACTACCTTCAAATTGCGCGCCacttatc              >  1:2443005/1‑92 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc            >  1:147838/1‑100 (MQ=255)
                                                                                  ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg           <  2:335378/100‑1 (MQ=255)
                                                                                  ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg           <  2:2151289/100‑1 (MQ=255)
                                                                                   tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc          >  2:632021/1‑100 (MQ=255)
                                                                                     taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc        >  2:2182888/1‑100 (MQ=255)
                                                                                      aCATAACAGTCGTTTTTTAAGTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGACGCGca       >  1:1587859/1‑100 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTTTTTTATATAATAATTCAAAGAAAGCAAGACAATAGCGCTGTAAAGCACGTAGTGAAAATTGCGCTCCAATCAtcg             >  1:811654/1‑91 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATATAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac      >  1:575928/1‑100 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac      <  2:2308338/100‑1 (MQ=255)
                                                                                         tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaat     <  2:2280600/99‑3 (MQ=255)
                                                                                           aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc  <  2:278965/100‑1 (MQ=255)
                                                                                           aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc  <  2:2149080/100‑1 (MQ=255)
                                                                                                  |                                                                                             
ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGC  >  NC_000913/2470306‑2470497

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