Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,935,581 C→G 8.4% intergenic (‑169/+109) PP_3467 ← / ← PP_3468 sigma‑54 dependent transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,935,5810CG 8.4% 176.1 / 3.3 56intergenic (‑169/+109)PP_3467/PP_3468sigma‑54 dependent transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base G (2/3);  total (28/28)
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

TGACGGGGTACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCCCAGCCAGTTGC  >  NC_002947/3935440‑3935726
                                                                                                                                             |                                                                                                                                                 
tgACGGGGTACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCaa                                                                                                                                            >  2:256026/1‑149 (MQ=255)
tgACGGGGTACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCaa                                                                                                                                            >  5:35164/1‑149 (MQ=255)
 gACGGGGTACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTGGCACCAAc                                                                                                                                           >  3:148037/1‑149 (MQ=255)
      ggTACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGa                                                                                                                                      <  2:50635/149‑1 (MQ=255)
       gTACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTGGTGCGAACCCTggc                                                                                                                                     <  7:144217/149‑3 (MQ=255)
        tACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAgg                                                                                                                                    <  8:108664/149‑1 (MQ=255)
              cTTGGGAATTTTCCGCAGTGTGCGGGTGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTgcg                                                                                                                              <  3:224987/149‑1 (MQ=255)
                  ggAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTg                                                                                                                           <  8:73666/148‑1 (MQ=255)
                    aaTTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTAc                                                                                                                        <  3:246675/149‑1 (MQ=255)
                        ttCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCtg                                                                                                                     <  1:256026/148‑1 (MQ=255)
                          ccGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTgg                                                                                                                  >  6:131932/1‑149 (MQ=255)
                           cGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTgg                                                                                                                  <  8:313131/148‑1 (MQ=255)
                                tgtgCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCgg                                                                                                            <  3:251932/149‑1 (MQ=255)
                                  tgCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCggg                                                                                                           <  4:254412/148‑1 (MQ=255)
                                  tgCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGGGGAAGCTGTCTTGCCAGGGTTCGAACCAACCCTGAGGACTGCGCTGTACCTGGTGTAGCGGGc                                                                                                          <  8:208158/149‑1 (MQ=255)
                                     gggAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAg                                                                                                       >  4:199394/1‑149 (MQ=255)
                                      ggAGAAATGGCTGGGCGGCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTGGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGc                                                                                                      <  2:97386/149‑1 (MQ=255)
                                           aaTGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTCCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCcgc                                                                                                  >  4:41875/1‑148 (MQ=255)
                                               gCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAaga                                                                                             <  1:53393/149‑1 (MQ=255)
                                                cTGGGCGTCATCTGGTAGGGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAagag                                                                                            <  3:303675/149‑1 (MQ=255)
                                                   ggCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGcc                                                                                         >  6:156603/1‑149 (MQ=255)
                                                         aTCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGa                                                                                    <  5:156603/148‑1 (MQ=255)
                                                          tCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGa                                                                                                        <  6:268510/127‑1 (MQ=255)
                                                          tCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGa                                                                                                        >  5:268510/1‑127 (MQ=255)
                                                           cTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGt                                                                                 >  8:103492/1‑149 (MQ=255)
                                                              gTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTCCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGt                                                                                 <  3:41875/146‑1 (MQ=255)
                                                                    gatggatTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAgg                                                                         >  2:85544/1‑148 (MQ=255)
                                                                    gatggatTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAgg                                                                         <  6:35164/148‑1 (MQ=255)
                                                                    gatggatTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACAAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAgg                                                                         >  5:238395/1‑148 (MQ=255)
                                                                      tggatTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGa                                                                       >  6:148903/1‑148 (MQ=255)
                                                                          ttGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTgcg                                                                                                                              >  4:151944/1‑89 (MQ=255)
                                                                          ttGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTgcg                                                                                                                              <  3:151944/89‑1 (MQ=255)
                                                                          ttGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATatc                                                                   >  5:173701/1‑148 (MQ=255)
                                                                             gggCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATatcgatc                                                               <  3:199394/149‑1 (MQ=255)
                                                                              tgCGTAGGTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGGCTGGGTTCGCACCAACCCTGAGGAGTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATatcgatcg                                                              <  5:131932/148‑1 (MQ=255)
                                                                               gCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATAtcgatcga                                                             <  1:85544/149‑1 (MQ=255)
                                                                                    gCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGat                                                                      <  6:54615/135‑1 (MQ=255)
                                                                                    gCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGat                                                                      >  5:54615/1‑135 (MQ=255)
                                                                                          cAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATa                                                  <  6:235097/149‑1 (MQ=255)
                                                                                                           gcgcAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAgg                                  >  6:140065/1‑148 (MQ=255)
                                                                                                           gcgcAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAg                                   <  5:148903/147‑1 (MQ=255)
                                                                                                              cAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGcccc                                        >  7:1739/1‑139 (MQ=255)
                                                                                                              cAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGcccc                                        <  8:1739/139‑1 (MQ=255)
                                                                                                               aaTCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTc                             <  7:103492/149‑1 (MQ=255)
                                                                                                                  cccTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCAGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCTATAGCCAAGCCCCGTCAGGCACTCGAc                          >  6:85096/1‑149 (MQ=255)
                                                                                                                   ccTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGAcg                         >  3:110408/1‑149 (MQ=255)
                                                                                                                   ccTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGAcg                         >  2:99906/1‑149 (MQ=255)
                                                                                                                      gtgGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGAcgccg                      >  3:166526/1‑149 (MQ=255)
                                                                                                                                 cTTGCCAGGGTTGGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCcca           <  4:148037/149‑1 (MQ=255)
                                                                                                                                  ttGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCccag          >  6:50016/1‑149 (MQ=255)
                                                                                                                                    gCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCccagc         >  5:304994/1‑148 (MQ=255)
                                                                                                                                     tcAGGGTTGGTGCGAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCccagc         >  3:298213/2‑147 (MQ=255)
                                                                                                                                     ccAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCccagc         >  7:202478/1‑147 (MQ=255)
                                                                                                                                      cAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCccagccag      >  2:60654/1‑149 (MQ=255)
                                                                                                                                          gTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCCCAGCCAGTtgc  >  8:255093/1‑149 (MQ=255)
                                                                                                                                          gTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCCCAGCCAGTtgc  >  3:32295/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
TGACGGGGTACGGGCTTGGGAATTTTCCGCAGTGTGCGGGAGAAATGGCTGGGCGTCATCTGGTAGTGGATGGATTGGGGCGTAGCTGGCCAGGTTAATTGAAGCTGGCGCAATCCCTGTGGAAGCTGGCTTGCCAGGGTTCGCACCAACCCTGAGGACTGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCCTCGAGGTCAATCAGGGATATCGATCGATATCACCAATAGCCAAGCCCCGTCAGGCACTCGACGCCGCCTTCCGCCCAGCCAGTTGC  >  NC_002947/3935440‑3935726

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: