Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,223,222 A→C 10.7% intergenic (‑529/+31) PP_2823 ← / ← PP_2824 methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,223,2220AC10.7% 158.7 / 10.2 56intergenic (‑529/+31)PP_2823/PP_2824methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (20/30);  new base C (4/2);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

AGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGC  >  NC_002947/3223074‑3223365
                                                                                                                                                    |                                                                                                                                               
aGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCa                                                                                                                                                 <  8:135553/149‑1 (MQ=255)
 gACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGGTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAg                                                                                                                                                <  3:64158/149‑1 (MQ=255)
   cGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAggg                                                                                                                                              <  5:278590/149‑1 (MQ=255)
               agagCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGt                                                                                                                                  <  5:72056/149‑1 (MQ=255)
                   ccTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTAc                                                                                                                              >  7:108929/1‑149 (MQ=255)
                      gCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                           >  8:265052/1‑149 (MQ=255)
                      gCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                           >  3:48782/1‑149 (MQ=255)
                       cAAATAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGTAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                           <  6:246168/148‑1 (MQ=255)
                       cAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                           <  5:335851/148‑1 (MQ=255)
                       cAAAAAGCAGGGCAGGCGATATCGGGGCCTGGTGCTTGGGGTGGATCGGCAAGGCGCTGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                           <  7:300650/148‑1 (MQ=255)
                        aaaaaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGGTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                           >  7:301281/1‑147 (MQ=255)
                           aaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGc                                                                                                                      >  3:162829/1‑149 (MQ=255)
                               agggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACt                                                                                                                  >  1:35445/1‑149 (MQ=255)
                                gggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCTGGAAGgc                                                                                                                                         >  5:315188/1‑125 (MQ=255)
                                gggcaggCGCTATCGGGGCCTGGGGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCTGGAAGgc                                                                                                                                         <  6:315188/125‑1 (MQ=255)
                                 ggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACt                                                                                                                  >  4:106730/1‑147 (MQ=255)
                                  gcaggCGCTATCGGGGCCGGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAg                                                                                                               <  5:87543/149‑1 (MQ=255)
                                    aggCGCTATCGGGGCCGGGGGCTTGGGTGTGCGCGGCAAGGCTCGGGGCGGCAGTTGCGCGACGACGGAAGGCGTCGAAGACACGATCGGTTGTGTGGGGGGGCCGCCTTCCAGGGAAGGCGACGGGTCTACCGGGAGGCGACTTAGGa                                                                                                             <  2:209663/149‑1 (MQ=255)
                                      gcgcTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACt                                                                                                           >  6:267714/1‑149 (MQ=255)
                                        gcTATCGGGGCGTGGTGCTGGGGGTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACTCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTg                                                                                                         <  4:195290/149‑1 (MQ=255)
                                         cTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGATGCGACTTAGGACTTgg                                                                                                        <  7:265052/149‑1 (MQ=255)
                                               gggCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGttt                                                                                                  >  4:323628/1‑149 (MQ=255)
                                                 gCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTg                                                                                                <  2:314367/149‑1 (MQ=255)
                                                  ccTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                               <  8:108929/149‑1 (MQ=255)
                                                       tGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGga                                                                                          <  2:223262/149‑1 (MQ=255)
                                                         cTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGgag                                                                                         >  7:299464/1‑148 (MQ=255)
                                                          ttGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGgaga                                                                                        >  6:213416/1‑148 (MQ=255)
                                                             ggTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATc                                                                                    <  3:323628/149‑1 (MQ=255)
                                                              gTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCa                                                                                   <  4:252774/149‑1 (MQ=255)
                                                               ttGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                  <  8:75933/149‑1 (MQ=255)
                                                                    tCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                               <  2:256723/131‑1 (MQ=255)
                                                                    tCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                               >  1:256723/1‑131 (MQ=255)
                                                                       gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTg                                                                                                         >  2:94432/1‑118 (MQ=255)
                                                                       gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTg                                                                                                         <  1:94432/118‑1 (MQ=255)
                                                                         aaGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAg                                                                        <  4:48782/149‑1 (MQ=255)
                                                                           ggCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGc                                                                      <  5:267714/149‑1 (MQ=255)
                                                                                     ggcAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAg                                                            <  2:35445/149‑1 (MQ=255)
                                                                                                      ggAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTgc                                           >  1:228261/1‑149 (MQ=255)
                                                                                                           gcggcgAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTaa                                      <  8:299464/149‑1 (MQ=255)
                                                                                                              gcgAAGACACGTGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAACGCGCCGGGTCTACCGGGATGCGATTTATGACTTGGCTGTTTTGGGTGGAGAAATCAAGACAATGAACGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACgg                                   <  2:259966/149‑1 (MQ=255)
                                                                                                                   gACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCg                              <  8:301281/149‑1 (MQ=255)
                                                                                                                     cacGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGa                            <  2:228261/149‑1 (MQ=255)
                                                                                                                           cGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGa                       >  1:27121/1‑148 (MQ=255)
                                                                                                                             gTTTGGTGGGGGGGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAg                    <  6:132352/149‑1 (MQ=255)
                                                                                                                              tttGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAgg                   >  5:331960/1‑149 (MQ=255)
                                                                                                                                tggtggGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGgc                  <  4:162829/148‑1 (MQ=255)
                                                                                                                                       ggTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat          >  6:31164/1‑149 (MQ=255)
                                                                                                                                       ggTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat          <  5:213416/149‑1 (MQ=255)
                                                                                                                                        ccggcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat          >  7:269033/5‑148 (MQ=255)
                                                                                                                                        ccggcGCCTTCCCTGGAAGACGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat          >  2:180478/5‑148 (MQ=255)
                                                                                                                                         tGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGa        <  5:252486/149‑1 (MQ=255)
                                                                                                                                             gCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGga                                                                                          >  2:261756/1‑63 (MQ=255)
                                                                                                                                             gCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGga                                                                                          <  1:261756/63‑1 (MQ=255)
                                                                                                                                                ttCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagccagc  >  5:306690/1‑148 (MQ=255)
                                                                                                                                                ttCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagccagc  >  3:167252/1‑148 (MQ=255)
                                                                                                                                                   cAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                  <  2:251552/65‑1 (MQ=255)
                                                                                                                                                   cAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                  >  1:251552/1‑65 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGC  >  NC_002947/3223074‑3223365

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