Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,443 C→G intergenic (+577/+179) rtcR → / ← glpG DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4430CG97.1% 88.4 / ‑5.3 34intergenic (+577/+179)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (11/22);  minor base T (0/1);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCG  >  NC_000913/3560208‑3560688
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
ttcttTTCCTATCAACGAGTTAGCAACGTAGAGCTTGTTGGTGACAATGTGCAAATTGCTGTGATTGCGAAGTGCGTGCGCTACAGCTTCCGGCGTGGTGCCGATATCGCTCAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGCTGCGCTCTTTTTCTTCGGCCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCTGCGCAGCAGCGCAATGACGGcgc                                                                                                                                                                                                                                          >  2:168148/1‑248 (MQ=255)
       cccATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAg                                                                                                                                                                                                                                  <  1:93996/250‑1 (MQ=255)
          aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGTGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCATCGCCATGATGGCGCAGGATCAGGtt                                                                                                                                                                                                                               <  2:127690/251‑1 (MQ=255)
           acAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGAGGGCGCAGGATCAGGttt                                                                                                                                                                                                                              <  2:65167/250‑1 (MQ=255)
                 tgtTAGCAACGTTGAGAGTGTTGGTGACAATGCGCAAATTGCTCTGATTGAGCACAGCGTGCGCTACCGCGTCCAGCGTGGTACCGATATCGATAAACAGCGTCGCGCCATTGGGGATTTGCTCCACCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGTCGCAGGATCAGGTTTTGcgcc                                                                                                                                                                                                                        >  2:43705/1‑248 (MQ=255)
                               ggaTTGTTGGGGACACTGCGCAGATTGCTGTGAGTGAGCCGTGCGTGCTCTACCGCCTCCGGCGTGGCGCCGTTAACGTTCACCTTAGTCTTGCCATTCTGGATTTGCTTCTCCTCTTTGCGGGCGTTGCGCTCTTTTTCTTCTGTCTGTGTGGCCTGGCGATCGTGCCACGGCGTTTTAACCGATCTGGAAGGCTTCGCCGCAGCTCCATGTTGTCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGat                                                                                                                                                                                                          <  1:54927/250‑2 (MQ=255)
                               agaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTGGCTTCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                                                                                                                          <  1:54924/251‑1 (MQ=255)
                                gaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                                                                                                                                         <  1:15788/251‑1 (MQ=255)
                                 aTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGATACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGACTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGt                                                                                                                                                                                                        >  2:105640/1‑251 (MQ=255)
                                             aaTGCGCAAAATGCGGTGATTGAGCGGTGCGTGCCCTACCGCTTCCGGCTTGGGGCCGATAGCGATAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTTCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGTCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                                                                                                                            <  2:332157/251‑1 (MQ=255)
                                                    aaaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                     >  1:312713/1‑251 (MQ=255)
                                                           tgtgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCAGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                              >  2:328564/1‑251 (MQ=255)
                                                                                        gcgggCGTGGGCCCGATATCCATAAACACTGGCGAGCCCTTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGTCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGCGGCGAATAGTGTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGtcc                                                                                                                                                  <  2:288723/247‑3 (MQ=255)
                                                                                          cGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                                                                                                                >  1:223136/1‑250 (MQ=255)
                                                                                                   gCCGATATCGAAAACCAGCGGCGAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGAGGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACCct                                                                                                                                      <  2:43731/251‑1 (MQ=255)
                                                                                                                   agtgTCGAGCCATTGGGGATTGGTTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                      <  2:223136/248‑1 (MQ=255)
                                                                                                                   aGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                      <  1:158644/251‑1 (MQ=255)
                                                                                                                    gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                                                                                                     <  1:279426/251‑1 (MQ=255)
                                                                                                                      gTCGTGCCATTGGGGATTTGCCCCGCCACTTTGCGGGCGAGGCGCGCTTTTTCTGCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCTCCGCAGCGCCATGATGGCGCAGGATCGGGTTTTGCTCCGCCTGCTCATTGAGTTCGCGGCGAATATTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                   <  1:200622/251‑1 (MQ=255)
                                                                                                                         gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                                                                                                                             >  1:328192/1‑174 (MQ=255)
                                                                                                                         gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                                                                                                                             <  2:328192/174‑1 (MQ=255)
                                                                                                                            ccattgtggctttggaacgcaaCTTTGCGGGCGATGCGATCTTTTTCTTCGGTCTGGGTCGACTTGCGATCGTGCCCTGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGAGAATAATCTGCGGGCAGACGGCGAACTGCTCCACCAGCTCTTAGGAAATGACATAACCCTGCTGTTTAACAGGTTCGATAATAcc                                                                                                             >  2:25977/22‑251 (MQ=255)
                                                                                                                                    gATTTGCTCCGCCACTTTGCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCAGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTTACATAACCCTGCTGTTTAAACAGATCGCTAATACCGTTGAGAc                                                                                                     >  2:270511/1‑251 (MQ=255)
                                                                                                                                          cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg                                                                                                         >  2:153392/1‑241 (MQ=255)
                                                                                                                                          cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg                                                                                                         <  1:153392/241‑1 (MQ=255)
                                                                                                                                             cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                            <  1:180735/251‑1 (MQ=255)
                                                                                                                                              gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                           <  1:146546/251‑1 (MQ=255)
                                                                                                                                                 ccTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGGCGTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCTCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt                                                                                        <  2:272493/250‑1 (MQ=255)
                                                                                                                                                               gcgcTCTTTTTCGTCGGTCGGGGTGGCCTTGCGATCGTGCCGCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                           <  2:141756/250‑1 (MQ=255)
                                                                                                                                                                                    ggTGGCCGTGCGATCGTGCCACGGCGTGTTATCCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTTAGTTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGTTCGATAATAACGTTGTGACGTTGTGTTTGTTTCATTTATAAATACCTGGAATTATTTTCGTTTTcgc                                                     <  1:199758/251‑1 (MQ=255)
                                                                                                                                                                                            ggCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTg                                               <  2:123588/248‑1 (MQ=255)
                                                                                                                                                                                                  cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCTATAATTCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTGTCGCGCATTGAGCGAAt                                        <  1:219787/250‑1 (MQ=255)
                                                                                                                                                                                                                           ggAATGCAGCGCCGCCGCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCAc               <  1:43705/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cgcAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGACATCGCTAAACCCACGGCTAACCCGGCg  >  1:197317/1‑251 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCG  >  NC_000913/3560208‑3560688

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