Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,103,656 Δ1 bp coding (1320/1374 nt) cpxA ← sensory histidine kinase in two‑component regulatory system with CpxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,103,6540G.98.2% 242.9 / ‑3.3 56coding (1322/1374 nt)cpxAsensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (25/30);  minor base C (1/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGA  >  NC_000913/4103512‑4103791
                                                                                                                                              |                                                                                                                                         
ccTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTCTACAGCGGCAACCAAATCACCAGCCCTCACCCGCCCAGCg                                                                                                                                           >  2:843374/1‑143 (MQ=255)
   gATAAGACGCTATCAGCGTCGCATCAGGTATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggct                                                                                                                                        <  1:452863/143‑3 (MQ=255)
    aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCcgcgcctg                                                                                                                                       >  1:900505/1‑139 (MQ=255)
       aGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGtct                                                                                                                                   <  2:19217/143‑1 (MQ=255)
         aCGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTAAAGTCAAGTTCGAAGTTTAACCCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGtcttc                                                                                                                                 <  1:1092917/143‑1 (MQ=37)
           tcTATCTGCGTCGCTTCAGGCATCCTGCCCAAGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCAAATGCGAAGCTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTg                                                                                                                               <  2:643467/142‑1 (MQ=37)
             tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGCTGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAACTCACCCGCCGTACACCGCCCCGCGGCCGCCTTCTGcc                                                                                                                             >  1:301015/1‑143 (MQ=37)
             tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGcc                                                                                                                             >  2:543494/1‑143 (MQ=255)
             tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGcc                                                                                                                             <  2:516349/143‑1 (MQ=255)
              aTCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCt                                                                                                                            >  1:178804/1‑143 (MQ=255)
               tCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCtt                                                                                                                           >  1:1098557/1‑143 (MQ=255)
               tCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCtt                                                                                                                           <  1:18558/143‑1 (MQ=255)
                    gTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCcccc                                                                                                                      >  1:883540/1‑143 (MQ=255)
                    gTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCAccc                                                                                                                      <  2:1004822/143‑1 (MQ=255)
                       gCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGc                                                                                                                   <  1:861168/143‑1 (MQ=255)
                            aGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCCGCGGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                              >  1:79585/1‑143 (MQ=255)
                            aGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                              <  1:610847/143‑1 (MQ=255)
                            aGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                              <  1:71808/143‑1 (MQ=255)
                              gCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATCCAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTCCTGCCTTCACCCAGCCACGAtg                                                                                                            >  2:839811/1‑143 (MQ=255)
                                   cTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATgctgct                                                                                                       <  1:835848/143‑1 (MQ=255)
                                    tGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGctgctg                                                                                                      >  1:1026463/1‑143 (MQ=255)
                                     gCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGa                                                                                                     >  1:82700/1‑143 (MQ=255)
                                     gCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGa                                                                                                     <  1:845438/143‑1 (MQ=255)
                                     gCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCACTTATTCAGCGACAACCAAATCACCAGACGTAAAATCATCACCGGCTATCATCTGCCTTCACCACGCCACTATACTGCTGa                                                                                                     >  2:650688/1‑143 (MQ=21)
                                        cAAACGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATg                                                                                                  <  1:978756/143‑1 (MQ=255)
                                          aaTGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGc                                                                                                >  2:306943/1‑143 (MQ=255)
                                               cGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTACCTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGttt                                                                                           <  2:346389/143‑1 (MQ=255)
                                               cGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATGCAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGttt                                                                                           <  2:299062/143‑1 (MQ=255)
                                                  atccggcgTCACCGCCTTATCCCGCCTGCAAATGCGACGTTTACCCCCGCTTCTACAGCGGCAACCAACTCCCCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTcaa                                                                                        <  1:62112/140‑1 (MQ=255)
                                                  aTGCGGCGTAAACGCCTTAGCCTGCCTGCAAATGCGAAGTTCAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTcaa                                                                                        <  1:1002549/143‑1 (MQ=255)
                                                      ggcgTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCCGCCACGATGCTGCTGAATGGCGGTTTCAACAAt                                                                                    >  1:598113/1‑143 (MQ=255)
                                                       gcgTAAACGCCTTATCCTGCCTGCAACTGCGAAGTTTAACTCCGCTTATACAGCGGCCACCAAATCACAAGCCGTAAACCGCCCCGCGGCGGTCTGCTGCCATCACCCAGCAGCGATGCTGCTGACTGCCGGTGTCAACAATc                                                                                   >  2:441961/1‑143 (MQ=255)
                                                       gcgTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATc                                                                                   <  2:96820/143‑1 (MQ=255)
                                                       gcgTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATc                                                                                   <  1:306943/143‑1 (MQ=255)
                                                                    acgctgcctgcAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAAc                                                                      <  2:797310/140‑1 (MQ=255)
                                                                    aTCCTGCCTGCAAATTCGACGTTGAACTCCGCTTACACAGCGGCAACCAAATCGCCAGCCGTAAACCGCCCCGCGGCTGTCTTCTGCCTTCACCCCGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAAc                                                                      <  2:805573/143‑1 (MQ=255)
                                                                     tCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAGcc                                                                     >  1:641141/1‑143 (MQ=255)
                                                                               aaaTGCGAAGGTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCa                                                           <  1:774020/143‑1 (MQ=255)
                                                                                           tAACTCCGCTTATACAGCGGCAACCAAATCACCACCCGTACACCGCCCAGCGCCCCTCTTCTGCCTTCACCCACCCACGATGCTGCTGAATGGCGGTCTTACCCATCGCCACCCCCAACCCTCTCACGCCAGACCCACCATcg                                               >  2:392879/1‑143 (MQ=255)
                                                                                               tCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATcgcggt                                           >  1:368023/1‑141 (MQ=255)
                                                                                               tCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACCACcgcgca                                           >  2:562834/1‑141 (MQ=255)
                                                                                                 cGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCt                                         <  1:285574/143‑1 (MQ=255)
                                                                                                 cGCTTAGACAGCGGCGCGCGAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTGCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCt                                         <  1:843374/143‑1 (MQ=255)
                                                                                                          aGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTa                                >  1:817542/1‑143 (MQ=255)
                                                                                                          aGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCCGAATGGCGGTTTCAACAATCCCCAGCCCCACACCAGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTa                                <  1:220300/143‑1 (MQ=255)
                                                                                                             ggCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGa                             <  1:277390/143‑1 (MQ=255)
                                                                                                              gCAACCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTAGCGCCAGCTTCACGATCGCGTGCTCCATCGGCCCGAt                            >  2:508482/1‑143 (MQ=255)
                                                                                                                aaCCAAATCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCCGTCTCAACAATCGCCAGCCCCAAAGCTGGACCGCCAGACTCACGATCGCGTCCTCCCTCGGCACGCTac                          >  2:159438/1‑142 (MQ=255)
                                                                                                                    aaaTCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAacg                      <  1:1035624/143‑1 (MQ=255)
                                                                                                                     aaTCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAacgg                     >  1:985418/1‑143 (MQ=255)
                                                                                                                     aaTCACCAGCCGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATACTGCTCGATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGACAGATTCACGATCGCACGCTTCCTCGGTACGATAGAacgc                     >  2:130684/1‑142 (MQ=255)
                                                                                                                               cGTAAACCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATc           <  1:1040193/143‑1 (MQ=255)
                                                                                                                                  aaaCCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGt        <  2:631227/143‑1 (MQ=255)
                                                                                                                                    aCCGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGCGGGAGAATCGGtgc      >  1:675791/1‑141 (MQ=255)
                                                                                                                                     ccGCCCAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTcg     >  1:1005995/1‑143 (MQ=255)
                                                                                                                                      cGCCCAGCGGCTGTCCTCTGCCCCCACCCAGCCACGATGCTGCTGAATGGCTGTTTCACCGACCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTcgc    <  1:405697/143‑1 (MQ=255)
                                                                                                                                        cccAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGa  <  1:439668/143‑1 (MQ=255)
                                                                                                                                        cccAGCGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGa  >  1:644420/1‑143 (MQ=255)
                                                                                                                                        cccAGCGGCTGGCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTGCCACAATCGCCAGCCCCAATCCTGTACCGCCAGCTTAACGATCGTGATTGTCATCGGTTCGATAGCACGGGCGGAGCATCTGTGcgcgg  >  1:40687/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                         
CCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGA  >  NC_000913/4103512‑4103791

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