Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,492,267 T→G 6.7% intergenic (+28/+67) cysNC → / ← imm sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,492,2670TG 6.7% 173.5 / 2.2 56intergenic (+28/+67)cysNC/immsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein
Reads supporting (aligned to +/- strand):  ref base T (29/23);  new base G (2/2);  total (31/25)
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.00e+00

CCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTAC  >  NC_002947/1492129‑1492410
                                                                                                                                          |                                                                                                                                               
ccGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGt                                                                                                                                        <  3:36216/148‑1 (MQ=255)
   ttCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGttttt                                                                                                                                    >  6:136997/1‑149 (MQ=255)
     cGATGTGCCGCTGGATGCTGATCTGGTCATCGACACCAAGGCCACACGGGTGTACTATGGCGTGAAACTGGTGTTGCTCGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGttttcc                                                                                                                                   <  1:280377/148‑3 (MQ=255)
      gATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                  >  6:45994/1‑148 (MQ=255)
       aTGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACACGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTa                                                                                                                                >  6:151156/1‑149 (MQ=255)
          tgCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTAtt                                                                                                                              <  5:136997/148‑1 (MQ=255)
           gCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGc                                                                                                                            >  7:263382/1‑149 (MQ=255)
            ccGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCt                                                                                                                           <  5:45994/149‑1 (MQ=255)
                      cTGATGTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTTCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGTGGTTTTTTCTATTGCTCAGTTCggg                                                                                                                  <  4:180549/148‑1 (MQ=255)
                       tGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGTCGTGAAACAGGTGCTGGACGTGCTCCGTCAGCGTGGCGCGATCTAAGCGCTAGCCAGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTATAATGCTCAGTTCGCGt                                                                                                                 >  1:74965/1‑148 (MQ=255)
                            ttGTGATCGACACCCAGGCCACCACCGTGGACGAAGCCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGAGCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTATTATCTATTGCTCAGTTCGGGAGAAAAg                                                                                                           >  3:194558/1‑149 (MQ=255)
                                   cGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCt                                                                                                      >  2:68375/1‑147 (MQ=255)
                                     acacCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCtgaac                                                                                                  >  2:321668/1‑145 (MQ=255)
                                         ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCgg                                                                                                                   <  3:152217/128‑1 (MQ=255)
                                         ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCgg                                                                                                                   >  4:152217/1‑128 (MQ=255)
                                         ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                               >  3:196452/1‑148 (MQ=255)
                                          cAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                          >  3:48715/1‑104 (MQ=255)
                                          cAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                          <  4:48715/104‑1 (MQ=255)
                                            ggCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAgg                                                                                           <  4:196452/149‑1 (MQ=255)
                                                accaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCt                                                                                       >  2:308691/1‑149 (MQ=255)
                                                     cGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTc                                                                                  >  3:154810/1‑149 (MQ=255)
                                                         gACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgtgt                                                                              <  8:263382/149‑1 (MQ=255)
                                                               ggCGTGATAGAGGTGCTGGTCGTGCTGCGTCAGCGTGGCGCGATCTAAGGGCTAGCCTGCAAGGAATAGCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCTGAGTTCAGGGCCTTGCTCGTGTTAGCt                                                                         <  4:194558/148‑1 (MQ=255)
                                                                gCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTgg                                                                       <  2:165820/149‑1 (MQ=255)
                                                                 cGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTggg                                                                      <  1:68375/149‑1 (MQ=255)
                                                                  gTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGc                                                                     >  8:252847/1‑149 (MQ=255)
                                                                          ggTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAAcc                                                             <  5:151156/149‑1 (MQ=255)
                                                                           gTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCg                                                            >  8:150681/1‑149 (MQ=255)
                                                                                 gACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGt                                                           <  5:130290/144‑1 (MQ=255)
                                                                                 gACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGt                                                           >  6:130290/1‑144 (MQ=255)
                                                                                   cGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTc                                                                                              <  1:87739/107‑1 (MQ=255)
                                                                                   cGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTc                                                                                              >  2:87739/1‑107 (MQ=255)
                                                                                    gTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCg                                                   <  7:150681/149‑1 (MQ=255)
                                                                                            tCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAAc                                                              >  6:243260/1‑130 (MQ=255)
                                                                                            tCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAAc                                                              <  5:243260/130‑1 (MQ=255)
                                                                                                      gcgATCTAAGCGCTAGCCTGCAAGGAAAACCCCACTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGc                                 >  1:73140/1‑149 (MQ=255)
                                                                                                       cgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCt                                <  5:117211/149‑1 (MQ=255)
                                                                                                                 gcTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTc                      >  7:143264/1‑149 (MQ=255)
                                                                                                                 gcTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTc                      <  1:308691/149‑1 (MQ=255)
                                                                                                                    aGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTa                    <  4:154810/148‑1 (MQ=255)
                                                                                                                     gCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAg                   >  8:72577/1‑148 (MQ=255)
                                                                                                                         atagaaaaaaCCCCGCCGAAGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACtt              >  7:110140/7‑149 (MQ=255)
                                                                                                                          cAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCAtt                                   >  6:161095/1‑127 (MQ=255)
                                                                                                                          cAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCAtt                                   <  5:161095/127‑1 (MQ=255)
                                                                                                                           aaGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACtttt            >  3:271493/1‑149 (MQ=255)
                                                                                                                            gaaaaaaCCCCGCCGAAGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTTGCTAGCCGCCATTGCTTGa                             <  4:164609/128‑1 (MQ=255)
                                                                                                                            gaaaaaaCCCCGCCGAAGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTTGCTAGCCGCCATTGCTTGa                             >  3:164609/4‑131 (MQ=255)
                                                                                                                            aGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTa                                            <  6:340854/116‑1 (MQ=255)
                                                                                                                            aGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTa                                            >  5:340854/1‑116 (MQ=255)
                                                                                                                             ggAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTc          <  7:252847/149‑1 (MQ=255)
                                                                                                                              gAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttttt           >  6:87006/1‑147 (MQ=255)
                                                                                                                              gAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttttt           >  7:22977/1‑147 (MQ=255)
                                                                                                                                   ccccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTc          <  5:245397/143‑1 (MQ=255)
                                                                                                                                   ccccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTc          >  6:245397/1‑143 (MQ=255)
                                                                                                                                     ccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGCCTATTTCTAGTACTTTTTCAGGTTTAc  >  2:53063/1‑149 (MQ=255)
                                                                                                                                     ccGCTTCGGCGGGGTTTTTTCTAATGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTa   >  7:162901/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                               
CCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTAC  >  NC_002947/1492129‑1492410

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: