Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,741,014 C→T 6.9% G307S (GGC→AGC)  PP_3305 ← TerC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,741,0140CT 6.9% 168.9 / 2.8 58G307S (GGC→AGC) PP_3305TerC family membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base T (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

GGCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAATGCAAAATAGAGC  >  NC_002947/3740874‑3741162
                                                                                                                                            |                                                                                                                                                    
ggCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTaa                                                                                                                                              >  7:129809/1‑149 (MQ=255)
ggCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTaa                                                                                                                                              <  7:225299/149‑1 (MQ=255)
   gACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGc                                                                                                                                           >  8:256538/1‑149 (MQ=255)
   gACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGc                                                                                                                                           >  2:241992/1‑149 (MQ=255)
      cGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGTGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcg                                                                                                                                              <  8:270505/143‑4 (MQ=255)
      cGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcggaga                                                                                                                                          >  7:270505/1‑140 (MQ=255)
         tttCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcg                                                                                                                                              >  1:130186/1‑137 (MQ=255)
         tttCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcg                                                                                                                                              <  2:130186/140‑4 (MQ=255)
           tCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGa                                                                                                                                    >  8:133070/1‑148 (MQ=255)
               gCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAAt                                                                                                                                             <  8:36333/135‑1 (MQ=255)
               gCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAAt                                                                                                                                             >  7:36333/1‑135 (MQ=255)
                    gAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGCCCCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCATACAGCAGGCCGATTGTAATCATCAACGATATTTCgcat                                                                                                                          >  2:10146/1‑146 (MQ=255)
                    gAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCgg                                                                                                                            <  6:15221/147‑1 (MQ=255)
                      aTCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTg                                                                                                                         <  7:133070/148‑1 (MQ=255)
                            ggCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACtt                                                                                                                   <  3:55893/148‑1 (MQ=255)
                               gATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTAcc                                                                                                               <  3:144402/149‑1 (MQ=255)
                               gATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTAcc                                                                                                               <  8:129809/149‑1 (MQ=255)
                                     gTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAAt                                                                                                         >  6:283592/1‑149 (MQ=255)
                                          cgcgCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCca                                                                                                             >  3:158371/1‑140 (MQ=255)
                                          cgcgCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCca                                                                                                             <  4:158371/140‑1 (MQ=255)
                                           gcgcTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCAt                                                                                                    >  4:251052/1‑148 (MQ=255)
                                             gcTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGaa                                                                                                 <  3:251052/149‑1 (MQ=255)
                                                   gcagcaGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACACACAATGTAATGAAGGaaaa                                                                                           >  8:268308/1‑149 (MQ=255)
                                                        agccagaGTCCGTAATTATTCAAGAACCATTAGCGGAAGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCATCAGGCCGATTGTAATCCTCAACTATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTc                                                                                      <  2:220786/149‑1 (MQ=255)
                                                          ccagcGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTGTACCCACAATGTCATGAAGGAAAATTTTc                                                                                      >  6:35627/1‑147 (MQ=255)
                                                             gcGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGcc                                                                                 >  1:189589/1‑149 (MQ=255)
                                                                gCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGAt                                                                              <  5:31028/149‑1 (MQ=255)
                                                                gCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGAt                                                                              <  5:35755/149‑1 (MQ=255)
                                                                       ttattCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACa                                                                       <  7:256538/149‑1 (MQ=255)
                                                                             aagaagCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCAcc                                                                 <  2:61361/149‑1 (MQ=255)
                                                                             aagaagCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCAc                                                                  >  3:275309/1‑148 (MQ=255)
                                                                             aagaagCATTAGCGGACGCTAGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCAcc                                                                 <  7:39615/149‑1 (MQ=255)
                                                                                     ttAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCATCGCCAg                                                         >  7:124818/1‑149 (MQ=255)
                                                                                           gACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCata                                                   <  1:10146/149‑1 (MQ=255)
                                                                                            aCGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCatat                                                  <  4:275309/149‑1 (MQ=255)
                                                                                            aCGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCatat                                                  >  6:250089/1‑149 (MQ=255)
                                                                                                tCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATAttt                                                >  5:293595/1‑147 (MQ=255)
                                                                                                tCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGAATAttt                                                >  8:129565/1‑147 (MQ=255)
                                                                                                  gggTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAgg                                            >  1:225319/1‑149 (MQ=255)
                                                                                                   ggTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGt                                           <  5:283592/149‑1 (MQ=255)
                                                                                                    gatttCCAAACCGAAAGCAACACCCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTa                                          <  8:124818/147‑1 (MQ=255)
                                                                                                    gTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTa                                          <  6:293595/149‑1 (MQ=255)
                                                                                                     ttttCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTa                                          >  2:10479/1‑148 (MQ=255)
                                                                                                           aaaGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAg                                   >  2:232451/1‑149 (MQ=255)
                                                                                                              gCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCa                                                                   <  8:197256/114‑1 (MQ=255)
                                                                                                              gCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCa                                                                   >  7:197256/1‑114 (MQ=255)
                                                                                                                gaaagaaaCAAGCCGCCAATCAGCAGGCCTATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGTAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGAc                               <  7:268308/148‑1 (MQ=255)
                                                                                                                gAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACt                              >  3:258913/1‑149 (MQ=255)
                                                                                                                   aGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGa                            >  3:261978/1‑148 (MQ=255)
                                                                                                                      aaCACGCCTCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGTCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAgg                                            <  5:75539/129‑1 (MQ=255)
                                                                                                                      aaCACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCAt                        <  6:24108/149‑1 (MQ=255)
                                                                                                                      aaCACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCAt                        >  8:82898/1‑149 (MQ=255)
                                                                                                                       acacGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATg                       >  5:108888/1‑149 (MQ=255)
                                                                                                                        cacGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAAc                                       >  6:51486/1‑132 (MQ=255)
                                                                                                                        cacGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAAc                                       <  5:51486/132‑1 (MQ=255)
                                                                                                                                cAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAAt              <  2:225319/149‑1 (MQ=255)
                                                                                                                                     agcagGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAg                                   <  6:15762/123‑1 (MQ=255)
                                                                                                                                     agcagGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAATGCaa          <  6:108888/148‑1 (MQ=255)
                                                                                                                                            cGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAATGCAAAATAGAgc  >  2:266880/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GGCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAATGCAAAATAGAGC  >  NC_002947/3740874‑3741162

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: