Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,855,020 G→T 7.3% A27A (GCC→GCA dnaX ← DNA polymerase III subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,855,0200GT 7.3% 170.7 / 3.0 55A27A (GCC→GCAdnaXDNA polymerase III subunit gamma
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base T (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 = 7.83e-01

AACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCAC  >  NC_002947/4854872‑4855168
                                                                                                                                                    |                                                                                                                                                    
aaCAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAg                                                                                                                                                      <  5:365979/149‑1 (MQ=255)
   aGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGcc                                                                                                                                                   >  6:223275/1‑149 (MQ=255)
     gTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCCTTGATCAATGCCtt                                                                                                                                                 <  5:223275/149‑1 (MQ=255)
     gTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCtt                                                                                                                                                 <  3:75998/149‑1 (MQ=255)
         cacaGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAg                                                                                                                                              >  8:278385/1‑148 (MQ=255)
             gggCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACAt                                                                                                                                         <  7:280172/149‑1 (MQ=255)
                  tGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTc                                                                                                                                    >  1:122168/1‑149 (MQ=255)
                            gCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCAAAGGCTTTGATCAAGGCATAGAGCACATGGTTCTGGCCGACCa                                                                                                                          >  8:121108/1‑149 (MQ=255)
                                gTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTc                                                                                                                      >  5:207702/1‑149 (MQ=255)
                                    cacaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTcgcg                                                                                                                   <  1:47040/148‑1 (MQ=255)
                                     acaGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCgg                                                                                                                  <  8:282849/148‑1 (MQ=255)
                                             ggCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAgcgc                                                                                                         <  3:1382/149‑1 (MQ=255)
                                                 cTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGAc                                                                                                     <  7:7430/149‑1 (MQ=255)
                                                  ttGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGCCCATTTCGCGGAACGAGCGCGGACg                                                                                                    >  2:304238/1‑149 (MQ=255)
                                                  ttGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACg                                                                                                    >  5:290923/1‑149 (MQ=255)
                                                        ataataCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCAttt                                                                                              >  5:332605/1‑149 (MQ=255)
                                                            taCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGt                                                                                          <  2:246119/149‑1 (MQ=255)
                                                            taCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGt                                                                                          <  8:21627/149‑1 (MQ=255)
                                                            taCGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGt                                                                                          >  7:21627/1‑149 (MQ=255)
                                                              cGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTg                                                                                         <  2:122168/148‑1 (MQ=255)
                                                                    aTGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAAc                                                                                  <  6:207702/149‑1 (MQ=255)
                                                                        tcgtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGa                                                                              <  7:278385/149‑1 (MQ=255)
                                                                          gtcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGAt                                                                             >  1:353468/1‑148 (MQ=255)
                                                                           tcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATaa                                                                           >  4:267988/1‑149 (MQ=255)
                                                                           tcTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTACAAGAACCTGATaa                                                                           >  1:128342/1‑149 (MQ=255)
                                                                                 ccACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATc                                                                     <  2:353468/149‑1 (MQ=255)
                                                                                       cgcgGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGaaaaa                                                               >  2:339942/1‑149 (MQ=255)
                                                                                           ggTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCAt                                                           >  4:303305/1‑149 (MQ=255)
                                                                                                 ggTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGc                                                     <  6:207689/149‑1 (MQ=255)
                                                                                                  gTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGcc                                                    >  4:165463/1‑149 (MQ=255)
                                                                                                           tAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGc                                           <  6:332605/149‑1 (MQ=255)
                                                                                                             ggCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCgg                                         <  3:267988/149‑1 (MQ=255)
                                                                                                               cATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGaa                                       <  3:303305/149‑1 (MQ=255)
                                                                                                                aTGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGAt                                                                             <  2:306744/110‑1 (MQ=255)
                                                                                                                   gTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTACAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATc                                   <  2:128342/149‑1 (MQ=255)
                                                                                                                      cAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCTTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGc                                 >  4:87426/1‑148 (MQ=255)
                                                                                                                      cAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGc                                 <  6:376917/148‑1 (MQ=255)
                                                                                                                        ggCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTaa                              >  4:351125/1‑149 (MQ=255)
                                                                                                                               gTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGCCGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagc                       >  8:279319/1‑149 (MQ=255)
                                                                                                                                ttGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcg                      <  7:121108/149‑1 (MQ=255)
                                                                                                                                 tgctcAAGGCCTTGATCAATGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGc                                 >  6:156530/5‑137 (MQ=255)
                                                                                                                                 tgctcAAGGCCTTGATCAATGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGc                                 <  5:156530/133‑1 (MQ=255)
                                                                                                                                  gTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcg                      >  2:318741/1‑147 (MQ=255)
                                                                                                                                   tCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACGCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcggag                   >  6:91569/1‑149 (MQ=255)
                                                                                                                                      aaGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcggagcgg                >  6:265980/1‑149 (MQ=255)
                                                                                                                                       aGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCggg               <  3:351125/149‑1 (MQ=255)
                                                                                                                                        ggCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCTCGTAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTTATCGAGCGGAAGATCGCTAATGCTAGCGGAGCgggg              <  7:174956/149‑1 (MQ=255)
                                                                                                                                           aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcggagcgg                >  7:29996/1‑144 (MQ=255)
                                                                                                                                           aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTagcggagcgg                <  8:29996/144‑1 (MQ=255)
                                                                                                                                           aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGc            >  2:112784/1‑148 (MQ=255)
                                                                                                                                           aTTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCCGAGCGGGGGc            >  2:316349/1‑148 (MQ=255)
                                                                                                                                            ttGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCa           >  3:181459/1‑148 (MQ=255)
                                                                                                                                               aTCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAAt       >  3:258249/1‑149 (MQ=255)
                                                                                                                                                 cAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTg     >  5:309633/1‑149 (MQ=255)
                                                                                                                                                   aGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCa   <  4:181459/149‑1 (MQ=255)
                                                                                                                                                   aGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCa   <  1:112784/149‑1 (MQ=255)
                                                                                                                                                    ggCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCAc  >  1:79658/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACAGGTGCCACAGGGCGTGGAGGTGATGCCGGTTTCACAGTTCAGGCACTTGGCGATAATACGGGCGATGGTCGTCTTGCCCACCCCGCGGGTACCGGTGAACAGGTAGGCATGGTGCAGGCGCTGGTTGTCCAAGGCATTGATCAAGGCCTTGAGCACATGGGTCTGGCCGACCATTTCGCGGAACGAGCGCGGACGCCATTTTCGTGCAAGAACCTGATAACTCATCGAAAAACCATCGTAGCCTGATCGAGCGGAAGATCGCTAATGCTAGCGGAGCGGGGGCAAAATTGCAC  >  NC_002947/4854872‑4855168

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: