Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 211,621 G→C 9.9% V5709V (GTG→GTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947211,6210GC 9.9% 173.1 / 10.5 61V5709V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base C (3/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

GACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGT  >  NC_002947/211472‑211770
                                                                                                                                                     |                                                                                                                                                     
gacgCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTg                                                                                                                                                       >  1:104796/1‑150 (MQ=255)
       gCCAGAGCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAg                                                                                                                                                >  2:1350703/1‑150 (MQ=255)
        ccAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAgg                                                                                                                                               >  1:1178141/1‑150 (MQ=255)
        ccAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAgg                                                                                                                                               >  1:991788/1‑150 (MQ=255)
          aGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggtg                                                                                                                                             <  1:262529/150‑1 (MQ=255)
           gATCAGCCTGGGCATCAACTCGGCGGAGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggtgg                                                                                                                                            <  2:184984/150‑1 (MQ=255)
            aTCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggtggt                                                                                                                                           <  2:1166468/150‑1 (MQ=255)
              cAGCCTGGGCATCAAATCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTga                                                                                                                                         <  2:991788/150‑1 (MQ=255)
                  cTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATc                                                                                                                                     >  1:764348/1‑150 (MQ=255)
                    gggCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCAc                                                                                                                                   <  2:613886/150‑1 (MQ=255)
                    gggCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGACTTCGGTGGCCGAAGGTTGTTAGATCAc                                                                                                                                   <  2:127721/150‑1 (MQ=255)
                       cATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTa                                                                                                                                >  1:1066156/1‑150 (MQ=255)
                        aTCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTac                                                                                                                               >  2:686690/1‑150 (MQ=255)
                        aTCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTCACCGAAg                                                                                                                                                <  2:1426297/133‑1 (MQ=255)
                        aTCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTCACCGAAg                                                                                                                                                >  1:1426297/1‑133 (MQ=255)
                                          gCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACcaac                                                                                                             >  1:695166/1‑150 (MQ=255)
                                                   gCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGCTCAAGGGACCGAAACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAATGTGGTGAAATAACCTACACCATCACCCTGACCAACAAAAATGGt                                                                                                    >  2:62074/1‑150 (MQ=17)
                                                    cGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTc                                                                                                   <  1:686690/150‑1 (MQ=255)
                                                       ccTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCAt                                                                                                                          >  1:131750/1‑124 (MQ=255)
                                                       ccTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCAt                                                                                                                          <  2:131750/124‑1 (MQ=255)
                                                                 ccTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                <  1:45860/140‑1 (MQ=255)
                                                                 ccTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                >  2:45860/1‑140 (MQ=255)
                                                                    cGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACca                                                                                   >  2:351639/1‑150 (MQ=255)
                                                                       gCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTAAACCATCCCCCTGACCACCAAAGATGGTCTGCCGATCAACAACCACAg                                                                                <  2:764348/150‑1 (MQ=255)
                                                                        cTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGt                                                                               <  2:1248345/150‑1 (MQ=255)
                                                                           ggCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCa                                                                            <  1:1352395/150‑1 (MQ=255)
                                                                            gCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCaacaac                                                                                     <  1:720672/140‑1 (MQ=255)
                                                                            gCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCaacaac                                                                                     >  2:720672/1‑140 (MQ=255)
                                                                                tGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACcaca                                                                                 >  1:704639/1‑140 (MQ=255)
                                                                                tGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACcaca                                                                                 <  2:704639/140‑1 (MQ=255)
                                                                                         ggTTCAAGTGAACGATACCAACGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTAACAACCAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCt                                                              <  2:200720/150‑1 (MQ=17)
                                                                                          gTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTg                                                             >  1:886995/1‑150 (MQ=255)
                                                                                             cAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGc                                                          <  2:1178141/150‑1 (MQ=255)
                                                                                                  gACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCAACCTGAGCGACgg                                                     >  1:1124957/1‑150 (MQ=255)
                                                                                                     cGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACcaa                                                                                                              <  2:42554/90‑1 (MQ=255)
                                                                                                     cGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACcaa                                                                                                              >  1:42554/1‑90 (MQ=255)
                                                                                                           caccGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGt                                            >  2:648138/1‑150 (MQ=255)
                                                                                                            accGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCa                                                                                                                        >  2:114853/1‑73 (MQ=255)
                                                                                                            accGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCa                                                                                                                        <  1:114853/73‑1 (MQ=255)
                                                                                                            accGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTATGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGtc                                           >  2:1332994/1‑150 (MQ=255)
                                                                                                                      tggtggCCAAGCTGACCGCGACGCCTTCGGGGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGc                                 >  1:37817/1‑150 (MQ=255)
                                                                                                                        gtggCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                                       >  2:361736/1‑110 (MQ=255)
                                                                                                                        gtggCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                                       <  1:361736/110‑1 (MQ=255)
                                                                                                                        gtggCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCg                               >  1:625763/1‑150 (MQ=255)
                                                                                                                            ccAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTa                           >  2:908939/1‑150 (MQ=255)
                                                                                                                             cAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTaa                          <  2:886995/150‑1 (MQ=255)
                                                                                                                             cAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTaa                          >  2:609351/1‑150 (MQ=255)
                                                                                                                              aaGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCa                                                                                                                                    >  1:653891/1‑43 (MQ=255)
                                                                                                                              aaGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCa                                                                                                                                    <  2:653891/43‑1 (MQ=255)
                                                                                                                              aaGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAAt                         <  1:1350703/150‑1 (MQ=255)
                                                                                                                                gCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATgg                       <  1:609351/150‑1 (MQ=255)
                                                                                                                                  tGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCa                     >  1:780352/1‑150 (MQ=255)
                                                                                                                                        cgACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTgg               <  1:648138/150‑1 (MQ=255)
                                                                                                                                         accaccTTCGGTCACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGt                                   <  2:939413/125‑1 (MQ=255)
                                                                                                                                         accaccTTCGGTCACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGt                                   >  1:939413/5‑129 (MQ=255)
                                                                                                                                              cTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCaaaa                                                <  2:176816/111‑1 (MQ=255)
                                                                                                                                              cTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCaaaa                                                >  1:176816/1‑111 (MQ=255)
                                                                                                                                               ttCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaac                                                                                        >  2:273495/1‑70 (MQ=255)
                                                                                                                                               ttCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaac                                                                                        <  1:273495/70‑1 (MQ=255)
                                                                                                                                               ttCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTg                >  2:1234831/1‑142 (MQ=255)
                                                                                                                                               ttCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTg                <  1:1234831/142‑1 (MQ=255)
                                                                                                                                                  ggTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCAc     >  1:611717/1‑150 (MQ=255)
                                                                                                                                                     gACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCAACCTGAGCGCCGGCACAACCGTCATCACCGTGCCGGCTAATGCCACCGCGAATACTGCCAcaat  <  2:1124957/150‑4 (MQ=14)
                                                                                                                                                     |                                                                                                                                                     
GACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGT  >  NC_002947/211472‑211770

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