Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,110,166 G→A 7.2% R1374C (CGC→TGC)  PP_1882 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,110,1660GA 7.2% 163.2 / 2.3 55R1374C (CGC→TGC) PP_1882membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

GTCTGCCAGCGATCGGCTTTGCCCGCTGTAGCCGATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCGCCATAGGCGGTATAACGGAC  >  NC_002947/2110020‑2110307
                                                                                                                                                  |                                                                                                                                             
gTCTGCCAGCGATCGGCTTTGCCCGCTGTAGCCGATGCTGTATTGACCGGTGGGGTCGCTGAAGGTTATTGGGTTGCCTTGGCAATAGCGGTAGTAGTTCAGCCCCCACGCGGCGAAGGGTCTGAGTGCGTCGGGGGTGTGAATgcgct                                                                                                                                             >  5:209415/1‑148 (MQ=255)
   tGCCAGCGATCGGCTTTGCCCGCTGTAGCCGATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCAt                                                                                                                                            >  3:54402/1‑147 (MQ=255)
     ccAGCGATCGGCTTTGCCCGCTGTAGCCGATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAg                                                                                                                                         >  8:217119/1‑148 (MQ=255)
                     cccGCTGTAGCCGATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCggg                                                                                                                                  <  6:23410/139‑1 (MQ=255)
                     cccGCTGTAGCCGATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCggg                                                                                                                                  >  5:23410/1‑139 (MQ=255)
                        gCTGTAGCCGATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGt                                                                                                                     <  4:158903/149‑1 (MQ=255)
                                 gATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTAc                                                                                                            >  7:44927/1‑149 (MQ=255)
                                        tAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTGGGCGCGGTAGCCGCTACCGAGCa                                                                                                      <  7:217119/148‑1 (MQ=255)
                                        tAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCa                                                                                                      <  8:44927/148‑1 (MQ=255)
                                        tAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGGCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAg                                                                                                     <  1:35213/149‑1 (MQ=255)
                                             cccGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTAc                                                                                                 >  2:135940/1‑148 (MQ=255)
                                              ccGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGTCAATAGCCTTAGTAGTTAAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTATGCGCGGTAGCCGCTACCGAGCAGGTACCa                                                                                               <  4:186906/149‑1 (MQ=255)
                                               cGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATc                                                                                                                                           <  4:136718/104‑1 (MQ=255)
                                               cGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATc                                                                                                                                           >  3:136718/1‑104 (MQ=255)
                                                     ggtcgCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAgcgc                                                                                                                                              <  7:84341/95‑1 (MQ=255)
                                                     ggtcgCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAgcgc                                                                                                                                              >  8:84341/1‑95 (MQ=255)
                                                      gtcgCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGc                                                                                                                  <  7:50496/122‑1 (MQ=255)
                                                      gtcgCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGc                                                                                                                  >  8:50496/1‑122 (MQ=255)
                                                          cgGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCgg                                                                                    >  7:169908/1‑148 (MQ=255)
                                                               ggTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGa                                                                              <  1:200205/149‑1 (MQ=255)
                                                                   aTCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCgg                                                                                    >  8:127035/1‑139 (MQ=255)
                                                                   aTCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCgg                                                                                    <  7:127035/139‑1 (MQ=255)
                                                                        gTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGc                                                                     <  3:95331/149‑1 (MQ=255)
                                                                        gTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGc                                                                     >  4:95331/1‑149 (MQ=255)
                                                                              ttGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTg                                                                >  8:88257/1‑148 (MQ=255)
                                                                              ttGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTg                                                                >  7:191921/1‑148 (MQ=255)
                                                                                ggCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGt                                                               <  4:54402/147‑1 (MQ=255)
                                                                                ggCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTa                                                              <  1:135940/148‑1 (MQ=255)
                                                                                     tAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGcc                                                                    >  2:6629/1‑137 (MQ=255)
                                                                                     tAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCACCCGCTGTCGGGATCGAGCGCTTCGcc                                                                    <  1:6629/137‑1 (MQ=255)
                                                                                       gCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAggg                                                      >  3:204728/1‑149 (MQ=255)
                                                                                          gtagtaGTTCAGCCCCCCCGCGCCGACGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGc                                                   <  8:169908/149‑1 (MQ=255)
                                                                                          gtagtaGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGc                                                   >  8:92639/1‑149 (MQ=255)
                                                                                          gtagtaGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGc                                                   <  2:32235/149‑1 (MQ=255)
                                                                                                   cAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGt                                           <  5:99604/148‑1 (MQ=255)
                                                                                                       ccccccGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGAtc                                      <  8:191921/149‑1 (MQ=255)
                                                                                                        cccccGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAg                                                                                                     >  7:113567/1‑85 (MQ=255)
                                                                                                        cccccGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAg                                                                                                     <  8:113567/85‑1 (MQ=255)
                                                                                                                 cGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAgc                                                                            >  5:46743/1‑101 (MQ=255)
                                                                                                                 cGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAgc                                                                            <  6:46743/101‑1 (MQ=255)
                                                                                                                 cGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGc                            >  2:242438/1‑149 (MQ=255)
                                                                                                                 cGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGc                            >  7:246338/1‑149 (MQ=255)
                                                                                                                  gAAGGGGCAGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCAACGGGTTGTAGGCGCGGTAGCCACTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCTTGcc                           >  2:161566/1‑149 (MQ=255)
                                                                                                                   aaGGGGCTGAGTGCGCCGGGGCTGTGAAAGCACATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCa                                             >  3:234084/1‑130 (MQ=255)
                                                                                                                   aaGGGGCTGAGTGCGCCGGGGCTGTGAAAGCACATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCa                                             >  3:237918/1‑130 (MQ=255)
                                                                                                                   aaGGGGCTGAGTGCGCCGGGGCTGTGAAAGCACATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCa                                             <  4:234084/130‑1 (MQ=255)
                                                                                                                   aaGGGGCTGAGTGCGCCGGGGCTGTGAAAGCACATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCa                                             <  4:237918/130‑1 (MQ=255)
                                                                                                                       ggCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCg                      >  4:115246/1‑149 (MQ=255)
                                                                                                                       ggCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCg                      <  4:204728/149‑1 (MQ=255)
                                                                                                                        gCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCGc                     <  7:92639/149‑1 (MQ=255)
                                                                                                                                       gCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGcc                                                  >  2:251279/1‑105 (MQ=255)
                                                                                                                                       gCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGcc                                                  <  1:251279/105‑1 (MQ=255)
                                                                                                                                       gCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCGCCATAGGCGGTATAac      >  2:121145/1‑149 (MQ=255)
                                                                                                                                       gCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCGCCATAGGCGGTATAac      <  1:242438/149‑1 (MQ=255)
                                                                                                                                       gCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGAGCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCGCCATAGGCGGTATAac      >  5:193925/1‑149 (MQ=255)
                                                                                                                                           tgAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCGCCATAGGCGGTATAacggac  >  2:45011/1‑149 (MQ=255)
                                                                                                                                            gAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCg                                                                                     <  7:60980/65‑1 (MQ=255)
                                                                                                                                            gAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCg                                                                                     >  8:60980/1‑65 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GTCTGCCAGCGATCGGCTTTGCCCGCTGTAGCCGATGCTGTAATGCCCGGTCGGGTCGCGGAAGGTTATCGGGTTGCCTTGGCAATAGCCGTAGTAGTTCAGCCCCCCCGCGCCGAAGGGGCTGAGTGCGTCGGGGCTGTGAAAGCGCATCAGCACCGGGTTGTAGGCGCGGTAGCCGCTACCGAGCAGGTACCAGCCGCTGTCGGGATCGAGCGCTTCGCCGTTGTAGCCCAGGGTGCCGAGCAGTGGATCGTCGGCGTGCCGCTCGCCATAGGCGGTATAACGGAC  >  NC_002947/2110020‑2110307

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